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and  asphalt.  But  for  damage  the  asphalt  exhibited  similar  results  as  on  desert  surface 
(Bee-6  and  Met-12).  Desert  Rock  displacement  and  damage  were  compared  with  the  desert 
surface.  Displacements  were  less  and  damage  was  equivalent  In  the  Desert  Rock  sector 
to  that  along  the  desert  line.  It  was  concluded  that  surface  hardness  contributed  to 
damage.  Desert  Rock  was  undisturbed  desert  surface  and  desert  line  was  disturbed 
producing  a fine,  loose,  sandy  surface. 

Once-ltgrated  (near-ideal)  and  twice- Iterated  (non-ideal)  I values  were  used  to  obtain 
scaled  I versus  range  curves.  The  data  on  these  curves  were  least  squares  fitted  for 
various  Combination  of  events  and  burst  heights.  These  curves  were  divided  Into  three 
categories,  Ideal /near-ideal , light  dust  and  heavy  dust  events.  The  Iq  data  from  asphalt 
surfaces  were  combined  with  I data  from  light  dust  events  and  good  correlation  was 
obtained.  The  data  from  these  curves  were  then  used  to  construct  two  HOB  Charts  - one 
for  light  dust/near-ideal  and  one  for  heavy  dust/near-ideal.  It  was  indicated  that  the 
surface  conditions  for  the  heavy  dust  events  were  unique  and- thus  produced  blast  con- 
ditions which  were  unique.  It  was  suqgested  that  little  Importance  be  given  to  the 
heavy  dust/near-ideal  HOB  Chart  and  perhaps  be  disregarded.  The  light  dust/near-ideal 
HOB  Chart  which  Includes  asphalt  data  Is  more  representative  of  expected  blast  conditions 
over  surfaces  where  precursors  would  be  formed.  It  was  also  felt  that  resultant  I 
values  for  these  events  (light  dust)  were  caused  primarily  by  the  precursor  action).  The 
conclusion  Is  that  the  light  dust/near-ideal  HOB  Chart  Is  the  preferred  chart  for 
general  applications. 
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angstrom 

atmosphere  (normal) 

bar 

barn 

British  thermal  unit 
(thermochemical ) 
calorie  (thermochemical) 
cal  (thermochemical )/cm2 
curie 

degree  (angle) 
degree  Fahrenheit 
electron  volt 
erg 

erg/second 
foot  . 

foot-pound-force 
gallon  (U.S.  liquid) 
inch 
jerk 

joule/kilogram  (J/kg) 
(radiation  dose  absorbed) 
kilotons 
kip  (1000  Ibf) 
kip/inch2  (ksi) 
ktap 

micron 
mi  1 

mile  (international) 
ounce 

pound-force 
pound-force 


1.000  000  X E -10 
1.013  25  X E +2 

1.000  000  X E +2 

1.000  000  X E -28 

1.054  350  X E +3 
4.184  000 
4.184  000  X E -2 
3.700  000  X E +1 
1.745  329  X E -2 
tk  = (t°f  + 459.67)/1.8 
1.602  19  X E -19 

1.000  000  X E -7 

1.000  000  X E -7 
3.048  000  X E -1 
1.355  818 
3.785  412  X E -3 
2.540  000  X E -2 

1.000  000  X E +9 

1.000  000 
4.183 

4.448  222  X E +3 
6.894  757  X E +3 

1.000  000  X E +2 

1.000  000  X E -6 
2.540  000  X E -5 
1.609  344  X E +3 
2.834  952  X E -2 


meters  (m) 
kilo  pascal  (kPa) 
kilo  pascal  (kPa) 
meter2  (m2) 

joule  (J) 
joule  (J) 

mega  joule/m2  (MJ/m2) 

*g1ga  becquerel  (GBq) 

radian  (rad) 

degree  kelvin  (K) 

joule  (J) 

joule  (J) 

watt  (W) 

meter  (m) , 

joule  (J) 

meter3  (m3) 

meter  (m) 

joule  (J) 

Gray  (Gy) 
terajoules 
newton  (N) 
kilo  pascal  (kPa) 
newton-second/m2 
(N-s/m2) 
meter  (m) 
meter  (m) 
meter  (m) 
kilogram  (kg) 
newton  (N) 
newton-meter  (N’m) 


(lbs  aviordupois)  4.448  222 
inch  1.129  848  X E -1 
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pound-force/inch 

1.751 

268  X E +2 

newton/meter  (N/m) 

pound-force/foot2 

4.788 

026  X E -2 

kilo  pascal  (kPa) 

pound-force/inch2  (psi) 

6.894 

757 

kilo  pascal  (kPa) 

pound-mass  (lbm  avoirdupois) 

4.535 

924  X E -1 

kilogram  (kg) 

pound-mass-foot2  (moment 

of  inertial ) 

4.214 

Oil  X E -2 

kilogram-meter2 

(kg’m2) 

pound-mass/foot3 

1.601 

846  X E +1 

kilogram/meter3 

(kg/m3) 

rad  (radiation  dose  absorbed) 

1.000 

000  X E -2 

**Gray  (Gy) 

roentgen 

2.579 

760  X E -4 

coulomb/kilogram 

(C/kg) 

shake 

1.000 

000  X E -8 

second  (s) 

slug 

1.459 

390  X E +1 

kilogram  (kg) 

torr  (mm  Hg,  0°C) 

1.333 

22  X E -1 

kilo  pascal  (kPa) 

*The  becquerel  (Bq)  is  the  SI  unit 

of  radioactivity; 

1 Bq  = 1 event/s. 

**Tha  Gray  (Gy)  is  the  SI  unit  of  absorbed  radiation 
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SECTION  1 

(This  Section  Is  Unclassified) 

INTRODUCTION 

1.1  GENERAL 

In  a previous  effort,  we  evaluated  dynamic  pressure  impulse,  Iq,  for 
Ideal /near- Ideal  blast  waves  by  showing  a correlation  with  wheeled  vehicle 
(1/4  ton  trucks)  displacement  (Reference  1).  We  found  that  there  is  a one- 
to-one  correspondence  between  Iq  and  vehicle  displacement  and  concluded  that 
from  known  vehicle  displacements  the  dynamic  pressure  Impulse  can  be  esti- 
mated with  reasonable  accuracy. 

Here,  we  present  the  results  for  evaluating  Iq,  not  only  for  Ideal/near- 
ideal  blast  waves  but  also  for  Iq  for  non-ideal  blast  waves.  We  combined  the 
data  from  our  first  evaluation  with  1/4  ton  truck  exposures  to  non-ideal 
blast  waves  along  with  that  of  2-1/2  ton  truck  exposures  to  both  non-ideal  and 
ideal /near-ideal  blast  environments. 

Additionally,  we  evaluated  damage  to  wheeled  vehicles  as  a function  of 
vehicle  displacement.  For  a given  displacement,  the  damage  can  then  be  corre- 
lated with  Iq.  The  scaling  of  damage  for  known  values  of  Iq  is  straightforward 
by  simply  using  Sachs  scaling.  The  damage  Itself  as  a function  of  range,  how- 
ever, cannot  be  scaled  by  use  of  a constant  exponent  of  the  weapon  yield  such 
as  W0,4.  The  exponent  of  W will  vary  depending  on  the  burst  height  of  the 
weapon. 

To  keep  this  report  from  being  too  cumbersome  we  do  not  Include  the 
damage  versus  displacement  here.  Instead  another  report  was  prepared  to 
present  this  Information.  We  also  feel  that  In  this  way  the  report  on  Iq  will 
be  of  special  Interest  to  groups  dealing  with  alrblast  phenomena  while  the 
report  on  damage  will  be  of  special  Interest  to  groups  dealing  with  targeting 
and  vulnerability/survivablllty.  This  report  on  Iq  presents  the  procedures 
for  correlating  Iq  versus  displacement,  the  effect  of  surface  conditions  on 
Iq,  the  construction  of  Iq  versus  range  curves  as  well  as  the  construction  of 
HOB  charts  for  iso-Iq. 

1.2  OBJECTIVES 

Our  main  objective  was  the  evaluation  of  dynamic  pressure  impulse  for 
ideal /near-ideal  and  non-ideal  blast  waves  using  wheeled  vehicle  response 
along  with  gage  measurements.  Other  objectives  were  to  correlate  damage 
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with  vehicle  response,  establish  the  Influence  of  diffraction  flow  on  dis- 
placement and  evaluate  the  peak  overpressure  effects  on  damage  to  wheeled 
vehicles. 

1.3  BACKGROUND 

The  exposure  of  wheeled  vehicles  (drag-type  targets)  on  past  nuclear  and 
HE  tests  served  several  purposes  such  as  troop  Indoctrination,  damage  evalua- 
tion and  verification  of  damage  predictions.  Different  types  of  vehicles  were 
exposed  under  a variety  of  conditions.  The  vehicles  included  1/4  ton,  3/4  ton, 
2-1/2  ton,  and  5 ton  trucks.  The  majority  of  vehicles  exposed  were  1/4  ton 
trucks  and  on  several  occasions  they  were  used  as  gages  to  evaluate  the  rela- 
tive effects  of  nuclear  blast  on  drag-type  targets.  The  Information  on 
wheeled  vehicle  exposures  has  been  compiled  and  reevaluated.  A report  has 
been  prepared  which  covers  the  exposure  conditions,  vehicle  damage,  and  dis- 
placements (Reference  2). 

Damage  to  and  displacements  of  drag-type  targets  under  some  nuclear  burst 
conditions  was  far  In  excess  of  damage  and  displacements  for  other  nuclear 
burst  conditions.  This  was  dramatically  demonstrated  on  events  of  Encore-9 
and  Grable-10,  Upshot/Knothole.  The  maximum  damage  of  vehicles  (1/4  ton 
trucks)  exposed  on  Encore-9. was  M.6  (Moderate  II)  and  the  maximum  displace- 
ment was  about  11  metres.  On  Grable-10,  however,  the  1/4  ton  trucks  were 
completely  dismembered  and  displacement  of  parts  (engine,  chassis,  wheels) 
was  greater  than  300  metres.  The  exposure  of  vehicles  on  these  two  events 
was  based  on  pec*  overpressure  values  and  it  was  assumed  that  the  other  blast 
parameters  could  be  calculated  by  using  the  Rankine-Hugoniot  relations.  But 
the  pressure  waveforms  measured  on  Encore-9  (Mach  Region)  were  primarily 
steep-rising  shock  fronts  (near-ideal),  while  the  pressure  waveforms  on 
Grable-10  were  somewhat  distorted  showing  rounded  fronts  with  high  frequency 
oscillations  and  In  some  cases  secondary  peaks  were  observed  (non-ideal). 

There  was  ne  correlation  between  overpressure  and  dynamic  pressure  on  Grable- 
10.  Th1r  difference  in  damage  and  displacement  was  then  attributed  to  the 
precursor  dust-laden  blast  wave  formed  on  Grabie-10.  The  dust  momentum 
Increased  the  total  loads  on  the  vehicles.  cim11ar  results  were  obtained 
on  other  precursor-forming  events.  Ti.us  one  of  the  criteria  used  for  dis- 
tinguishing between  ideal/near-loeal  and  non-ideal  blast  environments  was  the 
waveform  — r.teep  rising  shock  front  for  ideal/near-ideal  and  a rounded  front 
with  secondary  reaks  for  non-ideal.  The  non-ideai  could  be  further  divided 
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into  light  dust  and  heavy  dust  blast  environment.  One  of  the  measures  for 
this  was  that  dynamic  pressures  were  greater  on  heavy  dust  events  than  on 
light  dust  environments. 

When  vehicles  were  used  as  gages  they  were  located  adjacent  to  blast 
line  instrumentation  which  consisted  of  overpressure,  pitot-tube,  differential 
pressure,  total  head,  GREG  and  SNOB  gages.  In  a precursor  dust-laden  blast 
environment  the  dust  registry  coefficient  was  known  for  the  GREG  and  SNOB 
gages  only.  Many  of  the  measurements  for  deriving  dynamic  pressure  versus 
time  In  this  sort  of  an  environment  were  suspect.  A procedure  was  developed 
to  try  to  retrieve  the  suspect  data  (Reference  3).  The  technique  for  retriev- 
ing the  suspect  data  was  based  on  the  comparison  of  measured  data  from  two 
gages  — one  gage  with  a known  registry  coefficient  and  the  other  with  an 
unknown  coefficient.  This  procedure  was  not  successful  since  there  were  a 
limited  number  of  measurements  with  gages  having  a known  registry  coefficient. 
Wheeled  vehicles  on  the  other  hand  responded  by  translational  and  rotational 
motion  to  the  actual  blast  load  environments.  The  displacement  of  the  vehicle 
should,  therefore,  be  a measure  of  the  blast  loads  acting  on  the  target. 

Past  studies  have  resulted  In  calculatlonal  capability  to  compute 
rotational  and  translational  motion  of  drag-type  targets  (References  4 and  5). 
The  technique  of  Reference  4 calculates  either  translational  motion  or  rota- 
tional motion.  The  technique  of  Reference  5 which  Is  a multidegree-of- freedom 
code  calculates  both  translation  and  rotation.  However,  for  both  techniques 
the  calculation  for  rotation  Is  11mit»H  up  to  the  point  of  overturning.  The 
Important  blast  parameter  considered  ' both  techniques  Is  dynamic  pressure 
impulse.  For  fractional  KT  weapons,  diffraction  loading  due  to  overpressure 
influences  the  motion.  These  calculations  provide  an  insight  for  proper 
interpretation  of  the  vehicle  response  to  the  blast  loads. 

In  our  Initial  evaluation,  we  were  concerned  about  establishing  the 
values  of  dynamic  pressure  impulse  for  the  ideal/near-ideal  blast  waves  to 
be  correlated  with  vehicle  displacements.  Two  approaches  were  used  to 
establish  these  values.  One  approach  was  based  on  averaging  the  measured 
overpressure  and  positive  phase  duration  as  a function  of  range  and  then 
combining  this  with  the  Rankine-Hugonlot  relations  together  with  empirically 
derived  equations  to  compute  the  peak  dynamic  pressures  and  the  peak  dynamic 
pressure  Impulse,  respectively.  Another  approach  was  the  use  of  theoretical 
and  analytical  computations  for  static  pressure  impulse.  Ip,  and  Iq  combined 
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with  the  average  values  of  measured  peak  overpressure,  Ps,  and  static  pressure 
impulse.  The  latter  approach,  which  was  the  method  selected,  permitted  the 
construction  of  a set  of  curves  showing  the  ratio  of  Ip/ Iq  versus  Ps  for 
several  burst  heights.  The  Iq  was  determined  using  the  average  measured  values 
of  Ip  and  Ps  on  each  event.  Details  of  this  approach  are  given  In  Appendix  A. 


4 

UNCLASSIFIED 


UNCLASSIFIED 


SECTION  2 

(U)  BLAST  MEASUREMENTS  AND  VEHICLE  EXPOSURES  - NUCLEAR  AND  HE  EVENTS 

2.1  (U)  TECHNICAL  DATA  BASE 

2.1.  (U)  Opera  ions/Events 

(U)  Over  400  vehicles  wen-  ,ja|jos.2u  on  19  nuclear  and  HE  events.  Blast 

measurements  of  one  type  or  anathe"  were  also  obtained  from  these  events.  Of 
the  more  than  400  vehicles  exposed,  only  167  were  used  for  correlation  with 
Iq.  The  vehicles  excluded  for  this  correlation  were  those  located  on  the 
asphalt  and  water  surfaces  and  those  located  In  the  Desert  Rock  sector,  and 
also,  vehicles  which  were  constrained  to  rotate  without  sliding.  Blast  over- 
pressure measurements  were  obtained  on  the  majority  of  these  events.  Only  on 
two  events,  Encore-9  and  Wasp-1,  were  the  dynamic  pressure  measurements  un- 
available. 

(U)  The  listing  of  the  Operation/Events  are  given  In  Table  2.1.  The 

information  provided  includes  the  yield,  actual  burst  height,  scaled  burst 
height,  scaling  factors,  and  ambient  environment.  Ten  of  these  events  are 
additions  to  the  events  used  In  our  first  evaluation.  These  added  events  all 
took  place  at  Nevada  Test  Site;  most  produced  strong  precursor  actions  (non- 
ideal blast  waves).  The  nomenclature  given  for  each  event  refers  to  the  code 
name  and  to  the  shot  number  for  that  Operation. 

2.1.2  (U)  Blast  Wave  and  Vehicle  Displacement  Collations 

(U)  The  nuclear  and  HE  events,  designated  as  Ideal/near-ideal  blast 

waves  were  the  assigned  nine  events  (first  evaluation)  and  the  additional  two 
events:  Wasp-1  and  Wasp  Prime-9.  Our  initial  guidance  for  categorizing  the 

precursor  nuclear  events  was  based  on  the  information  given  In  Reference  6, 
i.e.,  whether  they  were  light  dust-laden  or  heavy  dust-laden  blast  waves. 

The  sorting  of  events  into  one  of  the  three  categories  is  as  follows: 

Ideal/Near-Ideal  Blast  Waves 

Upshot/Knothole:  Encore-9 

Castle:  Koon-3  and  Nectar-6 

Teapot:  Wasp-1  and  Wasp  Prime-9 

Redwing:  Lacrosse-1,  Zuni-3  and  Yuma-4 

HE:  Canadian  20-Ton,  100-Ton  and  Dice  Throw 
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Light  Dust-Laden  Blast  Waves 
Teapot:  Moth-2,  Turk-4,  Bee-6,  and  Apple  1-8 

Heavy  Dust-Laden  Blast  Waves 

Upshot/Knothole:  Grable-10 
Teapot:  Met-12  and  Apple  11-13 
Plumbbob,  Smoky- 15 

Some  of  the  events  within  the  light  dust  and  heavy  dust  categories,  however, 
were  Interchanged  as  will  be  described  later  In  this  report.  The  blast  mea- 
surements of  total  pressure  and  overpressure  as  a function  of  time  were  also 
gathered  and  categorized  accordingly. 

(U)  The  wheeled  vehicle  exposure  data  were  sorted  Into  these  three 

categories.  The  vehicle  data  pertinent  to  this  evaluation  were  the  location, 
orientations,  displacements  and  damage.  Vehicles  were  further  sorted  according 
to  surface  conditions,  l.e. , desert,  water,  asphalt. 

2.2  (U)  ANALYTICAL  APPROACH 

2.2.1  (u)  Blast  Data  Reduction  Procedures 
(U)  We  have  described  the  procedures  used  for  deriving  Iq  for  ideal/ 

near-ideal  blast  conditions  (see  Background  and  Appendix  A).  In  order  to  • 
obtain  the  Initial  best  estimates  of  Iq  for  the  non-ideal  blast  conditions  we 
reduced  the  data  Independently  of  the  previous  data  reduction  procedures.  The 
procedure  we  followed  was  to  Integrate  the  total  pressure-time  and  the  static 
pressure- time  waveforms  to  obtain  the  total  pressure  Impulse  and  the  static 
pressure  Impulse.  A check  was  made  for  certain  intervals  of  the  total  pressure- 
time and  the  static  pressure-time  curves  to  determine  the  range  of  Mach  numbers 
for  those  Intervals.  Where  appropriate  the  Mach  compressibility  corrections 
were  made  to  the  difference  between  the  total  pressure  impulse  and  static 
pressure  Impulse  for  those  Intervals.  The  sum  over  these  intervals  yielded 
the  dynamic  pressure  impulse.  The  details  of  this  procedure  are  presented  In 
Appendix  B. 

(U)  We  made  checks  of  this  procedure  by  comparing  the  Iq  values  we 

derived  with  those  reported  by  BRL  on  Plumbbob,  Smoky-15  (Reference  7)  and  on 
Dice  Throw  (Reference  8).  The  Mach  numbers  ranged  as  high  as  2.2  In  these 
comparisons.  Waveforms  on  Smoky- 15  displayed  a considerable  amount  of  oscil- 
lations while  the  waveforms  selected  from  Dico  Throw  were  reasonably  smooth. 

On  Smoky-15,  at  the  first  stations,  we  assu1  that  the  gage  was  plugged  since 
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the  total  pressure-time  curve  did  not  return  to  the  base  line.  For  computing 
the  total  pressure  impulse  we  terminated  the  integration  of  the  pressure-time 
curve  corresponding  to  the  end  of  the  positive-phase  on  the  static  pressure- 
time curve.  We  assumed  that  the  error  contribution  to  the  total  pressure 
impulse  would  be  small  at  the  tail-end  of  the  waveform.  The  comparisons  of 
Iq  are  shown  in  Table  2.2.  The  maximum  difference  between  these  two  values 
is  about  15  percent. 

2.2.2  (U)  Least  Squares  Fit  Procedures 

(U)  All  of  the  problems  evaluated  where  least  squares  fits  were  applied 

made  use  of  the  TI-59  Prograirmable  Calculator.  The  problems  evaluated  were: 

1)  Dynamic  Pressure  Impulse  vs.  Displacement  [subcases  included 
WWII  and  M38A1  (1/4  ton  trucks)  and  the  M35  and  Ml 35  (2-1/2  ton  trucks)  in 
both  side-on  and  face-on  orientations]. 

2)  Damage  vs.  Displacement  (several  subcases). 

3)  Displacement  vs.  Scaled  Ground  Range  (subcases  included 
side-on  and  face-on  for  1/4  ton  trucks  with  various  height  of  burst  conditions 
including  light  and  heavy  dust). 

4)  Dynamic  Pressure  Impulse  vs.  Scaled  Ground  Range  (subcases 
included  air  burst  and  surface  burst  as  well  as  light  and  heavy  dust  condi- 
tions). 

5)  Scaled  Displacement  vs.  Scaled  Ground  Range  (WWII,  1/4  ton 
trucks,  side-on  and  face-on,  surface  burst  conditions). 

6)  Scaled  Dynamic  Pressure  Impulse  vs.  Scaled  Ground  Range 
(subcases  included  air  burst  and  surface  burst  as  well  as  light  and  heavy 
dust  conditions). 

The  functional  forms  used  for  least  squares  fits  for  these 
preliminary  evaluations  were  as  follows: 

PROBLEM  1 : Iq  = B D"1  or  D = AIp£ 

where:  a and  b are  coefficients 
m and  n are  exponents 
D is  displacement  (metres) 

PROBLEM  2:  da  = aQ  + a-,D  + a2D2,  D < Dm 

da  = 10,  D > Dm 
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where:  da  = damage 

D = displacement 

Cfo  = value  of  displacement  for  which  curve  fit 
(first)  reaches  line  da  = 10 

2 

PROBLEMS  3 through  6:  £n  y = aQ  + a-jX  + a2x 

where:  y = ordinate.  Displacement  or  Iq 

x = independent  variable.  Scaled  Ground 
Range 

(U)  A detailed  description  of  the  least  squares  fit  procedures  is 

given  in  Appendix  C. 

As  stated  previously,  the  above  functional  forms  were  used  for 
the  preliminary  evaluations.  Later,  we  used  other  functional  forms  and  these 
will  be  defined  in  other  Sections  of  this  report. 

2.2.3  (U)  Displacement  - Dynamic  Pressure  Impulse  Least  Squares  Fits 
(U)  In  Appendix  A we  have  described  how  we  obtained  values  of  Iq  for 

the  ideal  shots  used  in  this  report.  In  Appendix  B we  have  described  in  con- 
siderable detail  how  we  obtained  values  of  Iq  from  gage  data  for  the  non-ideal 
shots  of  present  concern.  The  results  are  contained  in  Table  2.3.  Figures 
2.1  through  2.4  show  plots  of  these  values  of  Iq  versus  vehicle  displacement 
for  the  WWII  and  M38A1  one-quarter  ton  trucks  in  side-on  and  face-on  orienta- 
tions. 

(U)  We  obtained  least  squares  fits  for  both  Iq  versus  displacement  and 

for  disnlacement  versus  Iq  in  each  case  following  the  procedures  described  in 
Appendix  C.  The  results  are  given  in  Table  2.4.  Similar  results  are  given  for 
the  REO  M35  2-1/2  ton  trucks  in  Figures  2.5  and  2.6  and  Table  2.4.  (For  the 
REO  M3 5 improved  or  "first  iteration"  values  of  Iq  were  used  as  described  in 
Section  2.2.4. ) 

(U)  We  will  discuss  these  results  in  detail  in  Section  3.1.  Here  we 

note,  however,  that  while  the  data  exhibit  considerable  scatter,  the  curve 
fits  to  the  data  are  well  determined.  Thus,  on  the  average,  Iq  can  be  esti- 
mated with  fair  accuracy  from  a measurement  of  vehicle  displacement  — the 
fitted  curve  serving  as  a calibration  and  the  vehicle  serving  as  a gage.  The 
results*  in  Table  2.4  show  that  values  of  Iq  inferred  from  displacement  mea- 
surements are  of  comparable  accuracy  with  the  gage-reduced  data  for  the  non- 
ideal shots  discussed  in  Appendix  B. 


(U)  *Note  the  listed'  values  of  Ej. 
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2.2.4  (U)  Improved  or  "Iterated"  Values  of  I 

(U)  One  of  the  objectives  of  this  report  is  to  obtain  the  best,  most 

reliable  possible  values  of  dynamic  pressure  impulse  for  non-ideal  nuclear 
shots.  Since  vehicle  displacement  is  independent  of  the  detailed  cause  of 
the  impulse  to  which  the  vehicle  is  subjected,  i.e.,  whether  it  is  caused 
solely  by  airblast  or  by  an  air-dust  mixture,  and  since  displacement- inferred 
values  of  impulse  appear  to  be  comparable  in  accuracy  with  gage-inferred 
values,  it  appears  logical  to  use  displacement-inferred  results  to  improve 
our  estimates  of  dynamic  pressure  impulse.  Therefore,  for  each  instance  in 
this  report  where  we  have  one  or  more  1/4  ton  trucks  exposed  on  a non- ideal 
shot,  side-on  and/or  face-on,  at  the  same  ground  range  as  that  of  a gage  loca- 
tion — in  most  cases  gages  and  vehicles  were  placed  at  the  same  ground  range 
— we  have  averaged  the  values  of  dynamic  pressure  impulse  inferred  from  the 
measured  displacement  by  use  of  the  dynamic  pressure  impulse  versus  displace- 
ment data  fits.  (Note  that  these  fits  are  based  on  both  ideal/near-ideal  and 
non-ideal  nuclear  shot  data,  and  include  some  HE  shot  data  also.)  This  pro- 
vides, then,  a single,  mean  value  of  dynamic  pressure  impulse  inferred  from 
vehicle  displacements. 

(U)  Finally,  since  the  gage- inferred  and  vehicle  displacement-inferred 

values  of  dynamic  pressure  impulse  appear  to  be  of  comparable  accuracy,  we 
have  averaged  the  (average)  gage-inferred  value  with  the  (average)  displace- 
ment-inferred value  to  obtain  a "final"  average  or  "grand"  average  value  of 
dynamic  pressure  impulse.  We  note  that  in  most  cases  the  number  of  gage- 
inferred  values  producing  the  gage  average  result  is  2 or  3 (i.e.,  2 or  3 
of  the  4 possible  values  acceptable);  the  displacement-inferred  result  is 
usually  based  on  the  average  of  2 results,  one  side-on  and  one  face-on. 
However,  the  displacement  results  are  independent  whereas  the  gage-inferred 
results  are  partially  dependent  in  a complicated  way,  as  we  have  noted. 

(U)  The  improved  estimates  of  dynamic  pressure  impulse  obtained  for 

non-ideal  shots  by  this  procedure  are  given  in  Table  2.3.  These  results, 
along  with  our  original  estimates  for  Iq  for  ideal  and  near-ideal  shots,  are 
used  in  Sections  3 and  4 of  this  report. 

(U)  Having  obtained  the  above  results  for  he  non-ideal  shots  as 

described,  we  could  consider  replotting  displacement  versus  Iq  using  the 
improved  Iq  results  with  the  expectation  that  data  scatter  would  be  lessened 
(since  the  average  curve  for  the  originally  plotted  results  is  used  to  improve 
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the  estimates  of  some  of  the  points  now  being  replotted)!  However,  the  dimi- 
nution in  the  data  scatter  would  actually  be  quite  small:  the  scatter  of  the 

ideal /near-ideal  data  points  about  the  curves  is  as  great  as  for  th  non-ideal 
(and  there  are  more  ideal/near-ideal  data  points).  This  may  seem  -'"logical : 
gage  data  should  be  better  for  ideal  shots  since  the  gages  are  not  subject  to 
clogging  by  dust  and  pressure-time  waveforms  are  freer  of  irregularities.  It 
follows  from  this  that  the  ideal  gage  data  should  also  be  better  than  vehicle 
displacement-inferred  results  since  the  latter  were  found  to  be  comparable  in 
accuracy  to  the  gage-inferred  values  of  Iq  for  non-ideal  shots.  For  the  fact 
that  the  ideal  nuclear  shot  results  show  as  much  data  scatter  in  the  displace- 
ment-Iq  plots  as  do  the  non-ideal  we  offer  the  following  explanation: 

1)  there  is  a greater  variety  of  ideal/near-ideal  shots, 
especially  with  regard  to  the  range  of  weapon  yields  in  the  data  base; 

2)  our  calculation  of  four  gage-derived  values  of  Iq,  largely 
independent  from  the  standpoint  of  data  reduction  procedure  (see  Appendix  B), 
and  intercomparison  of  the  calculated  values,  after  culling  obviously  bad 
gage  data,  resulted  in  elimination  of  most  of  the  poor  data  (whether  caused 
by  a gage  result  or  the  data  reduction  process)  before  averaging  to  obtain 
the  accepted  gage-derived  result;  (the'  latter  was  then  averaged  with  the 
average  displacement-derived  result  to  obtain  our  improved  estimate). 

Whether  or  not  the  above  explanation  accounts  entirely  for  it,  the  fact  is 
that  ideal  and  non-ideal  data  have  about  the  same  scatter  on  the  displacement 
versus  Iq  plots.  Thus,  while  initially  it  seemed  reasonable  to  improve  our 
estimates  of  Iq  for  the  non-ideal  shots  by  using  the  vehicle  displacement  data 
— the  "calibration"  of  the  vehicle  as  a gage  being  determined  by  both  the 
ideal/near-ideal  data  along  with  our  initial  estimates  of  non-ideal  data  (gage 
-inferred)  — it  now  appears  equally  reasonable  to  use  the  vehicle  displacement 
data  to  improve  Iq  estimates  for  the  ideal /near-ideal  shots:  the  results  for 

the  latter  are  apparently  not  as  accurate  as  might  have  been  supposed.  The 
same  rationale  justifies  the  procedure:  the  central  curve  in  the  displacement 

-Iq  plots  is  well  determined;  the  data  point  scatter  about  this  curve  is  fairly 
large  but  quite  uniform.  So  while  displacement-inferred  values  of  Iq  have 
considerable  scatter  (there  being  several  uncontrolled  variables  which  affect 
displacement,  (see  Section  3.1),  their  accuracy  is  comparable  with  that  of  the 
gage-inferred  and/or  calculated  results  for  the  ideal/near-ideal  shots. 
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(U)  7ol lowing  the  above  reasoning  we  again  used  displacement- inferred 

values  of  1^  to  improve  our  estimates  in  order  to  obtain  the  best  possible 
estimates  of  Iq  for  construction  of  scaled  HOB  contours.  To  distinguish 
these  results  from  the  previous  results,  we  now  refer  to  the  first  such 
results  (Table  2.3),  applicable  only  to  non-ideal  shots,  as  "first  iteration" 
results.  Our  second  use  of  this  procedure,  which  applies  to  both  ideal  and 
non-ideal  shots  and  which  starts  with  the  best  previous  estimate  (i.e.,  initial 
estimate  for  ideal /near- ideal  shots,  "improved"  or  "first  iteration"  results 
for  non-ideal  shots)  — we  now  refer  to  as  the  "second  iteration".  We  note 
here  the  changes  produced  in  this  iteration  were  quite  modest.  However,  as 
! will  be  seen  in  Section  4.1,  it  is  quite  difficult  to  determine  accr  itely 

the  shapes  of  scaled  HOB  1^  contours  so  that  the  best  possible  data  should 
be  used  for  this  purpose.  Section  4.1  is  the  only  part  of  this  report  in 
we  use  "second  iteration"  values  rf  1^.  The  second  iteration  i.e., 
best  estimates  of  1^  are  given  in  Tables  4.1  through  4.4. 
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TABLE  2.2  (U)  Comparison  of  Dynamic  Pressure  Impulse 

UNCLASSIFIED 


Plumbbob,  Smoky- 15 


RANGE 

Calculated 

BRL  (Report  WT-1407) 

(Metres) 

psi-sec 

kPa-sec 

psi-sec 

kPa-sec 

457 

20.02 

138 

21.99 

152 

789 

3.44 

23.7 

3.85 

26.5 

897 

2.29 

15.8 

2.02 

13.9 

1038 

2.28 

15.7 

2.26 

15.6 

Dice  throw 

- 

Calculated 

BRL  (Report  DNA4377P-1) 

In 

I 

RANGE 

q0 

qo 

(Metres) 

psi-sec  kPa-sec 

psi-sec 

kPa-sec 

134 

1.56  10.79 

1.45 

9.98 

195 

0.59  4.07 

0.53 

3.67 

375 

0.096  0.66 

0.082 

0.57 
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TABLE  2.4  (U)  LEAST  SQUARES  FIT  RESULTS 
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FIGURE  2.1-1  (U)  Dynamic  Pressure  Impulse  versus  Displacement  - 1/4  Ton 

Trucks,  WWII,  Side-on,  Desert  Data  Curves,  Non-Desert 
Data  Points 
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FIGURE  2.2  (U)  Dynamic  Pressure  Impulse  versus  Displacement  - 1/4 
Ton  Truck,  M38A1 , Side-on 
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FIGURE  2.3  (U)  Dynamic  Pressure  Impulse  versus  Displacement  - 1/4  Ton 

Trucks,  WWII,  Face-on 
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SECTION  3 

(U)  RESULTS  AND  DISCUSSION 

3.1  (U)  ANALYSIS  OF  VEHICLE  MOTION 

3.1.1  (U)  Dynamic  Pressure  Impulse  versus  Displacement  Results 

(U)  The  results  for  Dynamic  Pressure  Impulse  versus  Displacement  are 

given  in  Figures  2.1  through  2.6  and  Table  2.4. 

(U)  The  log-log  plots  indicate  a rather  uniform  scatter  of  the  data 

about  the  straight  line  fits.  The  values  of  r,  the  coefficient  of  correlation, 
are  given  in  the  tables.  While  r is  a measure  of  the  goodness  of  fit  of  the 
selected  fitting  function  with  respect  to  the  data  (see  Appendix  C,  "Least 
Squares  Fitting  Procedures"),  it  appears  that  better  values  would  not  be  expec- 
ted with  a different  fitting  function;  the  data  scatter  preclude  that. 

(U)  There  is  more  data  for  the  WWII  than  for  the  M38A1  vehicle.  The 

fact  that  there  is  more  scatter  in  the  WWII  data  stems  from  the  fact  that  a 
largernumbfr  and  variety  of  shots  give  rise  to  the  data  for  the  WWII.  Also 
for  some  or  the  greater  displacements,  large  resultant  errors  are  introduced 
because  of  the  dismemberment  of  the  vehicles,  particularly  for  the  side-on 
orientations.  There  is  some  tendency  for  points  from  given  shots  to  be  above 
or  below  the  fitted  curve.  Since  this  tendency  has  superimposed  upon  it  the 
data  point  scatter  associated  with  each  shot,  the  total  scatter  is  greater 
than  for  the  M38A1  vehicle.  Thus,  for  the  WWII,  side-on  plot.  Figure  2.1,  the 
data  corresponding  to  shots  Koon-3,  Nectar-6,  Canadian  20-Ton  and  100-Ton  have 
somewhat  smaller  displacements  for  a given  value  of  Iq  than  the  curve  fit  while 
data  from  shots  Wasp  Prime-9,  Yuma-4,  Bee-6  and  Apple  1-8  have  somewhat  higher 
displacements.  Similarly  for  the  WWII,  face-on  plot.  Figure  2.3,  (except  that 
here  no  Canadian  20-Ton  and  100-Ton  points  are  shown  — see  next  paragraph). 
Thus,  a small  HOB  effect  is  present  in  the  data.  The  data  points  for  the  M38A1 , 
with  very  limited  shot  variety,  result  in  a rather  uniform  scatter  of  the  data 
for  individual  shots  about  the  curve  fit  as  best  we  can  judge  this  from  the 
small  number  of  data  points. 

3 . 1 . 1 . 1 ( U ) Effect  of 'Diffraction 

In  the  side-on  orientation  both  vehicles  usually  turned 
over,  thus  causing  displacements,  as  measured  — before  and  after  position  of 
vehicle  center  of  mass  — to  be  larger  for  side-on  than  for  face-on  orienta- 
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tlon,  except  possibly  for  very  large  displacements  where  they  appear  to  be 
equal  for  the  WWII  vehicle.  In  the  face-on  orientation  the  vehicle  turned  over 
infrequently  in  the  case  of  the  WWII  and  not  at  all  for  the  M38A1.  A few 
points  corresponding  to  low  yield  shots  are  not  shown  on  the  plots  nor  are 
they  Included  In  the  curve  fits.  These  show  the  effects  of  diffraction.  Data 
for  very  low  yield  shots  {10  tons)  exhibit  much  higher  values  of  uisplacement 
for  small  values  of  dynamic  pressure  Impulse  than  the  plotted  data.  This  Is 
because  the  blast  wave  is  so  short  that  the  vehicle  received  almost  all  of  the 
dynamic  pressure  Impulse  on  Its  struck  side  while  also  subjected  to  a "static" 
overpressure  before  any  significant  counterpressure  has  affected  the  opposite 
side.  In  the  face-on  data  this  effect  Is  evident  also  for  the  Canadian  20-Ton, 
100-Ton  and  Yuma-4  (fractional  KT,  see  Table  2.1)  shots.  The  effect  for  these 
shots  is  not  evident  in  the  side-on  data.  We  believe  the  reason  for  this  is 
that  In  the  side-on  case  vehicle  turnover  caused  comparable  displacements  for 
low  values  of  Iq  for  the  higher  yield  shots,  thus  concealing  the  diffraction 
effect. 

3. 1.1. 2 (U)  Partial  Interpretation  of  Results 

(U)  When  a vehicle  is  struck  by  a blast  wave,  it  will  acquire 

a velocity  proportional  to  the  impulse  received  (If  we  ignore  the  effect  of 
the  motion-resisting  force  during  the  time  of  .delivery  of  the  impulse).  There- 
fore, we  expected  to  find  a good  correlation  between  displacement  and  Iq.  In 
later  discussion  we  will  note  that  damage  correlates  with  displacement  wehreas 
in  a wide  variety  of  studies  involving  blast  damage  to  equipment,  damage  seems 
to  depend  upon  both  pressure  and  Impulse.  We,  therefore,  made  plots  of  dis- 
placement versus  each  of  the  other  blast  field  variables,  namely,  overpressure, 
dynamic  pressure,  and  overpressure  impulse.  As  expected,  we  found  in  each  case 
a very  large  point  scatter  about  a trend  In  which  displacement  Increases  with 
increases  In  any  of  the  blast  field  variables.  We  Interpret  this  as  follows: 
displacement  Is  causally  dependent  upon  Iq;  correlation  between  displacement 
and  the  other  blast  field  variables  (with  very  large  scatter)  is  accidental, 
existing  only  because  these  variables  have  a tendency  (in  the  data  base)  to 
have  high  values  when  Iq  has  a high  value  (and  low  values  when  Iq  has  a low 
value). 

(U)  Upon  analysis  of  the  results  for  the  two  vehicles  and 

two  orientations  we  find  that:  for  fairly  small  values  of  Iq  turnover  occurs 

in  the  side-on  orientation,  but  not  in  the  face-on,  causing  measured  displace- 
ments to  Ke  larger  side-on;  In  the  data  base  there  are  no  large  displacements 
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for  the  M38A1 , especially  face-on,  nor  small  displacements  for  the  WWII, 
especially  side-on;  over  most  of  the  range  of  significant  overlap  of  the  data 
the  disparity  between  side-on  and  face-on  displacements  (for  a given  value  of 
Iq)  Is  about  the  same  for  each  vehicle.  However,  the  WWII  vehicles  were  ex- 
posed to  higher  values  of  Iq  than  the  M38Al's  and  consequently  suffered  larger 
displacements.  For  large  values  of  Iq  (and  hence  displacement)  the  side-on 
and  face-on  displacements  are  close  In  value.  This  Is  because  turnover  of  the 
vehicles  (In  both  orientations)  takes  place,  almost  certainly,  early  during 
the  blast  pulse.  This  tends  to  randomize  vehicle  orientation,  thus  making  the 
blast  force  on  the  vehicle  over  the  remainder  of  the  blast  wave  and  the  effec- 
tive resisting  force  over  almost  all  of  the  trajectory  approximately  the  same 
for  the  two  (Initial)  orientations.  (We  expect  that  this  would  be  true  for 
any  initial  orientation.) 

3.1.2  (u)  Calculations  of  Vehicle  Motion 

3. 1.2.1  (U)  Translational  and  Rotational  Velocities 
(U)  We  made  several  calculations  based  on  vehicle  character- 

istics in  order  to  explain  the  observations  and  differences  between  the 
responses  of  the  two  vehicles.  Most  of  the  vehicle  characteristics  we  require 
are  given  in  Reference  2.  These  are  given  In  the  following  table  along  with 
our  estimates  of  the  area.  A,  presented  to  the  blast  by  each  vehicle  (based 
on  overall  sketches,  not  detailed  drawings  of  the  vehicle),  and  our  calcula- 
tion of  the  amount  of  Inertia  about  the  overturning  axis  for  the  side-on  orien- 
tation, I0T,  (obtained  by  use  of  the  parallel  axis  theorem  and  the  listed  values 
of  the  moment,  I^M,  about  the  center  of  mass). 


VEHICLE/ 

ORIENTA- 

TION 

HEIGHT 
OF  CM 
(cm) 

HALF  WIDTH 
OF  VEHICLE 
(cm) 

MASS  OF 
VEHICLE, M 
(gms) 

rCM  „ 
(qm  cnr) 

‘OX  , * 

(gm  cnr)  (cnr) 

u. 

COEFF.  OF 
FRICTION 

WWII,  SO 

64 

70.2 

1 . 09x1 06 

2, 54x1 09 

1 .25x10  0 2 . 8x1 04 

0.70 

WWII,  FO 

64 

1.09xl06 

1 . 4x1 04 

0.25 

M38A1 , SO 

67.6 

71.6 

1 . 1 9x1 06 

3.00xl09 

1.46xl010  3.0xl04 

0.70 

M38A1 , FO 

67.6 

1 . 1 9x1 06 

1.5xl04 

0.25 

(U) 

The  listed  characteristics  Indicate  that  only  small 

differences  are  to  be  expected  between  the  responses  of  the  two  vehicles.  For 
example,  the  M38A1  has  a slightly  larger  mass  and  overturning  moment  about  the 
center  of  mass  than  the  WWII:  but  It  also  has  a slightly  larger  area  presented 
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to  the  blast  direction  and  a slightly  larger  overturning  moment  arm.  There- 
fore, the  differences  In  the  two  vehicles  tend  to  cancel  one  another  both  for 
translational  and  rotational  response. 

(U)  Using  the  above  data  we  first  performed  a few  very  simple 

calculations  applicable  to  both  vehicles. 

(1)  Taking  Iq  - 80  kPa-sec  * 8 x 10^  dyne  sec/cm^  we  find,  using 

MV0  * IqA  where  M ■ vehicle  mass  and  V0  * initial  velocity 

V0  3 2 x 104  cm/sec  3 450  miles/hour. 

(This  Is  for  a side-on  orientation.)  In  the  experiments  subjection  of  the 
vehicles  to  blast  with  Iq  considerably  less  than  80  kPa-sec  caused  the  vehicles 
to  be  demolished,  as  would  be  expected  with  such  a value  of  the  (calculated) 
Initial  velocity. 

(2)  Taking  Iq  * 10  kPa-sec,  we  obtain 

V0  3 2.5  x 103  cm/sec  3 57  miles/hour. 

If  we  assume  that  most  of  this  Impulse  is  delivered  In  0.1  sec,  the  vehicle 
acceleration  will  be  2.5  x 10^  cm/sec^  and  it  will  travel  a distance  of  1.3 
metres . 

(U)  In  general,  such  calculations  show  that  displacement  of 

the  vehicle  during  the  Intense  part  of  the  blast  wave  Is  fairly  small  compared 
with  the  measured  displacement  of  the  vehicle. 

(3)  The  average  force  acting  on  the  vehicle  during  delivery  of  the 
impulse  in  (2)  is: 

1°5  dyne  sec/cm2  3 x 104  cn,2  . 3 x 1010  d 
0.1  sec 

This  is  much  larger  than  the  force  of  friction,  which  is 

Ffrl-C  3 1.1  x 106  gms  x 980  dynes/gm  x 0.7 
- 7.8  x 108  dynes. 

(4)  If  we  take  Iq  * 1 kPa-sec  delivered  in  0.1  sec,  the  average  force 
on  the  vehicle  is  3 x 109  dynes  which  is  still  considerably  larger  than  the 
opposing  frictional  force.  Thus  even  for  the  smallest  blast  impulses  in  the 
experiments  the  force  exerted  by  the  blast  wave  on  the  vehicle  is  greater 
than  the  opposing  frictional  force.  (For  the  face-on  orientation  the  vehicle 
are?  is  on  bout  one-half  as  large  but  the  coefficient  of  friction  is  only 
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slightly  over  one-third  as  large  as  In  the  side-on  case,  so  the  statement  Is 
true  In  both  vehicle  orientations.) 

(5)  We  made  similar  calculations  for  overturning.  Here  we  assumed  that 
the  effective  moment  arm  for  the  force  on  the  vehicle  due  to  the  blast  Is  the 
height  of  the  center  of  mass.  Then  with  Iq  - 10  kPa-sec,  time  of  Impulse 
delivery  0.1  sec,  we  find  the  angular  acceleration,  a,  to  be 

o 

a 3 150  radians/sec 

After  0.1  sec,  the  angular  velocity,  w.  Is 

to  - 15  rad/sec  3 2.4  revolutions/sec 
and  the  angular  displacement,  6,  Is 

6 - 1/2  at2  « 0.76  radian  3 43°. 

Here  we  have  Ignored  the  torque  owing  to  gravity,  Tg,  tending  to  prevent  over- 
turning. Again  using  the  tabular  data,  we  find 

Tq  3 l.l  x 106  gms  x 980  dynes/gm  x 70  cm 
3 7.5  x 10^®  dyne  cm 

The  torque  applied  by  the  blast  wave  Is  [see  (3)  above]. 

3 x 1010  dynes  x 70  cm  3 2 x 1012  dyne  cm, 
which  Is  much  larger  than  Tg. 

(U)  In  general,  the  angular  displacement  during  the  blast 

pulse  Is  much  smaller  than  the  angular  displacement  well  after  the  pulse  (as 
measured  wherever  the  vehicle  Is  found  overturned)  just  as  was  found  true  for 
linear  displacement;  when  the  bulk  of  the  blast  Impulse  Is  received  over  a 
substantially  larger  time  period,  0.5  seconds,  for  example,  the  linear  and 
angular  displacements  during  the  pulse  time  become  quite  significant. 

(6)  From  the  calculation  in  (5)  we  see  that  for  Iq  * 10  kPa-sec  the 
vehicle  certainly  overturns  (side-on);  for  Iq  * a few  kPa-sec,  It  certainly 
overturns  but  for  Iq  ~ 1 kPa-sec  It  appears  marginal  and  we  cannot  draw  any 
conclusions  from  such  a calculation.  This  is  because  as  the  vehicle  starts 
to  rotate  a larger  area  is  exposed  to  the  blast  wave  and  the  effective  moment 
arm  also  increases.  Also,  the  moment  arm  associated  with  the  force  of  gravity 
decreases  and  when  the  vehicle  has  rotated  approximately  45°,  the  torque  due 
to  gravity  reverses  direction  and  helps  to  overturn  the  vehicle.  With  Iq  = 

1 kPa-sec,  delivery  time  0.1  sec,  the  angular  displacement  at  0.1  sec  is 
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4.3°  and  the  angular  velocity  is  86°  per  second,  with  the  restoring  torque 
owing  to  gravity  rapidly  diminishing. 

(U)  In  the  experiments,  vehicles  exposed  side-on  to  a dynamic 

pressure  impulse  ~ 1 kPa-sec  sometimes  overturned  and  sometimes  did  not. 

(U)  The  above  simple  calculations  are  applicable  to  both  the 

WWII  and  the  M38A1  vehicles,  there  being  very  little  difference  between  the 
two. 

3. 1.2. 2 (U)  Models  of  Vehicle  Motion 
(U)  Since  the  Iq  versus  displacement  data  fits  for  the  two 

vehicles  appear  to  be  somewhat  different  and  only  partially  accounted  for  in 
the  previous  discussion  we  attempted  a further  calculation.  The  forces  we 
should  expect  to  act  on  the  vehicle  during  Its  motion,  other  than  the  blast 
force,  are:  a frictional  force  — during  the  part  of  the  motion  in  which 

the  vehicle  is  sliding  over  {or  gouging)  the  earth;  an  air  drag  force  which 
at  least  In  the  later  stages  of  the  motion  would  resemble  the  usual  type  of 
air  drag  — vehicle  penetrating  and  imparting  some  of  its  momentum  to  the  air 
in  its  path.  Earlier  in  the  motion  when  the  long  tail  of  the  blast  wave 
involves  air  going  past  the  vehicle  this  would  not  seem  to  apply  and  ?s  time 
progresses  the  air  density  in  front  of  the  vehicle  will  vary  considerably  as 
the  air  refills  the  region  partially  emptied  of  air  by  the-blast.  The  Iq 
versus  displacement  plots  contain  data  from  a considerable  variety  of  shots; 
the  air  density  in  front  of  the  vehicle  depends  not  only  on  tl‘  * vehicle's 
(instantaneous)  position  but  also  on  the  weapon  yield,  height  of  burst,  and 
ambient  conditions.  Notwithstanding  these  difficulties  we  attempted  to  write 
a (composite)  "equation  of  motion"  for  a vehicle  and  determine  an  effective 
resistive  force  in  agreement  with  the  displacement  versus  Iq  curve  (for  each 
vehicle,  in  each  orientation). 

(U)  We  assume  the  vehicle  responds  truly  impulsively  to  the 

blast  wave;  that  is,  the  blast  wave  imparts  a velocity  to  the  vehicle  but  the 
vehicle  moves  only  negligibly  during  the  blast  wave.  We  also  assume  that 
during  this  time  the  force  resisting  vehicle  motion  is  negligible.  While 
somewhat  idealized,  these  assumptions  are  in  general  agreement  with  the  simple 
calculations  already  presented.  We  further  assume  that  the  resisting  force 
(after  the  dynamic  pressure  impulse  has  been  received)  is,  on  the  average,  a 
function  only  of  the  vehicle's  instantaneous  velocity,  i.e.,  we  ignore  fluc- 
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tuations  associated  with  the  vehicle's  orientation  (relative  to  both  its  direc- 
tion of  motion  and  the  ground).  This  appears  reasonable  from  one  standpoint: 
the  vehicle's  instantaneous  velocity  Is  seemingly  the  only  continuous  variable 
during  the  motion  to  which  the  effective  resistive  force  can  be  related.  Then 
with  M,  V,  A representing  the  vehicle's  mass,  velocity,  and  Initial  projected 
area  on  the  blast  direction,  we  have 


$ ■ ■ 

where  the  resisting  force  or  equivalently 
Then 


-F(V)  « - 


iTvT  (1) 

function  to  be  determined. 


- /J  g(V)dV  * j where  D - displacement  (at  end  of  motion)  (2) 


But 


D - ClJ 


from  our  curve  fits* 


(3) 


and  MV0  - iqA 


(4) 


* Here  l Is  dimensionless  since 

in  D * In  C + l In  Iq 


so  that 


In  CGS  units,  C has  dimensions: 

[C]  - [cm/  + 1 [gm]^  [sec/  . 


To  obtain  the  values  of  C for  the  two  vehicles/orientations  from  the  fits 
listed  in  Table  2.4,  where  the  fits  are  of  the  form  D * Alq  with  D in  metres 
and  Iq  in  kPa-rec,  we  have 


C = 100A 

(104/ 

since  lm  = 100cm  and  1 kPa-sec  = 104  ^n-e-| - ec 

cnr 
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where  VQ  is  the  vehicle's  velocity  immediately  after  receipt  of  impulse  Iq. 


So 


M 


D - C(^  V0) 


(5) 


Differentiating  the  integral  with  respect  to  its  lower  limit,  V0,  and  substi- 
tuting this  relation  for  D on  the  right  hand  side  yields: 


"A  9{v)dv  m dfc  [/S  a w:  lS]  " m (R} 


(6) 


(7) 


C ,M> 


IV1' 


1 


So  g(V)  = ft  (J) 

and  the  resisting  force  is  given  by  the  right  hand  side  of  the  equation 


dt 


V 

'gWT 


i^ 


(8) 

(9) 


Finally,  = MV^-  = " (pf)  V (10) 

• * » 

(U)  Comparison  of  the  resistive  force,  r.h.s.  of  Equation  (10), 

2 

which  contains  V with  the  values  of  Z given  in  Table  2.4^  we  see  that  the 
velocity  dependence  of  the  effective  resistive  force  varies  greatly  and  that 
2-1,  the  exponent  of  V,  can  be  either  positive  or  negative  to  accord  with  the 
data  fits.  Thus  for  low  values  of  V,  the  resistive  force  approaches  zero  for 
2-1  positive,  infinity  for  2-1  negative,  neither  of  which  is  physically  real- 
istic. This  sort  of  occurrence  can  perhaps  be  expected  with  empirical  data 
fits;  it  does  not  mean  that  the  effective  resistive  force  derived  is  necessarily 
invalid  over  substantial  portions  of  the  vehicle's  motion.  With  the  resistive 
force  derived  we  can  calculate  velocity  and  displacement  of  the  vehicle  as  a 
function  of  time  by  integrating  the  "equation  of  motion".  However,  we  have  no 
physical  Interpretation  for  the  derived  resistive  force  so  we  have  only  used 
it  as  a partial  check  on  physically  motivated  models. 

(U)  The  displacement  versus  dynamic  pressure  impulse  curve  fits 

contain  several  phenomena.  For  small  enough  values  of  Iq  the  vehicles  will 
slide  across  the  ground  opposed  by  a conventional  frictional  force.  This 
occurred  for  both  vehicles  in  the  face-on  orientation  and  for  the  M38A1  side-on. 
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The  WWII  in  side-on  orientation  was  not  exposed  to  blasts  with  such  small  values 
of  Iq.  For  larger  values  of  Iq  vehicle  turnover  occurs  quite  readily  in  the 
side-on  orientation,  less  readily  face-on.  The  M38A1  was  not  exposed  to  any 
blasts  with  sufficiently  high  Iq  for  turnover  face-on.  As  Iq  is  raised  and 
the  vehicle  initial  velocity  Increases  air  drag  can  become  Important,  at  least 
in  the  later  phases  of  the  motion  when  the  vehicle  pushes  air  which  has  moved 
back  In  toward  ground  zero.  At  sufficiently  high  values  of  Iq  the  vehicle  can 
be  partly  or  completely  lifted  off  the  ground,  so  that  there  Is  no  ground  fric- 
tion term.  If  there  is  no  air  drag  term  either,  the  vehicle  Is  in  free  flight, 
probably  tumbling  in  some  fashion,  but  is  eventually  pulled  down  and  strikes 
the  earth.  Thus  sliding,  rolling,  bouncing  across  the  ground  along  with  goug- 
ing of  the  ground  all  can  take  place  at  one  or  another  time  during  the  motion. 
This,  along  with  the  variety  of  weapon  yields,  burst  heights,  ground  ranges, 
and  ambient  conditions,  is  the  reason  that  the  best  that  can  be  hoped  for  is 
a "composite  equation  of  motion". 

(U)  With  these  considerations  in  mind  we  consider  two  cases: 

(1)  low  values  of  Iq  such  that  only  sliding  takes  place;  (2)  higher  values 
where  the  other  mentioned  phenomena  also  take  place  in  which  we  include  a 
constant  frictional  resistance  along  with  a conventional  type  air  drag  resis- 
tance. There  are  clearly  some  portions  of  the  motion  for  which  these  assump- 
tions are  inadequate.  Consider  first  the  case  of  small  dynamic  pressure  impulse 
so  that  only  ground  friction  is  important.  We  take  the  equation  of  motion  to 
be  (it  is  a legitimate  equation  for  this  very  restricted  case): 


€ - «£ 


where  K is  the  constant 


frictional  resisting  force.  Integrating  this  equation  yields 

D = MVq/2K  . ( 

The  coefficient  of  friction,  p,  is  then 

u = K/Mg  . ( 

For  the  M38A1  in  side-on  orientation,  we  have  for  Iq  = 1 , D = 46.6  cm;  with 
values  of  Iq  £ 1 kPa-sec,  the  M38A1  did  not  turn  over. 

(U)  Using  the  tabulated  data  of  vehicle  characteristics  (A 

and  M values)  we  then  obtain 

p = 0.71 
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Such  good  agreement  with  the  measured  value  of  0.7  is  probably  fortuitous. 

(In  fact  over  the  range  of  the  data  0.597  <_  Iq  £ 1.04  kPa-sec  where  vehicle 
turnover  does  not  play  a role,  such  calculations  provide  values  of  y between 
0.68  and  0.83.) 

(U)  For  higher  values  of  Iq  we  include  air  drag  in  a "compo- 

site equation  of  motion"  in  order  to  check  the  result  against  experiment.  Our 
starting  equation  is  then: 

MHT  " MV3x  * + k2^)»  which  upon  integration  and  a (14) 

little  simplification  yields 

D * tn  (1  + ■£-  Vq1" ) . (15) 

Since  MV0  = IqA,  Equation  (4),  we  can  simply  take  two  values  of  Iq  in  the 
desired  range,  calculate  V0  and  find  D from  our  curve  fits,  then  use  Equation 
(15)  to  determine  K]  and  <2-  However,  the  results  are  quite  sensitive  to  the 
choices  of  Iq  values.  We  found,  however,  that  by  fitting  the  data  (for  a given 
vehicle/orientation)  over  a limited  range  to  the  expression  for  D (somewhat 
crudely  — we  have  no  convenient  least  squares  program  for  such  a functional 
dependence)  we  can  obtain  reasonable  results.  [That  this  expression  for 'D  • 
cannot  fit  the  data  over  a large  range  is  clear  since  it  is-  incompatible  with 

t 

the  fit  D = Alq ; the  latter  fit  corresponds  to  the  resistive  force  previously 
derived,  F(V)  = ^ (j^)'*'  V^"^.  Equation  (10).] 

(U)  We  note  that  for  small  V0,  i.e.,  small  Iq,  according  to 

Equation  (15),  v o o 

D *T  v°  = 2*7  v°  • 

the  same  result  achieved  with  no  air  drag.  For  large  V0, 

K 2 

0 - ln  (K7  V0  ) a 2 In  V0 

or  D becomes  proportional  to  In  Iq.  So  Equation  (15)  can  only  be  fitted  to 
very  limited  ranges  of  data. 

(U)  We  fitted  Equation  (15)  to  data  over  a limited  region 

for  a few  cases.  Taking  the  WWII  vehicle,  face-on,  over  the  region  1 £ Iq  <_ 

8 kPa-sec  we  determined  reasonable  values  of  K-|  and  l<2  [such  that  Equation  (15) 
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predicts  displacement  to  within  15%]: 


K-]  = 2.35  x 108 


K2  = 850 


M- ft  = 0.22 


which  is  in  reasonable  agreement  with  the  measured  value,  y = 0.25. 

(U)  Computing  <2V02  we  find  that  for  Iq  = 8,  K2V^,  the  air 

drag  term  in  the  (composite)  model,  is  almost  4 times  as  great  as  the  fric- 
tional term  (initially).  For  Iq  = 1 , ^Vq2  « so  that  here  the  air  drag 
term  is  negligible  throughout  the  motion. 

(u)  We  also  used  Equation  (15)  for  the  M38A1,  face-on,  over 

the  range  1 1 Iq  £ 4 kPa-sec  (about  the  range  of  tl*  data  for  this  case). 
Following  the  same  procedure  we  obtained  a value  for  the  coefficient  of 
friction  y = 0.30  which  is  to  be  compared  with  the  experimental  value  0.25. 

In  these  cases  there  were  no  vehicle  turnovers. 

(U)  Finally  we  used  Equation  (15)  for  the  WWII,  side-on,  over 

the  range  1 < Iq  £ 8.  We  obtained  the  result  y = 0.15. 

(U)  This  is  much  smaller  than  the  measured'  value  of  0.70, 

corresponding  to  the  fact  that  all  of  the  vehicles  turned  over,  almost  all 
180°;  thus  the  vehicles  rolled,  bounced,  and  may  have  partially  or  completely 

lost  ground  contact  for  part  of  the  motion.  Thus  the  measured  coefficient  of 

friction  is  much  greater  than  the  physically  effective  value. 

(U)  Again,  calculating  K2V02  we  find  that  for  Iq  = 1 the  air 

drag  force  is  (initially)  a little  less  than  one-half  the  frictional  force 

while  for  Iq  = 8 the  force  attributed  to  air  drag  in  this  model  is  initially 
27  times  as  large  as  the  frictional  force  — clearly  an  important  effective 
force  over  a significant  portion  of  the  vehicle  motion,  even  if  not  a conven- 
tional air  drag  force. 

(U)  Finally,  each  of  the  "composite  equations  of  motion", 

Equations  (10),  (11),  and  (14),  can  be  integrated  enabling  us  to  obtain 
vehicle  velocity  and  displacement  as  functions  of  time.  (However,  we  have  no 
experimental  data  with  which  to  compare  such  model -predicted  values.)  In 
general,  we  find  that  the  vehicle  motion  takes  place  over  a time  period  - 
one  to  a few  seconds.  Thus  in  most  cases,  all  three  models  predict  that 
vehicle  motion  persists  considerably  longer  than  the  intense  portion  of  the 
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blast  wave,  but  in  some  cases  not  longer,  or  not  much  longer,  than  the  posi- 
tive phase  duration.  Therefore,  in  some  cases  the  air  drag  in  our  model 
corresponds  at  least  roughly  to  a conventional  air  drag  term  (with  varying 
air  density)  while  in  other  cases  it  may  bear  little  resemblance  to  a con- 
ventional air  drag  force. 

(U)  In  sum,  the  various  calculations  described  provide  some 

insight  into  the  meaning  of  the  displacement  versus  Iq  curves.  Further,  they 
•indicate  that  differences  in  results  for  the  two  vehicles  are  not  explained 
to  any  significant  degree  by  differences  in  the  physical  characteristics  of 
the  WWII  and  M38A1  vehicles. 

3. I. 2. 3 (U)  Surface  Effects  and  Explanation  of  Remaining  Discrepancies 
Between  Results  for  the  WWII  and  M38A1  Vehicles 

(U)  Most  of  the  differences  in  results  for  the  two  1/4  ton 

trucks  have  already  been  explained,  the  causes  for  the  differences  being: 

exposure  to  different  ranges  of  I.  values;  relatively  small  regions  of  over- 

M 

lap  of  the  data  for  the  two  vehicles  compared  with  the  total  data  range  for 
each  vehicle;  fairly  large  scatter  of  the  data  points,  especially  for  the  WVill 
vehicle,  this  in  turn  being  caused  by  the  variety  of  shots  which  give  rise  to 
the  data  for  this  vehicle  — there  being  a small,  but  not  completely  negligible 
height  of  burst  effect  present  in  the  data.  We  note  that  the  curve  fits,  good 
straight  lines  on  the  log- log  plots  over  the  range  of  data  for  each  vehicle/ 
orientation,  would  not  necessarily  be  well  represented  by  such  straight  lines 
if  the  data  extended  over  a larger  range  — as  it  does  when  we  combine  the 
data  for  the  two  vehicles,  which  it  now  seems  reasonable  to  do  since  differ- 
ences in  the  data  for  the  two  vehicles  appear  to  result  from  causes  other 
than  differences  in  the  physical  characteristics  of  the  two  vehicles. 

(U)  The  previously  described  calculations  and  data  compari- 

sons led  us  to  look  carefully  at  the  displacement  versus  Iq  data  for  the  two 
vehicles  in  the  region  of  overlap  of  the  respective  data  ranges  (in  each 
orientation)  and  to  associate  the  vehicle  overturn  data  with  the  results 
while  making  the  comparisons.  For  the  two  vehicles  in  face-on  orientation 
no  vehicle  overturns  occur  in  the  region  of  overlap  of  the  data.  In  this 
region  there  are  only  a few  data  points  for  each  vehicle,  10  for  the  WWII 
and  7 for  the  M38A1.  There  is  considerable  scatter.  The  combination  of 
few  data  points  and  large  scatter  causes  considerable  difference  in  least 
squares  fits  to  the  data  for  the  two  vehicles  in  their  common  region. 
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However,  the  data  for  the  two  vehicles  appear  to  be  quite  compatible  in  this 
region;  there  certainly  is  no  large  difference.  (The  overlap  region  only 
spans  1-3  £ Iq  £ 4 kPa-sec.) 

(U)  In  side-on  orientation  there  are  many  more  data  points 

for  each  of  the  vehicles  in  the  region  of  overlap  of  the  data.  The  range  of 
overlap  is  1 £ Iq  £ 8 kPa-sec.  The  slopes  of  curves  for  the  two  vehicles  in 
this  region  are  quite  different.  However,  the  only  significant  difference  in 
the  data  sets  lies  in  the  region  1 £ Iq  £ 3 kPa-sec.  Now  as  we  have  seen  Iq 
~ 1 kPa-sec  is  marginal  for  vehicle  turnover.  The  WWII  vehicle  was  not  exposed 
to  any  blasts  with  Iq  < 1.  The  M38A1  data  contain  5 points  with  Iq  £ 1.04; 
in  4 of  these  cases  the  vehicle  did  not  turn  over.  The  5th  instance,  Iq  = 

1.04,  D = Q.3.u  appears  peculiar  in  that  a 90°  turnover  occurred  — and  a 90° 
turnover  results  in  a measured  displacement  of  about  1.4m.  In  most  cases 
examination  of  our  data  shows,  for  small  and  moderate  displacements,  that  the 
displacement  associated  with  turnover  (90°  or  180°)  accounts  for  a large 
part  of  the  measured  displacement,  the  latter  almost  always  being  somewhat 
larger  than  the  turnover-caused  fraction  of  the  measured  displacement.  (This 
statement  applies  to  side-on  orientation;  there  were  no  turnovers  face-on  with 
only  small  or  moderate  displacement.;  when  a very  large  displacement  occurs,  in 
either  orientation,  the  vehicle  may  turn  over  360°  or  more. which  is  not 
evident  in  our  data  — or  it  may  be  demolished,  which  is  evident.) 

(U)  In  tracing  back  the  difference  in  responses  of  the  two 

vehicles  in  the  region  1 £ Iq  £ 3 kPa-sec,  we  found  first  that  the  M38A1 
turned  over  90°  whereas  when  the  WWII  turned  over  it  was  turned  180°.  There  were  a 
couple  of  instances  with  Iq  - 3 kPa-sec  in  which  the  WWII  vehicle  did  not 
turn  over;  each  occurred  with  a low  yield  shot  so  that,  because  of  the  diffrac- 
tion effect  previously  described,  the  displacement  was  considerably  larger 
than  occurs  for  a similar  Iq  value  associated  with  a higher  yield  shot  (longer 
blast  wave).  A 180°  overturn  causes  twice  the  displacement  as  a 90°  overturn 
and,  since  most  of  the  measured  displacement  is  caused  by  overturning,  the 
WWII  exhibited  larger  displacements  in  the  region  1 £ Iq  £ 3 kPa-sec  than  the 
M38A1  in  all  cases  (owing  either  to  180°  overturn,  as  opposed  to  90°  for  the 
M38A1 , or  to  a low  yield  diffraction  effect).  Upon  further  inquiry  we  found 
that  the  WWII  vehicle  had  been  exposed  to  the  blast  while  standing  on  loose 
soil  while  the  M38A1  stood  on  firm  ground.  We  believe  this  accounts  for  the 
180°  WWII  overturns  compared  with  the  90°  M38A1  overturns  for  the  following 
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reason.  As  the  vehicle  starts  to  slide  (and  rotate  simultaneously)  the  loose 
dirt  will  be  pushed  ahead  of  the  tires.  When  the  vehicle  has  rotated  90°  the 
top  of  the  side  of  the  vehicle  farthest  from  ground  zero  strikes  loose  dirt 
pushing  it  in  somewhat.  Thus  the  vehicle  actually  rotates  somewhat  more  than 
90°  in  a single  overturn,  the  tires  elevated  a little  with  a pile  of  displaced 
loose  dirt  under  them  and  the  top  of  the  vehicle  side  farthest  from  ground 
zero  depressed  a little.  This  makes  it  easier  for  the  vehicle  to  overturn 
again  than  if  it  were  lying  horizontal.  Therefore,  over  a small  range  of  Iq 
just  above  Iq  ~ 1 kPa-sec,  we  should  expect  the  WWII,  placed  on  loose  dirt,  to 
be  likely  to  overturn  180°  while  the  M38A1,  placed  on  firm  soil,  is  likely  to 
overturn  only  90°. 

(U)  This  explanation,  qualitatively,  appears  capable  of 

accounting  for  the  only  as  yet  unexplained  difference  in  the  data  for  the 
two  1/4  ton  trucks.  We  also  note  in  this  connection  that  the  results  for  the 
2-1/2  ton  trucks.  Figures  2.5  and  2.6,  and  Table  2.4  an  very  similar  to  those 
for  the  M38A1  1/4  ton  trucks.  The  2-1/2  ton  trucks  also  stood  on  firm  soil. 

(U)  As  we  noted  in  a previous  section,  the  displacement 

versus  Iq  plots  and  fits  are  based  on  gage-derived  Iq  values  for  both  near- 
ideal and  non-ideal  shots.  As  discussed  elsewhere  in  this  report  we  have  been 
able  to  improve  upon  the  Iq  estimates  for  non-ideal  shots  by  use  of  displace- 
ment data  (the  vehicle  in  essence  playing  the  role  of  a gage,  after  being 
properly  calibrated).  If  we  were  to  use  the  best  Iq  estimates  we  have  obtained 
for  plots  of  displacement  versus  Iq  we  would  find  not  only  a rather  modest 
improvement  but  also  that  the  discrepancy  between  data  for  the  two  vehicles 
is  somewhat  diminished  where  the  data  scatter  is  large.  This  tends  to 
corroborate  the  finding  that  there  is  no  significant  difference  between  the 
two  vehicles  and  that  our  explanation  of  the  differences  In  the  response  data 
is  correct. 

(U)  Finally,  since  we  have  found  that  the  differences  in 

data  for  the  two  vehicles  are  due  to  factors  other  than  differences  in  the 
physical  characteristics  of  the  two  vehicles,  it  Is  meaningful  to  combine 
the  data  for  the  two  1/4  ton  trucks,  in  each  orientation,  and  to  obtain  curve 
fits  for  the  combined  data. 

3. 1.2. 4 (U)  Final  Dynamic  Pressure  Impulse  versus  Displacement 

Least  Squares  Fits 

(U)  We  combined  data  for  the  two  1/4  ton  trucks  and  obtained 
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fits  for  Iq  versus  displacement  and  displacement  versus  Iq  as  shown  in 
Figures  3.1  and  3.2.  When  the  cata  are  combined  the  ranges  of  displacement 
and  Iq  (in  each  orientation)  are  greater  than  for  either  vehicle  alone  and  a 
straight  line  log-log  fit  is  not  as  good  as  a quadratic;  there  is  a signifi- 
cant curvature.  The  statistical  measures  associated  with  the  fits  are  given 
in  Table  3.1. 

(u)  We  likewise  combined  the  data  for  the  2-1/2  ton  trucks, 

REO  M35  and  GMC  Ml 35 . There  is  insufficient  data  for  a side-on  fit  for  the 
GMC  Ml 35  and  no  face-on  data  for  this  vehicle.  Therefore,  we  obtained  data 
fits  for  the  REO  M35,  side-on  and  face-on,  and  a combined  fit  to  the  side-on 
for  the  two  vehicles.  (The  few  data  points  for  the  GMC  Ml 35  do  appear  to 
have  slightly  higher  1^  values  for  a given  displacement  than  Is  the  case  for 
the  REO  M35. ) 

(u)  We  obtained  both  linear  (log-log)  and  quadratic  fits  to 

the  combined  data;  the  latter  is  the  better  and  is  shown  In  Figure  3.3.  The 
results  are  shown  in  Table  3.1.  The  results  for  REO  M35  alone,  side-on  and 
face-on,  are  given  in  Figures  2.5  and  2.6  and  Table  2.4. 

(U)  For  the  combined  displacement-Iq  fits  for  the  two  1/4 

ton  trucks,  and  for  all  of  the  displacement-Iq  fits  for  the  2-1/2  ton  trucks, 
we  used  the  iterated  (improved)  Iq  values  previously  discussed  for  the  non- 
ideal  shots. 

3.2  (U)  DISPLACEMENT  VERSUS  SCALED  GROUND  RANGE.  Results  are  given  in 

Figures  3.4  through  3.11  and  Table  3.2. 

3.2.1  (L!)  Scaling 


(U)  Seal  ing  ground  ranges  used  here  and  elsewhere  in  this  report  are 

given  by 


x = SjR  (Modified  Sachs  scaling) 

(16a) 

or  by 

x = S^R  . 

(16b) 

Here 

P0  1/3  1 

Sd  = (TO77}  wT/3 

(17a) 

and 

Po  1/3  1 
Sd  (T4T)  wn 

(17b) 

where  P0  = ambient  pressure, 

n = a number  used  in  the  assumed  scaling  law,  *nd 
W = weapon  yield  (KT). 
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(U)  We  tried  several  values  of  n (in  addition  to  1/3)  for  each  subcase: 

n = 0.37,  0.4,  and  0.43,  and  in  some  instances,  additional  values  of  n.  In  most 
of  the  subcases  in  this  section  the  optimum  value  of  n was  found  to  be  near  0.4 
(U)  In  general,  the  data  have  a rather  large  scatter  (since  there  are 

many  uncontrolled  variables  involved).  Thus  fairly  large  errors  can  be  expec- 
ted in  inferring  displacement  from  scaled  ground  range.  (See  Table  3.2.) 

(U)  We  also  attempted  to  determine  "scaled  displacement  versus  scaled 

ground  range"  (analagous  to  "scaled  impulse  vs.  scaled  ground  range",  Section 
3.3).  To  do  this  we  proceeded  as  follows.  For  the  scaled  impulse  case  we 
have 

Si Iq  - f(SdR) 

T (°C)  + 273  1/2  u , 2/3 
where  Si  = [— 288 ^ 

(.loci ried  Sachs  Scaling), 

T0  = ambient  temperature 
and  f indicates  a functional  relationship. 

Then  iq=^-f(SdR)  (20) 

and  since  displacement,  D = AIq  (see  Section  3. 1.2. 2 and/or  Appendix  C)  we 
1 1 

have  D = A[^j-  f(SdR)]  . This  can  be  rewritten 

sfD  « A[f(SdR)]*  , (21) 

where  the  left  hand  side  now  represents  scaled  displacement.  In  individual  cases, 
i.e.,  specified  vehic1e  and  orientation,  we  know  A and  l from  our  data  fits  of 
"Displacement  versus  Dynamic  Pressure  Impulse".  From  out  data  fits  of  "Scaled 
Impulse  versus  Scaled  Ground  Range"  we  also  have  a functional  form  f containing 
coefficients  a0,  a-],  and  whose  values  we  have  determined.  If  we  use  these 
results,  scaled  displacement  may  be  regarded  as  known  and  no  fitting  of  data  is 
required.  However,  there  is  considerable  scatter  in  the  data  determining  A 
and  l (i.e.,  the  data  for  "Displacement  versus  Dynamic  Pressure  Impulse")  and 
there  is  also  scatter  in  the  data  for  determining  f(SdR)  (i.e.,  the  data  for 


W 


1 

T73 


(18) 

(19) 
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"Scaled  Dynamic  Impulse  versus  Scaled  Ground  Range").  Therefore,  we  simply 
took  sfo  * G(S<jR)  where  G(SjjR)  is  a function  to  be  determined  by  the  data  and 
It  not  required  that  G(S<jR)  « A[f (S^R) (the  right  hand  side  being  already 
determined  by  the  above  mentioned  data  fits).  This  expresses  "scaled  displace- 
ment" as  a function  of  scaled  ground  range. 

(U)  We  fitted  scaled  displacement,  s|d,  to  the  functional  form  of 

Equation  (2),  Appendix  C,  with  X]  * x ■ scaled  ground  range  and  X2  * x2  (In  two 
cases  where  we  have  sufficient  data  for  such  a fit:  WWII  Vehicle,  Ideal  Sur- 
face Burst  with  Yield  > 1 KT,  side-on  and  face-on  vehicle  orientations).  We 
regarded  l as  fixed.  The  resulting  fits  were  not  quite  as  good  as  the  fits 
obtained  without  scaling  displacements.  It  Is,  of  course,  possible  that  by 
also  varying  l the  results  could  be  Improved.  This,  however,  appears  unwarranted: 
if  the  value  of  l to  be  used  In  an  application  of  these  results  cannot  be  that 
obtained  from  "Displacement  versus  Scaled  Dynamic  Pressure  Impulse"  results, 
scaling  the  displacement  has  no  utility.  Therefore,  we  pursued  this  no  further. 

3.2.2  (U)  Effect  of  Dust 

(U)  Having  obtained  plots  and  data  fits  of  displacement  versus  scaled 

ground  range  for  various  burst  conditions  we  compared  the  results.  (The  group- 
ing of  the  shots  according  to  burst  conditions  will  be  discussed  in  the  next 
section.)  Comparing  Figure  3.5  with  3.4  (surface  burst)  we  see  that  the  effect 
of  dust  over  a limited  span  of  scaled  ground  ranges  is  to  greatly  enhance  the 
displacement.  (Since  n » 0.4  for  both  of  these  figures  scaled  ground  ranges 
are  directly  comparable.  Both  figures  pertain  to  the  WWII  vehicle  side-on.) 

At  very  short  range  displacement  Is  enhanced  by  a factor  of  6.  At  long  range 
the  dust-laden  blast  wave  is  weaker  than  the  ideal  since  energy  is  sapped  from 
the  wave  by  the  dust  and  the  dust  ultimately  falls  behind  the  wave,  no  longer 
contributing  significantly  to  the  dynamic  pressure  impulse. 

(U)  Similarly  comparing  Figure  3.7a  with  3.6  (air  burst)  we  again  see 

that  the  effect  of  dust  is  to  increase  displacement  by  up  to  a factor  of  6 
over  a limited  spatial  region.  In  this  case  n = 0.43  for  both  figures  and 
the  vehicle  involved  is  the  M38A1  in  side-on  orientation. 

(U)  Finally,  comparing  Figure  3.11  (light  dust)  and  3.10a  (near-ideal 

air  burst)  we  obtain  the  same  result  for  dust  enhancement  of  displacement, 
this  time  for  2-1/2  ton  trucks  (with  n = 0.43  for  both  figures). 
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3.3  (U)  DYNAMIC  PRESSURE  IMPULSE  AND  SCALED  DYNAMIC  PRESSURE  IMPULSE 

VERSUS  SCALED  GROUND  RANGE 

(U)  Results  are  given  In  Figures  3.12  through  3.21  and  Tables  3.3  and  3.4. 

3.3.1  (U)  Surface  Burst 

(U)  We  first  discuss  the  surface  burst  data.  Ideal /near- Ideal  blast 

data  spanning  the  yield  range  of  20  tons  - 3.5  MT.  Here,  as  In  the  case  of 
displacement  versus  scaled  ground  range,  we  consider  both  scaled  and  unsealed 
dynamic  pressure  Impulse  and  use  both 

I 

x * S^R  and  x * S^R  . (16a, b) 

In  the  latter  case  we  again  took  values  of  n (In  addition  to  1/3)  of  0.37,  0.4, 
0.43,  and  In  some  cases  additional  values  — somewhat  higher  values  when  data 
over  a very  large  range  of  yields  were  fitted. 

(U)  For  both  scaled  and  unsealed  dynamic  pressure  Impulse  we  obtained 

two  sets  of  fits  to  the  experimental  data,  one  using  only  nuclear  shots  ranging 
in  yield  from  39.5  KT  to  3500  KT  and  one  Including  these  shots  along  with  HE 
shots  extending  'u  y ■ M range  downward  to  20  tons.  We  consider  first  of  all 
the  nuclear  case. 

3 o.l.l  (J)  Nuclear  Events 

(U)  For  both  scaled  and  unsealed  impulse  we  have  two  fits,  one 

extending  over  the  entire  scaled  ground  range  and  one  over  the  large  scaled 
ground  range  domain  of  the  data.  See  Figures  3.12a,  3.12b,  and  3.13.  In 
each  case  the  data  fit  over  the  entire  range  reaches  a minimum  near  the  end 
of  the  data  range  (and  would  thus  be  incapable  of  extrapolation).  Comparing 
the  fits  over  the  entire  data  range  we  see  that  (Tables  3.3  and  3.4)  the  statis- 
tical measures  of  goodness  of  fit  are  a little  better  for  scaled  than  for 
unsealed  impulse  even  with  optimum  value  of  n,  0.425  in  this  case.  [The 
ground  range  scaling  factor,  S^,  depends  upon  n.  See  Equation  (17b).]  The 
results  with  n = 0.4,  however,  are  not  greatly  inferior  to  those  with  n = 0.425. 
Comparing  the  statistical  measures  for  the  scaled  and  unsealed  impulse  fits 
we  have:  scaled,  o ^ j q = 0.1637,  a}  ^ = 0.1611  and  Ej  = 15.11%;  unsealed 
corresponding  values  0.1756,  0.1774,  and  16.64%,  respectively.  We  note  that 
this  is  in  spite  of  the  fact  that  1-r2  is  0.0120  for  the  unsealed  case  and 
0.0173  for  the  scaled  case.  Thus  R^  cannot  be  used  to  determine  which  is  the 
better  of  two  fits  when  one  or  both  variables  are  different  in  one  fit  than  in 
the  other. 
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(U)  Comparing  the  two  fits  over  the  large  scaled  ground  range 

end  of  the  data  range,  again  with  n 3 0.425  determining  S^,  we  see  that  In 
this  case  the  statistical  measures  of  goodness  of  fit  are  a little  better  for 
unsealed  than  for  scaled  Impulse.  (See  Tables  3.3  and  3.4.)  In  each  case 
the  minimum  of  the  fitted  function  here  Is  well  beyond  the  range  of  the  data 
so  that  reasonable  extrapolation  Is  possible. 

(U)  In  obtaining  data  fits  over  the  large  x end  of  the  data 

region  there  Is  some  arbitrariness  concerning  just  how  much  of  the  data  range 
to  Include.  There  Is  In  each  case  substantially  more  data  In  the  low  and 
moderate  x than  In  the  large  x portion  of  the  data  range.  In  each  case  we 
obtained  additional  fits  covering  a greater  or  smaller  range  than  Indicated 
In  Figures  3.12a,  3.12b,  and  3.13.  We  selected  what  we  regarded  as  the  best 
fit  on  the  basis  of  accuracy  of  fit  to  the  data  points  used,  number  of  data 
points  used  and  a subjective  criterion:  that  when  the  fit  obtained  Is  coupled 

with  the  portion  of  the  total  data  fit  covering  the  remainder  (small  x)  part 
of  the  scaled  ground  range  domain  — the  result  Is  a curve  over  the  entire 
range  essentially  coinciding  with  a good  eye-drawn  curve. 

(U)  In  further  analyzing  these  results  we  can  say  that  the 

prediction  capability  afforded  Is  quite  good.  There  is  much  less  data  scatter 
than  in  corresponding  displacement  plots.  This  Is  partly  due  to  the  fact  that 
more  uncontrolled  variables  affect  displacement  than  dynamic  pressure  impulse, 
but  Is  also  partly  due  to  the  fact  that  In  the  data  reduction  process  (by 
many  persons  over  a long  time  period)  there  has  Inevitably  been  some  smooth- 
ing of  data  for  the  blast  field  variables  as  a function  of  ground  range  for 
each  shot.  Also,  some  of  the  latter  results  are  based,  at  least  in  part,  on 
calculations.  (This  Is  not  the  case  for  displacement  values.) 

(U)  Finally,  we  remark  that  In  neither  the  scaled  nor  unsealed 

Impulse  plots  Is  the  scaling  (only  ground  range  is  scaled  in  the  latter  case) 
perfect;  l.e.,  the  data  for  the  4 shots  do  not  coalesce  perfectly.  In  both 
cases  the  Nectar-6  points  lie  a little  above  the  Zunl-3  (and  curve  slopes 
based  on  each  set  separately  would  differ  slightly).  These  points  make  up 
the  upper  end  of  the  curves.  Similarly  lower  down  on  the  curve  the  Koon-3 
points  lie  a little  above  the  Lacrosse-1  points  In  their  (small)  region  of 
overlap. 

(U)  Overall  there  is  not  much  difference  between  the  relia- 

bility of  the  two  sets  of  fits.  Both  are  quite  good;  the  scaled  or  unsealed 
Impulse  results  may  be  used  with  equal  confidence. 
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3. 3. 1.2  (U)  Nuclear  Plus  HE  Events 

(U)  We  now  consider  the  corresponding  sets  of  scaled  and 

unsealed  impulse  data  curve  fits  over  the  entire  yield  range  of  20  tons  to 
3500  KT.  For  the  unsealed  impulse  data  we  found  that  for  the  scaled  ground 
range  x * S^R,  the  optimum  value  of  n (upon  which  depends)  is  about  0.46, 
the  larger  value  of  n in  this  case  being  required  to  coalesce  the  low  yield 
(non-nuclear)  data  on  the  same  curve  as  the  higher  yield  shot  data  — causing 
somewhat  more  separation  of  the  Nectar-6  and  Zuni-3  points  at  the  upper  end  of 
the  curve.  See  Figures  3.14  and  3.15. 

(U)  The  comparisons  for  the  curve  fits  for  scaled  and  unsealed 

impulse  are  very  similar  to  those  just  discussed  over  the  smaller  yield  range: 
the  statistical  measures  are  a little  better  for  the  scaled  impulse  fits  for 
the  whole  data  range  fits  but  a little  better  for  the  unsealed  impulse  curve 
fits  over  the  large  x end  of  the  data.  (See  Tables  3.3  and  3.4.)  Again  the 
whole  data  range  fits  have  a minimum  near  the  end  of  the  data  range;  the 
large  x region  fits  have  minima  well  beyond  the  data  range  and  are  capable  of 
reasonable  extrapolation  beyond  the  data  range.  Once  again,  the  data  for  the 
several  shots  are  not  completely  coalesced  under  the  scalings  used.  In  addi- 
tion to  a moderate  separation  of  the  Nectar-6  and  Zuni-3  points  at  the  upper 
end  of  the  curves,  the  Dice  Throw  (1  KT,  HE)  points  are  somewhat  above  the 
curve  except  at  large  x.  (The  scaled  impulse  curve  contains  one  more  Dice 
Throw  point  than  the  unsealed  total  data  fit.  As  can  be  seen  from  Figure  3.14, 
the  Dice  Throw  point  at  x = 574.7  cannot  reasonably  be  included  in  the  curve 
fit  over  the  total  data  range  without  considerable  worsening  of  the  fit. 

Thus  the  total  data  fit  does  not  extend  to  this  scaled  ground  range,  though 
the  fit  to  the  large  scaled  ground  range  data  includes  this  point.  We  have 
no  reason  to  suppose  this  to  be  a bad  data  point;  the  functional  form 
selected  for  fitting  the  data,  however,  does  not  provide  as  accurate  and 
reliable  a predictive  capability  when  this  point  is  included  in  the  fit  over 
the  entire  data  range.)  At  the  large  x end  of  the  data  the  Koon-3  points  are 
somewhat  higher  than  the  other  data  points.  Deviations  from  the  curves  are 
modest,  however,  and  considering  the  large  yield  range  encompassed  by  the 
data,  the  scalings  used  coalesce  the  data  quite  well.  Again,  good  eye-drawn 
curves  would  coincide  rather  well  with  the  combined  curves  formed  by  the 
large  scaled  ground  region  fits  and  the  non-overlapping  regions  of  the  total 
data  curve  fits.  (At  the  low  ground  range  end  where  the  curves  are  very 
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steep,  eye-drawn  curves  could  Improve  on  the  fitted  curves  somewhat  more  here 
than  In  the  all-nuclear  shot  case  previously  discussed  where  the  yield  range 
covered  is  much  smaller.) 

(U)  Finally,  the  predictive  accuracy  is  somewhat  diminished 

as  compared  with  the  all-nuclear,  smaller  yield  range  data  fits.  (See  Table 
3,3  and  3.4.)  For  example,  the  RMS  errors,  E,  for  the  scaled  impulse  versus 
scaled  ground  range  fits  over  the  whole  data  range  and  over  the  yield  ranges 
(20  tons  - 3500  KT)  and  (39.5  KT  to  3500  KT)  are  19.24%  and  15.11%,  respec- 
tively. The  corresponding  results  for  (unsealed)  impulse  versus  scaled  ground 
range  are:  20.22%  and  16.64%,  respectively.  Similar  results  hold  for  the 
other  statistical  measures. 

(U)  Once  more  from  an  overall  viewpoint  there  is  little 

difference  in  reliability  between  the  two  sets  of  curve  fits  to  the  data. 

Both  are  good.  The  scaled  or  unsealed  impulse  versus  scaled  ground  range 
curve  fits  may  be  used  with  equal  confidence. 

3. 3. 1.3  (U)  Comparison  of  Scaled  and  Unsealed  Dynamic  Pressure 

Impulse  Results 

(U)  To  show  that  this  is  true  we  compared  values  of  Iq 

obtained  by  one  of  the  scaled  and  unsealed  fits  for  a matrix  of  values  of 

yield  and  ground  range.  (We  ignored  effects  of  ambient  conditions  here 

which  are  small  and  have  a negligible  effect  on  the  comparison.)  We  chose 

weapon  yields  W,  of  0.1,  10,  and  1000  KT  and  values  of  R/W^/3  of  150,  200, 

300,  400  and  600  metres.  Here  R is  the  ground  range  and  the  specification 

of  R/W1/3  was  simply  for  convenience  in  the  computation  using  the  scaled  Iq 

versus  scaled  ground  range  fits.  The  actual  domain  of  ground  range  R in  the 

matrix  of  R,  W values  is  69.62  < R < 6000  metres.  We  used  Figures  3.14  and 

3.15  and  the  data  fits  listed  thereon.  For  the  selected  matrix  of  values 

we  thus  calculated  15  values  of  Iq  over  the  very  large  span  of  yields  and 

ground  ranges  using  both  the  scaled  and  unsealed  Iq  data  fits.  (The  range 

of  R extended  into  the  extrapolated  region  of  the  data  fits  - large  x 

region.)  The  average  percent  difference  between  the  two  sets  of  results 

was  13.6%. 

3.3.2  (U)  Non-Surface  Bursts  (HOB  and  Surface  Effects) 

(U)  In  section  3.2.2  we  compared  displacements  for  near-ideal  bursts 

with  those  for  bursts  which  gave  rise  to  appreciable  dust  in  the  blast  wave. 
Prior  to  plotting  and  fitting  this  data,  however,  it  was  necessary  to  group 
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the  shots.  In  order  to  do  this  we  must  nave  some  criterion  since  the  amount 
of  dust  will  vary  with  each  shot  for  each  grc.'nd  range.  Since  a primary 
objective  of  the  present  study  Is  to  obtain  goc-1  estimates  of  Iq,  especially 
for  non-ideal  shots,  we  first  plotted  scaled  Iq  versus  scaled  ground  range  for 
the  various  shots.  We  made  this  choice  because  the  scaled  Iq  plots  involve 
W]/3  fn  scaled  ground  range  — whereas  (unsealed)  Iq  plots  involve  Wn  with  n 
variable,  at  least  over  a limited  range  of  values.  Having  done  this,  however, 
we  find  the  results  obtained  to  be  consistent  for  scaled  and  unsealed  Iq  and 
also  for  displacement  versus  scaled  ground  range.  The  resulting  grouping 
Is  as  follows:  "Light  Dust"  - Grable-10,  Moth-2,  Turk-4,  Bee-6,  Apple  1-8, 

and  Smoky-15;  "Heavy  Dust"  - Met-12  and  Apple  11-13;  "Non-Surface  Near-Ideal"  - 
Encore-9,  Wasp  Prime-9,  Yuma-4  (and  Wasp-1,  although  we  have  no  Iq  data  for 
this  shot  except  what  can  be  derived  from  Iq  versus  displacement;  but  all 
vehicle  displacements  were  very  small  for  Wasp-1,  hence  Inference  of  Iq  is 
very  inaccurate).  This- grouping  is  somewhat  coarse  since  the  data  do  not 
warrant  fine  distinctions  as  to  degree  of  dust.  See,  however.  Section  4. 

(U)  This  grouping  parallels  the  scaled  HOB  grouping  for  the  shots 

— see  Table  2.1,  with  one  exception:  by  virtue  of  its  scaled  HOB  we  would 

expect  Turk-4  to  be  in  the  "Heavy  Du$t"  group  with  Met-12  and  Apple  11-13. 
However,  Turk-4  was  exploded  over  a well-vegetated  region,  thus  creating  less 
dust  than  shots  over  desert  at  the  same  selected  HOB.  As  far  as  dynamic 
pressure  Impulse  is  concerned,  the  dust,  not  the  scaled  HOB,  is  the  relevant 
factor. 

(U)  Figures  3.16a,  3.16b,  and  3.17  show  that  the  data  for  the  various 

shots,  grouped  as  indicated,  fall  into  a reasonably  compact  set.  The  Bee-6 
points  are  a little  high,  perhaps  indicating  somewhat  more  dust  than  other 
shots  in  the  group.  These  figures  are  ror  "Light  Dust".  Figures  3.18  and  3.19 
for  "Heavy  Dust"  exhibit  somewhat  more  data  scatter.  See  Section  4.  (In  the 
case  of  Figure  3.18  the  value  of  n is  almost  irrelevant  since  Met-12  and  Apple 
11-13  have  scaling  factors  which  are  close  in  value;  thus  changing  n produces 
very  little  change  in  the  relative  positions  of  the  data  points  for  the  two 
shots. ) 

(U)  Finally,  we  wish  to  evaluate  the  effect  of  dust  on  Iq.  Comparing 

Figures  3.16a  and  3.12b  (each  of  which  has  n = 0.4)  we  see  that  light  dust 
increases  Iq  by  a factor  of  about  1.2  to  1.6  over  the  region  of  scaled  ground 
range  overlap  of  the  data.  Likewise  comparing  Figures  3.18  and  3.12b  (n  = 0.4 
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for  both)  we  see  that  heavy  dust  enhances  Iq  by  a factor  varying  from  about 
1.5  to  3 over  the  scaled  ground  range  region  common  to  the  two  figures. 

(U)  The  effect  of  dust  on  Iq  can  also  be  seen  by  comparing  the  scaled 

Iq  versus  scaled  ground  range  plots.  Thus  comparing  Figures  3.17  and  3.13 
(or  3.15)  we  see  that  over  the  common  region  of  scaled  ground  range  light  dust 
increases  (scaled)  Iq  by  a factor  varying  from  about  1 to  3.  Again  for  larger 
ground  range  values  (scaled)  Iq  falls  off  much  more  rapidly  when  dust  Is  present 
than  when  It  is  not.  Comparing  Figures  3.19  and  3.13  (or  3.15)  we  see  that 
heavy  dust  enhances  (scaled)  Iq  by  a factor  of  about  3 to  5 over  a limited 
spatial  region;  thereafter  (scaled)  Iq  falls  off  very  rapidly  with  Increasing 
ground  range. 

(U)  Generally  comparing  Figures  3.21  and  3.13  (or  3.15)  we  see  that 

at  large  scaled  ground  ranges  Iq  falls  off  more  slowly  for  a near-ideal  air 
burst  than  for  an  ideal  surface  burst.  For  small  scaled  ground  range  Iq  is 
much  greater  for  a near-ideal  surface  burst  than  for  a near-ideal  air  burst; 
however,  our  data  does  not  extend  to  very  small  ground  ranges.  We  note  also 
that  the  scaled  HOB'S  for  the  shots  we  have  grouped  in  the  near-ideal  air 
burst  category  differ  considerably  from  one  another,  which  probably  explains 
at  least  partly  the  lack  of  good  coalescence  of  the  data  for  the  separate  shots 
in  Figures  3.20a,  3.20b  and  3.21. 

3. 3. 2.1  (U)  Comparison  of  Scaled  and  Unsealed  Dynamic  Pressure 

Impulse  Results 

(U)  We  used  both  the  scaled  and  unsealed  Iq  versus  scaled 

ground  range  plots  (separately)  in  obtaining  preliminary  scaled  Iq  HOB  contours. 
(In  Section  4 we  describe  the  method  we  finally  used  In  obtaining  these  con- 
tours; however,  the  effort  described  here  provides  a good  comparison  of  the  use 
of  the  scaled  and  unsealed  dynamic  pressure  Impulse  results.)  Fxcept  for  the 
surface  bursts  each  shot  provided  data  points  on  the  contours  at  it.  scaled  HOB, 
the  latter  being  different  for  each  non-sirface  shot.  In  using  the  figures 
just  discussed  we  followed  the  fitted  curve  for  the  shot  group  in  general,  but 
allowed  for  deviations  of  the  individual  shots  from  the  group  curve  (since  the 
shots  do  not  coalesce  perfectly  on  the  group  curve).  That  is,  we  assume  that 
there  are  real  physical  deviations  of  individual  shots  from  the  group  curves 
(certainly  there  are  various  degrees  of  dust)  — the  deviations  are  not  due 
only  to  inrccuracy  of  our  knowledge  of  Iq.  In  any  given  case  we  cannot  know 
to  what  extent  our  inaccurate  knowledge  of  Iq  and  to  what  extent  real  physical 
differences  account  for  deviations  from  the  group  curve.  However,  the  group 
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curve  served  as  an  excellent  guide  In  all  cases;  In  many  cases  there  would  not 
be  sufficient  data  to  determine  a reliable  curve  for  a single  shot.  The  re- 
sults obtained  In  this  way,  using  separately  the  scaled  and  unsealed  Iq  versus 
scaled  ground  range  fits,  for  obtaining  values  of  S-jIq  for  the  HOB  contours 
were  very  close  — nearly  always  within  about  2%  of  one  another. 

3.4  (U)  BEST  VALUE  OF  n 

(U)  When  we  started  fitting  displacement  and  Iq  versus  scaled  ground  range, 
the  first  several  cases  we  tried  seemed  to  have  their  optimum  values  of  n In 
the  vicinity  of  0.4.  As  we  have  seen,  comparisons  of  figures  can  be  made 
directly  only  when  the  values  of  n are  the  same.  So  we  adopted  n 3 0.4  as  a 
"standard"  for  comparison  purposes.  For  our  scaled  Iq  versus  scaled  ground 
range  fits,  however,  the  value  of  n is  fixed  at  n » 1/3  (modified  Sachs  scal- 
ing). There  Is  a theoretical  basis  for  this,  although  effect  of  dust  Is  not 
Included;  we  have  seen,  however,  that  It  Is  not  useful  to  scale  displacement. 
Therefore,  having  Improved  our  estimates  of  1^  and  obtained  fits  for  a con- 
siderable variety  of  cases  we  wish  to  establish  our  "optimum"  overall  value  of 
n — one  which  could  be  used  for  cases  not  treated  herein  with  the  expectation 
that  reasonably  accurate  results  would  be  achieved. 

(U)  In  tables  3.2  and  3.3  we  list  the  results  for  many  cases  we  have  fitted. 

The  values  of  n are  optimum  within  approximately  + 0.01,  when  only  one 
value  of  n is  listed.  (In  a few  cases,  as  mentioned  previously,  changing  n 
d'-es  not  shift  the  points  very  much  relative  to  one  another  so  that  the  re- 
sults are  almost  independent  of  n.  This  can  happen  when  only  two  or  three 
shots  supply  data  for  a curve  to  be  fitted*;  generally  when  there  are  three 
or  more  shots,  however,  the  results  vary  quite  significantly  with  n.)  In 
several  cases  we  list  the  results  for  the  optimum  value  of  n and  also  for  n 
k 0.4.  Comparison  then  shows  hew  much  better  the  optimum  Is  than  that  ob- 
tained wlf  n 3 0.4.  A graphical  comparison  is  obtained  by  comparing  Figures 
3.7a,  O.V’b;  3.8a,  3.8b;  3.12a,  3.12b;  3.16a,  3.16b**;  3.20a,  3.2Cb.  The  two 

only  two  shots  contribute  data  it  may  also  happen  that  their  data 
ranges  don't  overlap  in  1 or  scaled  I (or  overlap  is  based  on  one  or  two 
data  points)  so  that  the  data  can  be  "doalesc.ed"  by  a single  curve  within  a 
wide  range  of  n (including  n = i/3). 

(U)**  The  fits  in  Figures  3,16a  and  3.i6b  illustrate  the  case  in  which  one  statis 
ticul  measure  of  goodness  of  fit  is  slightly  better  for  n = 0.4,  the  other 
two  slightly  better  for  n = 0.37.  See  Table  3.3.  The  optimum  value  of  n 
in  this  case  is  propably  between  thes?  t^o  values. 
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figures  in  each  pair  differ  only  in  their  values  of  n.  Inspection  of  the 
entire  set  of  results.  Tables  3.2  and  3.3,  shows  that  the  optimum- value  of 
overall,  is  abcut  n = 0.42  or  0.43.  Most  of  the  cases  for  which  the  data  are 
best  and  there  are  the  most  data  points  cluster  about  such  a value.  In  fact 
all  of  the  "good"  cases  and  almost  all  of  the  "fair"  and  "not  so  good"  ca:  es 
lie  within  the  range  0.40  < n < 0.46.  Perhaps  the  two  best  cases,  3.12a  and 
3.14,  all  nuclear  near-ideal  surface  burst  and  nuclear  + HE  near-ideal  surface 
burst  Iq  versus  scaled  ground  range,  have  optimum  values  of  n of  0.425  a « J 
0.46,  respectively,  the  latter  comparatively  high  value  of  n being  needed  to 
coalesce  the  data  over  the  very  large  yield  range  of  0.04  £ Y < 3500  KT.  As 
would  be  expected  with  many  data  fits,  a goodly  fraction  of  which  have  rather 
few  data  points  with  considerable  scatter,  values  of  optimum  n will  vary  some- 
what more  than  is  physically  realistic.  Figures  3.9,  3.10b,  3.20a  and  3.20b 
may  well  be  examples  of  this. 

(l!)  So  while  we  have  no  criterion  which  allows  us  to  determine  "the"  optimum 
overall  value  of  n,  a value  of  n = 0.42  or  n = 0.43  appears  to  be  a very  good 

choice  and  one  which  can  be  used  with  some  confidence  for  a broader  class  of 

data  than  encompassed  herein.  References  4 and  9 show  similar  ”alues  of  n 
for  blast  damage  to  drag- type  targets. 

3.5  (U)  EFFECTS  OF  NON-DESERT  SURFACE  (WATER,  ASPHALT,  DESERT  ROCK) 

(U)  In  shot  Met-12  gages  and  vehicles  were  placed  not  only  on  desert  surface 
but  also  on  water  (over  desert),  asphalt  and  desert  rock;  on  shot  Bee-6 
similar  data  were  taken  over  asphalt.  We  did  not  use  any  of  this  data  in 
obtaining  the  fitted  curves.  We  now  wish  to  compare  this  data  with  the 

corresponding  desert  surface  results.  We  start  with  Iq  (and  scaled  Iq) 

versus  scaled  ground  range. 

3.5.1  (U)  Effect  on  Dynast 'c  Pressure  Impulse 

(U)  Placing  the  data  pi  +s  for  a water  surface  on  Figure  3.18  and 

3.19  we  see  that  the  Iq  values  for  water  are  substantially  lower  than  the 
desert  values  (except  for  the  lowest  point).  See  Figures  3.18-1  and  3.19-1. 

The  Iq  values  for  asphalt  art  lower  than  those  for  water  over  the  whole  range 
and  much  lower  than  desert,  a factor  of  2 or  more.  Placing  the  asphalt  data 
for  shot  Bee-6  on  Figure  3.16a  again  shows  the  asphalt  Iq  values  to  be  lower 
than  those  for  desert  (3  of  the  4 data  points  lower,  fourth  about  the  same). 

See  Figure  3.16a-l.  On  the  average  the  asphalt  Iq  values  are  lower  by  a factor 
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of  about  1.4.  Use  of  Figure  3.17  with  the  asphalt  data  produces  the  same 
result;  see  Figure  3.17-1. 

3.5.2  (U)  Effect  on  Displacement 

(U)  We  next  compare  displacements  for  the  desert  and  non-desert  data. 

Referring  to  Figure  3.5-1  (side-on)  we  see  that  displacements  for  vehicles  on 
a water  surface  are  about  the  same  as  for  desert  (although  the  Iq  values  for 
water  are  somewhat  lower  than  for  desert).  Use  of  Figure  3.9-1  provides  the 
same  result;  in  this  case  we  compare  face-on  displacements.  Using  the  same 
figures  we  see  that  vehicle  displacement  on  asphalt  is  much  smaller  than  on 
desert  surface.  This  is  consistent  with  the  considerably  lower  Iq  values 
found  over  asphalt. 

(U)  We  also  find  that  vehicles  placed  on  Desert  Rock  are  displaced 

about  the  same  as  on  asphalt  — if  we  assume  no  difference  between  the 
M38A1  and  WWII  vehicles.  This  is  a good  assumption  since  the  displacements 
involved  here  are  large  (see  Section  3. 1.2. 3). 

(U)  In  the  case  of  the  2-1/2  ton  trucks  placing  the  Met-12  (heavy 

dust  shot  over  the  desert  portion  of  the  terrain)  displacements  (Desert  Rock) 
on  the  light  oust  plot  for  2-1/2  ton  trucks.  Figure  3.11,  the  most  comparable 
data  we  have  — we  see  that  the  displacements  over  Desert  Rock  are  uniformly 
substantially  lower  than  over  desert  in  agreement  with  the -above  result  for 
1/4  ton  trucks.  See  Figure  3.11-1. 

3.5.3  (U)  Effect  on  Displacement  versus  Dynamic  Pressure  Impulse 

(U)  We  next  compare  displacement  versus  Iq  for  desert  and  non-desert 

cases.  Comparing  vehicle  displacements  over  water  with  those  over  desert, 
Figures  2.1-1  and  2.3-1,  we  see  that  the  differences  are  rather  small.  There 
is  not  sufficient  data  to  determine  whether  or  not  this  is  significant. 

(U)  Comparing  displacements  over  asphalt  with  those  over  desert,  ag.<in 

using  Figures  2.1-1  (side-on)  and  2.3-1  (face-on),  we  see  that  there  is  no 
difference  in  the  face-on  orientation,  while  the  only  significant  difference 
in  the  side-on  orientation  occurs  for  one  Met-12  point  with  Iq  value  in  the 
vicinity  of  1 kPa-sec,  i.e.,  the  marginal  region  for  turnover. 

(U)  Thus  as  far  as  displacement  is  concerned  for  a specified  value  of 

Iq,  there  is  not  much  difference  among  water,  asphalt  and  desert.  (See,  how- 
ever, Section  3. 1.2. 3 for  effect  of  surface  on  vehicle  turnover  in  cases  where 
Iq  value  is  marginal  for  vehicle  turnover.) 
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3.5.4  (U)  Effect  on  Damage 

(U)  Finally  we  come  to  a comparison  of  damage  between  desert  and 

non-desert  surfaces.  The  damage  data  have  inherently  large  scatter?  Refer- 
ring to  Figure  3.3-1,  we  see  that  the  results  for  water  differ  little  from 
those  for  desert,  though  one  point  appears  low  — this  is  in  face-on  orienta- 
tion. Using  Figure  3.1-1,  side-on,  we  see  that  two  of  the  water  points  lie 
far  below  the  curve  (while  the  other  3 are  on  the  curve  in  the  damage  =10 
region).  For  one  of  the  two  points,  the  vehicle  apparently  was  supported  by 
a dike  retaining  the  water;  in  any  case  it  did  not  turn  over  despite  a large 
displacement  (12  metres).  In  the  second  case,  78  metre  displacement,  the 
vehicle  probably  turned  over  360°  — but  was  not  damaged  as  it  would  have 
been  had  it  been  bouncing  and/or  scraping  the  earth. 

(U)  For  asphalt,  comparing  the  side-on  data  with  that  of  desert. 

Figure  3.1-1 , we  find  that  8 of  the  9 points  are  within  la  of  the  curve 
(mostly  above);  the  ninth  is  a damage  = 10  point  3o  above  the  curve,  yet 
not  very  far  (to  the  low  displacement  side)  from  other  (desert)  damage  = 10 
points.  In  the  face-on  case  comparison.  Figure  3.3-1,  6 of  the  9 points  are 
within  1 a of  the  curve  again  mostly  above.  Two  bf  the  remaining  3 are  damage 
= 10  points,  one  of  these  not  far  from  other  such  points.  The  other  was 
evaluated  as  damage  = 10  owing  to  fire  rather  than  blast  damage;  the  blast 
damage  was  moderate.  The  remaining  point  corresponds  to  a vehicle  for  which 
part  of  the  damage  was  caused  by  debris  rather  than  blast.  Overall,  then, 
damage  for  a given  displacement  appears  to  average  somewhat  higher  on  an 
asphalt  surface  than  on  desert  surface.  (As  we  have  seen,  for  a specified 
displacement,^  is  about  the  same  on  desert  and  asphalt  surfaces,  so  that 
any  dependence  of  damage  on  Iq  over  and  above  that  owing  to  displacement  — 


XU]  * The  damage  data  and  results  are  discussed  in  a separate  report.  Reference 
17;  figures  from  that  report  used  here  follow  Figure  3.21  in  this  report. 
Figures  2.1,  2.2,  2.3,  and  2.8  of  Reference  17  correspond  to  Figures  3.1-1, 
3.2-1,  3.3-1,  and  3.8-1,  respectively,  of  the  present  report.  The  curve 
fit  shown  is,  in  each  case  in  both  reports,  the  damage-displacement  least 
squares  data  fit  for  vehicles  on  a desert  surface.  However,  in  Reference 
17  the  (desert)  data  points  used  in  the  fit  are  shown  whereas  in  the  present 
report  the  data  points  shown  are  non-desert  data  points  (and  a^e  not  used 
in  obtaining  the  fitted  curve). 
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which  is  in  turn  caused  by  the  dynamic  pressure  impulse,  1^  — is  not  a 
factor  here.) 

(U)  Lastly,  comparing  damage  over  a Desert  Rock  surface  with  that 

over  desert  we  have  recourse  to  Figures  3.2-1  and  3.8-1.  From  Figure  3.2-1, 
we  see  that  there  are  not  sufficient  Desert  Rock  data  points  on  which  to 
base  a judgement;  however,  the  Desert  Rock  points  do  exhibit  somewhat  more 
damage  for  a given  displacement  than  the  desert  points.  In  the  case  of 
Figure  3.8-1  (2-1/2  ton  trucks)  the  result  is  more  clear  cut:  there  are 
8 desert  rock  data  points,  none  in  the  damage  = 10  region  and  all  above  both 
curve  fits,  averaging  more  than  2a  above  curve  fit  1 and  1.5a  above  fit  2. 
Thus,  it  is  clear  that  for  a given  displacement  vehicle  damage  is  greater  on 
a Desert  Rock  surface  than  on  desert. 

(U)  Thus  the  pattern  is  consistent:  for  a given  Iq  value,  displace- 

ment is,  at  least  approximately  independent  of  the  surface  on  which  the 
vehicle  is  placed  (provided,  for  the  side-on  orientation  case,  that  I is 
above  the  marginal  region  for  vehicle  turnover,  -1  kPa-sec);  but  for  a given 
displacement  damage  is  greater  on  the  harder  surfaces,  asphalt  and  Desert 
Rock,  than  on  desert  and  greater  on  desert  than  on  water  (over  desert),  i.e., 
a soft,  mushy  surface.  This  is  quite  plausible  since  bouncing,  scraping  and 
gouging  a hard  surface  would  be  expected  to  do  considerable  damage  to  a 
vehicle,  whereas  a soft,  mushy  surface  has  a cushioning  effect. 

3.5.5  (U)  Plumbbob,  Smoky-15  Results 

(U)  The  data  for  Plumbbob,  Smoky-15  differ  moderately  from  other 

desert  data: 

(1)  As  far  as  1^  and  scaled  1^  versus  scaled  ground  range 
are  concerned,  the  data  appear  consistent  with  other  shot  data.  See  Figures 
3.16a,  3.16b,  and  3.17. 

(2)  As  far  as  displacement  versus  1^  is  concerned,  the 
Smoky-15  data  are  within  the  data  scatter.  Smoky-15  data  are  included  in 
these  figures. 

(3)  For  displacement  versus  scaled  ground  range,  the  Smoky-15 
points  are  somewhat  low,  both  side-on  and  face-on;  we  did  not  use  them  in 
the  fits.  (They  are  comparable  with  the  asphalt  points  here.) 
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(4)  For  damage  versus  displacement*  we  did  not  use  the 
Smoky- 15  data  In  the  fits.  On  comparing  the  data  with  the  fits,  however, 
we  found  that  the  Smoky-15  points  averaged  somewhat  below  the  desert  data 
curve  In  side-on  orientation  and  above  the  desert  data  curve  In  face-on 
orientation. 

(U)  There  are  only  a few  points  In  each  Instance  so  that  It  Is 

difficult  to  draw  clear-cut  conclusions.  However,  It  appears  probable  that 
the  (modest)  discrepancies  can  be  accounted  for  by  combination  of  the  usual 
data  scatter  and  a small  Inaccuracy  In  our  estimates  of  dynamic  pressure 
Impulse.  (Our  estimate  of  dynamic  pressure  Impulse  may  be  a little  high.) 


(U)  * Damage  results  are  reported  in  Reference  17. 
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TABLE  3.1  (U)  LEAST  SQUARES  FIT  RESULTS  — 

Combined  Vehicles  - 1/4  Ton  and  2-1/2  Ton  Trucks  (Iterated  I Values) 

CONFIDENTIAL  - FORMERLY  RESTRICTED  DATA 
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(U)  * This  event  and  yield  in  combination  are  CONFIDENTIAL  - FORMERLY  RESTRICTED  DATA. 


TABLE  3.1  (U)  LEAST  SQUARES  FIT  RESULTS  - 

Combined  Vehicles  - 1/4  Ton  and  2-1/2  Ton  Trucks  (Iterated  1^  Values)  (Concluded) 

CONFIDENTIAL  - FORMERLY  RESTRICTED  DATA 
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TABLE  3.2  (U)  LEAST  SQURES  FIT  RESULTS 
Displacement  vs.  Scaled  Ground  Range 
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TABLE  3.3  (U)  LEAST  SQURES  FIT  RESULTS 
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TABLE  3.4  (U)  LEAST  SQUARES  FIT  RESULTS  — 

Scaled  Dynamic  Pressjre  Impulse  vs.  Scaled  Ground  Range  (Concluded) 
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UNCLASSIFIED 


VEHICLE 

TYPE 


WWII 

WHII 

WWII 

WWII 


WWII 

WWII 

WWII 

WWII 


M38A1 


WWII 

wwi: 

WWII 

WWII 

WWII 

WWII 

WWII 

WWII 

WWII 

WWII 


M38A1 

M38A1 

M38A1 

M38A1 

REO  M35 
REO  M35 
REO  M3  5 
REO  M35 
RED  H35 

GMC  MI  35 
GMC  Ml 35 
GMC  MI  35 


TABLE  3.5  (U)  Non-Desert  Surface  Data 
UNCLASSIFIED 


GROUND 

Iq 

Displacement 

RANGE  (M) 

(kPa-sec) 

IMJ 

da  SURFACE 

TEAPOT,  Bee-6,  Side-on 


549 

6.59 

20.1 

7 

Asphalt 

610 

4.16 

5.7 

3 

Asphalt 

701 

1.60 

3.6 

3 

Asphalt 

777 

1.02 

1.7 

3 

Asphalt 

TEAPOT,  Bee-6, 

Faca-on 

549 

6.59 

13.4 

3 

Asphalt 

610 

4.16 

1.6 

7 

Asphalt 

701 

1.60 

1.1 

2 

Asphalt 

777 

1.02 

0.5 

3 

Asphalt 

TEAPOT,  Apple  1-8,  Side-on 

869 

— 

2.7 

3 

Desert  Rock 

TEAPOT,  Met- 12 

, Side-on 

610 

28.1 

113  • 

10 

Water 

686 

19.6 

103 

10 

Water 

762 

13.1 

■157 

10 

Water 

838 

, 7.55 

78 

5 

Water 

914 

3.94 

12 

1 

Water 

610 

19.0 

68 

10 

Asphalt 

686 

13.7 

59 

10 

Asphalt 

762 

8.4 

23 

10 

Asphalt 

838 

5.58 

14 

6 

Asphalt 

914 

1.72 

1.0 

1 

Asphalt 

TEAPOT,  Met- 12 

, Side-on 
• 

610 

.. 

81.0 

10 

Desert  Rock 

686 

54 

10 

Desert  Rock 

762 

_ _ 

22 

9 

Desert  Rock 

838 

— 

11.8 

6 

Desert  Rock 

762 

„ 

15.3 

8 

Desert  Rock 

762 

_ — 

8.1 

9 

Desert  Rock 

762 

_ — 

16.9 

7 

Desert  Rock 

838 

_ _ 

10.3 

6 

Desert  Rock 

914 

— 

5.7 

6 

Desert  Rock 

686 

16.7 

9 

Desert  Rock 

762 

_ _ 

13.1 

8 

Desert  Rock 

838 

— 

6.7 

5 

Desert  Rock 
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TABLE  3.5  (U)  Non-Desert  Surface  Data  (Concluded) 
UNCLASSIFIED 


VEHICLE 

GROUND 

Iq 

Displacement 

TYPE 

RANGE  (Ml 

(kPa-sec) 

(M) 

da 

SURFACE 

TEAPOT,  Met- 

12,  Face-on 

UWII 

610 

28.1 

no 

10 

Water 

WWII 

686 

19.6 

91.4 

10 

Water 

UWII 

762 

13.1 

88 

10 

Water 

WWII 

838 

7.55 

8.8 

1 

Water 

WWII 

914 

3.94 

3.3 

1 

Water 

WWII 

610 

19.0 

71.3. 

.9 

Asphalt 

WWII 

686 

13.7 

41.4 

10 

Asphalt 

UWII 

762 

' 8.4 

20 

10 

Asphalt 

WWII 

838 

5.58 

4.C 

3 

Asphalt 

WWII 

914 

1.72 

0.5 

3 

Asphalt 
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FIGURE  3.2  (U)  Dynamic  Pressure  Impulse  versus  Displacement  - WWII  and 

M38A1 , Face-on 


Dynamic  Pressure  Impulse  versus  Displacement  - 2-1/2  Ton 
Trucks,  GMC  Ml 35  and  REO  M35,  Side-on 


UNCLASSIFIED 


FIT:  tn  D • 9.423  - 0.0195x  - 0.0000238  x2  DESERT  DATA 
• 0.925  T..  ^„  . ..... . - CURVE  FIT 


• Met-12,  Water 
9 Met-12,  Asphalt 
b Bee-6.  Asphalt 

Met-12,  Desert  Rock,  M38A1 
■ Apple  1-8,  Desert  Rock,  M38A1 


Scaled  Ground  Range  (Metres) 


UNCLASSIFIED 

FIGURE  3.5-1  ■( U)  Displacement  versus  Scaled  Ground  Range  - 1/4  Ton  Trucks, 
WWII,  Side-on,  Light  and  Heavy  Dust,  n = 0.4,  Desert  Data 
Curve,  Non-Desert  Data  Points 


UNCLASSIFIED 


K FIT:  In  D 

i •; 

» 13.07  - 0.0562X  + 0.0000473  x2  J 
R2  - 0.993 

TT-^  - 1 i.  : . : r 

. : . r : : -3Tr.-r.  . : 

L r : 

-.r.il-l-i': 

- - * • - ♦— * 

i ri  • Moth-2 

£ — * 

: : : : : : : rri  to  Turk-4 

r - • - • } » - • • r 

i— . - - . , . . . - . 

t :j  ♦ Grable-10 

i rrj.i'.fri.’r!  4rrrr-rT>  r t T?r:-?  ■— ~i 

in 


« 


ls.?_a3J 


_ Scaled  Ground  Range  (Metres) 

200  300  400  '500' 


UNCLASSIFIED 


FIGURE  3.7a  (U)  Displacement  versus  Scaled  Ground  Range  - 1/4  Ton  Trucks 

M38A1,  Side-on,  Light  Dust,  n = 0.43 

7? 


UNCLASSIFIED 


UNCLASSIFIED 


FIT:  U 0 • 9.S71  - O.OSttx  ♦ 0.0000773  x2 
ft2  - 0.176 


♦ Bmk3 
■ £mH-3 


jj-1  Sealed  Ground  Range  (Metres) 

0 TOO  200  300'  400 

UNCLASSIFIED 

FIGURE  3.8a  (U)  Displacement  versus  Scaled  Ground  Range  - 1/4  Ton  Trucks 
WWII,  Face-on,  Near-Ideal  Surface  Burst,  n = 0.43 
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DISPLACEMENT 


UNCLASSIFIED 


*pri  FIT:  Zri  D « 9.397  - 0.0S12x  + 0.0000587  x2  hfti'  l' 
•111  R2  • 0.7962 


lO^EEERET 


'i- ; E.-.i 


_i  o Encore-9,  REO  M35 
i Q Wasp-1 , REO  M35 
♦ Wasp-1,  GMC  Ml  35 
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10" 2 i ' Scaled  Ground  Range  (Metres) 

0 100  200  300  400  500  600 

UNCLASSIFIED 

FIGURE  3.10a  (U)  Displacement  versus  Scaled  Ground  Range  - 2-1/2  Ton  Trucks 
REO  M35  and  GMC  M135,  Side-on,  Near-Ideal  Air  Burst,  n = 0 
78 

UNCLASSIFIED 


(kPa-sec) 


UNCLASSIFIED 


-FIT  1:  in  Iq  - 6.466  - 0.0287x  + 0.0000280  x2 

R2  - 0.985 

FIT  2:  in  Iq  ■ 3.755  - 0.0135x  + 0.00000892  x2 

R2  - 0.968 


® Nectar-6 
■ Lacrosse-1 
-*■  Koon-3 
Q Zunl-3 


- \\ 


Scaled  Ground  Range  (Metres) 

300  400  500 ' 
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3.12b  (U)  Dynamic  Pressure  Impulse  versus  Scaled  Ground  Range  - Near-Ideal 

Surface  Burst,  n = 0.4 
83 
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Scaled  Ground  Range  (Metres) 


' UNCLASSIFIED 

FIGURE  3.13  (U)  Scaled  Dynamic  Pressure  Impulse  versus  Scaled  Ground  Range 

Near-Ideal  Surface  Burst 
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FIGURE  3.16a  (U)  Dynamic  Pressure  Impul se  -versus  Scaled  Ground  Range 

Light  Dust,  n = 0.4 
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HFIT:  in  Iq  « 5.732  - 0.0157x  - 0.000000764  x2  H ! i I 

R2  - 0.941 


G 

Bee-6 

a 

Moth-2 

+ 

Turk-4 

• 

Apple  1-8 

■ 

Smoky- 15 

«r 

Qrable-10 

■ • ■ Scaled  Ground  Range  (Metres)  ■ 

100  200  300  400  500 

UNCLASSIFIED 

FIGURE  3.16b  (U)  Dynamic  Pressure  Impulse  versus  Scaled  Ground  Range 
Light  Dust,  n = 0.37 
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FIT:  in  Iq  - 5.732  - 0.01 57x  - 0.000000764  x2 

R2  - 0.941  DESERT  DATA  CURVE  FIT 


KJ 


10  ' Scaled  Ground  Range  (Metres) 

'o  100  200  300  400  500  ' 600  700 

UNCLASSIFIED 

FIGURE  3. 1 6b-1  (U)  Dynamic  Pressure  Impulse  versus  Scaled  Ground  Range  - Light 

Dust,  n = 0.37,  Desert  Data  Curve,  Non-Desert  Data  Points, 
Bee-6,  Asphalt 
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Iq  (kPa-sec) 


'li 


UNCLASSIFIED 


FIT:  £«  Iq  ■ 4.244  - 0.00479x  - 0.0000552  x2  4. 
mi  ■ 0.909  DESERT  DATA  CURVE  FIT  \ 


rf 


iti  • Met-12,  Water 
e Met-12,  Asphalt 


Scaled  Ground  Range  (Metres) 

) 100  200  300  400  ' "500  "500"  ‘ ' 700 

UNCLASSIFIED 

FIGURE  3.18-1  (u)  Dynamic  Pressure  Impulse  versus  Scaled  Ground  Range  - 
Heavy  Dust,  n = 0.4,  Desert  Data  Curve,  Non-Desert  Data 
Points 
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FIGURE  3.19  (U)  Scaled  Dynamic  Pressure  Impulse  versus  Scaled  Ground  ^anqe 
Heavy  Dust 
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FIGURE  ,1.20b  (U)  Dynamic  Pressure  Impulse  versus  Scaled  Ground  Range 
Near-Ideal  Air  Burst,  n = 0.50 
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A'. 
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FIGURE  3.1-1  (U)  Damage  versus  Displacement  - 1/4  Ton  Trucks,  WWII,  Side-on 

(Reference  17)  Desert  Data  Curve,  Non-Desert  Data  Points  (Yield  > 1 KT) 
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(U)  SECTION  4 

(U)  HEIGHT  OF  BURST  CHARTS  FOR  Iq 

4.1  (U)  SCALED  DYNAMIC  PRESSURE  IMPULSE  VERSUS  SCALED  GROUND  RANGE  CURVES 

(U)  Prior  to  constructing  the  Height  of  Burst  (HOB)  Charts,  we  again 

went  through  an  iterative  procedure  to  obtain  Iq  values  as  a function  of  range 
for  each  nuclear  event  (HE  information  was  not  used  in  constructing  the  HOB 
Charts).  The  procedure  used  was  similar  to  that  described  in  Section  2.2.4, 
except  that  we  used  the  quadratic  least  squares  fits  (see  Section  3. 1.2. 4) 
for  the  combined  WWII  and  M38A1  (1/4  ton  trucks)  and  for  the  combined  M35  and 
Ml 35  (2-1/2  ton  trucks)  to  obtain  In  . The  data  for  each  nuclear  event 
are  shown  in  Tables  4.1  through  4.3  (Iq  derived  from  Desert  Rock  or  from 
water  surfaces  were  not  used  in  these  Tables).  Tables  4.1  and  4.2  show  the 
data  for  ideal /near- ideal  blast  environments  for  surface  bursts  (4  events) 
and  for  air  burst  (4  events),  respectively.  Table  4.3  shows  the  data  for 
non- ideal  blast  waves  (8  events).  Iqg  in  these  Tables  represents  the  original 
estimates  (non-iterated)  for  the  ideal/near-ideal  events  and  the  first  itera- 
tion for  the  non-ideal  events.  Iqavg  represents  the  calculated  values  from 
vehicle  displacements  (both  1/4  ton  and  the  2-1/2  ton  trucks).  The  Iq0  and 
Iqavg  were  then  averaged  to  produce  Iq  (actual).  Note  that  in  some  cases 
there  were  no  derived  Iqavg  values.  In  these  cases,  vehicles  were  not  loca- 
ted at  these  ranges  so  only  IqQ  was  used  and  these  were  then  labeled  as  Iq. 

The  Iq  along  with  range  were  scaled  to  1 KT  sea  level  conditions  using  modi- 
fied Sachs  scaling.  The  scaled  Iq  versus  scaled  range  data  for  each  event 
were  grouped  according  to  their  scaled  height  of  burst.  Our  procedures 
described  in  this  Section  on  grouping  and  applying  least  squares  fits  to  the 
data  for  various  burst  conditions  are  similar  but  a refinement  to  that 
described  in  Section  3.3. 

(U)  The  ideal /near- ideal  surface  bursts  were  comprised  of  four 

events:  Koon-3,  Nectar-6,  Lacrosse-1,  and  Zuni-3.  Least  squares  fits  were 

made  to  these  four  events.  Of  the  four  ideal/near-ideal  air  bursts,  there 
is  some  uncertainty  about  Yuma-4,  i.e.,  whether  it  was  ideal/near-ideal  or 
non-ideal.  We  previously  considered  this  to  be  an  ideal /near-ideal  blast 
event  since  the  waveforms  displayed  steep-rising  shock  fronts  (one  of  the 
measures  for  identifying  blast  wave  categories).  However,  in  our  final 
resolution  of  Iq  for  Yuma-4  these  values  are  all  above  our  initial  estimates 
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(see  Appendix  A),  implying  that  these  higher  values  of  IQ  could  Le  the  result 
of  precursor  action.  But  Yuma-4  was  a fractional  KT  device  and  infraction 
could  contribute  to  the  motion.  The  derived  Iqavg  from  vehicle  ^placements 
would  then  indicate  higher  values.  Data  on  Yuma-4  is  insufficient"  to  be  able 
to  distinguish  the  discrepancy  between  diffraction  effects  and  precursor 
action.  Another  questionable  event  is  Wasp-1.  Only  peak  pressure  values 
were  obtained  on  this  event  so  the  IqQ  values  for  Wasp-1  were  obtained  by 
scaling  Encore-9  to  the  actual  Wasp-1  conditions.  The  scaling  using  modified 
Sachs  is  straightforward  and  there  is  no  reason  to  doubt  the  validity  of  the 
IqQ  values  on  Wasp-1.  However,  there  are  two  factors  which  contribute  to  the 
questionable  data  on  Wasp-1.  First,  the  actual  ground  zero  was  somewhat  dis- 
tant from  the  intended  ground  zero  causing  the  orientation  of  the  vehicles  to 
be  at  some  angle  other  than  side-on  or  face-on.  Our  Iq  versus  displacement 
correlations  are  either  for  side-on  or  face-on  orientations.  The  second 
factor  was  that  the  displacements  were  fractions  of  a metre.  Both  factors 
would  contribute  to  large  errors  in  computing  Iqavg  from  displacements. 
Nevertheless,  least  squares  fits  were  made  where  we  grouped  Wasp  Prime-9  with 
Yuma-4  and  Encore-9  with  Wasp-1. 

(U)  There  were  eight  non-ideal  air  bursts  within  the  scaled  HOB  from 

about  40  to  73  metres.  Various  combinations  of  these  events  were  used  for 
obtaining  least  squares  fits.  We  combined  light  dust  events  with  heavy  dust 
events  in  an  attempt  to  evaluate  any  significant  differences  between  these  two 
categories.  After  establishing  which  of  the  events  were  light  dust  or  heavy 
dust  blast  waves  we  used  the  Ig  values  derived  from  the  asphalt  surfaces  for 
comparisons  (see  Table  4.4)  with  the  light  dust  events. 

(U)  A total  of  23  least  squares  fits  were  made  to  the  scaled  Iq  versus 

scaled  ground  range  data.  Table  4.5  lists  the  combination  of  events,  the 
burst  conditions,  the  number  of  data  points  used  in  the  fit,  the  derived 
functional  least  squares  equation,  the  associated  statistical  data,  and 
comments  regarding  what  was  used  in  the  fits  and  the  scaled  range  applicable 
to  the  fit.  Under  the  remarks  column  we  make  a distinction  between  Iq  and 
Iq  measured.  The  Iq  is,  of  course,  the  average  of  Iq  and  Iqavg  and  Iq 
measured  are  the  measured  values  only.  The  Iq  measured  were  used  as  addi- 
tional points  for  the  fits  as  well  as  to  establish  the  validity  of  the 
measurements.  In  one  case  we  fitted  the  data  using  only  Iq  measured  (Met-12 
and  Apple  11-13).  The  values  for  these  events  are  given  in  Table  4.6.  Also, 
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in  a couple  of  cases  we  fitted  the  curves  by  eye  because  of  lack  of  data 
points,  e.g.,  Bee-6  and  Wasp  Prime-9.  The  total  number  of  Figures  prepared 
for  the  various  combinations  of  events  was  18  and  the  results  are  shown  in  Figures 
4.1  through  4.18.  These  curves  are  shown  in  sequence  according  to  the  average 
burst  height. 

(U)  Two  least  squares  fits  were  made  for  the  surface  bursts  (Figure  4.1). 

The  reason  for  this  is  that  the  functional  form  selected  did  not  fit  the  data 
points  at  the  larger  ranges.  Fit  1 and  Fit  2 are  of  the  same  functional  form 
but  Fit  1 contains  all  data  points  used  and  Fit  2 contains  10  of  these  points 
which  were  arbitrarily  selected  to  fit  the  data  points  at  the  larger  ranges. 

We  place  high  confidence  on  these  fits  because  of  the  high  coefficient  of 
determination,  R2,  (see  Table  4.5)  and  the  low  values  for  the  $ ^tistical 
parameters  and  in  particular  the  percentage  errors,  Ej(%),  which  are  about 
13  percent.  Again,  these  scaled  Iq  values  represent  the  ideal/near-ideal 
blast  waves. 

(U)  The  next  set  of  curves.  Figures  4.2  through  4.6  rep>  ent  the 

scaled  Iq  for  non- ideal  blast  waves.  The  a/erage  range  for  the  HOB  is  from 
.41.5  metres  out  to  44.6  metres  — an  arbitrary  selection  for  the  HOB.  Note 
on  Figure  4.2  that  the  asphalt  data  on  Met-12  fits  reasonably  with  the  data 
of  Turk-4,  which  we  later  defined  to  be  light  dust.  The  asphalt  data  includes 
the  Iqavg  calculated  from  vehicle  displacements  over  asphalt  surface  using 
the  quadratic  fits  to  the  Iq  versus  displacement  which  were  obtained  for  the 
desert  surface  (Nevada  Test  Site)  and  for  the  coral  surface  (Pacific), 

Section  (3. 1.2. 4).  The  Iq  versus  displacement  for  the  asphalt  surface  are 
comparable  to  the  desert  and  coral  surfaces  (see  Section  3.5.3). 

(U)  Figures  4.3  and  4.4  combine  what  we  consider  to  be  light  dust  and 

heavy  dust-laden  blast  waves.  Met-12  Iq  values  are  generally  higher  than  the 
Turk-4  or  Apple  11-13  Iq  values  (Figure  4.3).  The  Turk-4  data  and  Apple  Il- 
ls data  appear  to  coalesce.  Matter  of  fact,  when  we  fit  just  the  Iq  values 
on  Turk-4  and  Apple  11-13  we  get  good  agreement  (Figure  4.4,  Fit  2).  This 
is  not  the  case  when  we  included  the  Iq  measured  values  (Figure  4.4,  Fit  1). 

The  trend  here  shows  an  increase  in  Iq  values  at  closer  in  ranges  for  Apple 
11-13  than  for  Turk-4.  We  did  note  that  as  reported  in  Reference  9 (WT-1123) 
the  damage  to  vehicles  (likewise  displacements)  on  Apple  11-13  was  not  as 
extensive  as  for  other  events.  The  reason  given  was  that  the  vehicles  were 
located  on  a surface  which  was  unusually  soft  and  sandy.  We  proceeded  to 
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compare  only  Apple  11-13  and  Met-12  wherein  the  surface  was  also  unusually 
soft  and  sandy,  particularly  along  the  desert  blast  line.  These  comparisons 
are  shown  in  Figures  4.5,  Fits  1 and  2 and  Figure  4.6,  Fit  3.  The  scatter  of 
data  is  about  the  same  for  Fits  1 and  3 and  somewhat  larger  for  Fit  2 (see 
Table  4.2).  When  deriving  the  Iq  values  for  Apple  11-13,  the  Iqavg  derived 
from  vehicle  displacements  and  averaged  with  IqQ  (measured)  tended  to  lower 
the  Iq  values  relative  to  IqQ.  The  opposite  of  this  tendency  was  the  case 
for  the  majority  of  Iq  on  other  events,  i.e.,  the  Iqavg  derived  from  vehicle 
displacements  averaged  with  Iqg  values  increased  Iq  relative  to  Iqg.  We, 
therefore,  compared  only  the  measured  Iq  on  Apple  11-13  with  Met-12,  Figure 
4.6.  Above  Iq  = 10  kPa-sec  the  comparison  is  poor  but  for  the  remaining  Iq 
values  there  is  good  agreement.  Note  that  the  percentage  difference  in  range 
for  the  three  Fits  is  small  for  iq  between  0.8  kPa-sec  and  10  kPa-sec.  If  we 
believe  the  measurements  (our  inclination  is  that  they  are  valid  and  reason- 
able) then  we  have  a dichotomy  with  regard  to  Apple  11-13. 

(U)  On  the  one  hand,  just  using  Iq  values  for  Apple  11-13  and  combining 

with  Turk-4  indicates  that  Apple  11-13  could  be  considered  a light  dust  event 
or  we  could  argue  that  Turk-4  is  a heavy  dust  event.  However,  Turk-4  Iq  data 
at  close  in  ranges  are  definitely  below  the  Met-12  data  (one  of  the  measures 
for  heavy  dust  and  light  dust)  so  Turk-4  was  classified  as  light  dust.  Apple 
11-13  when  combined  with  the  Met-12  data  could  be  classified  as  a heavy  dust 
event  (all  previous  evaluations  stress  the  fact  that  Met-12  was  a heavy  dust 
event)  since  there  is  some  agreement  when  the  data  from  these  two  events  are 
combined.  We  dwell  on  this  matter  to  point  out  that  the  data  per  se  for  a 
number  of  events  does  not  permit  a clear  distinction  between  light  dust  and 
heavy  dust  blast  waves  throughout  the  Iq  versus  range  curves.  All  one  can 
do  is  surmise  that  the  type  of  surface  conditions  will  produce  varying 
amounts  of  dust  content  in  the  blast  wave. 

(U)  Now  we  come  to  the  Iq  versus  range  curves.  Figures  4.7  to  4.13, 

which  we  classified  as  light  dust  events  since  all  the  Iq  data  on  these 
curves  are  lower,  particularly  at  the  close  in  range  than  the  heavy  dust 
events.  These  curves  are  for  HOB  ranging  from  46.8  metres  out  to  66.9 
metres.  If  one  were  to  superimpose  these  curves  on  top  of  one  another  they 
would  not  vary  by  more  than  13  percent  in  range  and  this  is  mainly  in  the 
region  of  Iq  from  0.6  kPa-sec  to  0.8  kPa-sec.  As  before,  in  this  set  of 
curves  we  made  use  of  the  Iq  measured  as  well  as  Iq.  It  is  to  be  noted  that 
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the  asphalt  results,  Met-12  have  also  been  included  with  Turk-4  and  Smoky-15. 
Here,  we  classify  Smoky-15  as  a light  dust  event.  The  asphalt  results  fit 
reasonably  well  with  the  other  data  (Figure  4.7).  We  observe  the  same  situa- 
tion when  we  Included  asphalt  results  from  Bee-6  with  Moth-2  and  Grable-10 
(Figure  4.13).  Again  the  asphalt  iq  values  from  Bee-6  include  the  Iqayg  from 
vehicle  displacements  over  the  asphalt  using  the  quadratic  fits  (Section 
3. 1.2. 4).  An  Interesting  observation  about  Grable-10  Is  that  this  event  was 
detonated  in  the  Frenchman’s  Flat  area  as  was  Met-12,  yet  we  consider  Grable- 
10  as  producing  light  dust-laden  blast  waves.  We  pointed  out  In  our  initial 
evaluation  (Reference  1)  that  on  Met-12  the  vehicle  response  along  the  desert 
line  (disturbed,  very  fine  loose  soil)  was  greater  than  the  vehicle  response 
along  the  Desert  Rock  sector  (same  soil  but  undisturbed).  For  Met-12,  the 
Iqavg  (calculated  from  vehicle  displacements)  was  greater  on  the  desert  line 
than  on  the  Desert  Rock  sector.  We  can  safely  assume,  therefore,  the  higher 
IqaVg  along  the  desert  line  was  the  result  of  the  disturbed  surface  (producing 
higher  dust  concentrations  in  the  blast  wave).  The  Frenchman's  Flat  soil  con- 
sists of  alluvial  deposits  which  are  fine  microsize  particles  similar  to  talcum 
powder.  The  desert  line  region  was  disturbed  by  construction  of  test  items, 
digging  up  cable  ditches,  and  high  concentration  of  daily  traffic  prior  to 
detonating  the  device.  These  are  somewhat  unique  conditions.  No  such  acti- 
vity took  place  in  the  Desert  Rock  sector.  The  vehicles  on  Grable-10  were 
located  in  the  Desert  Rock  sector  and  thus,  had  lower  values  of  Iq  equivalent 
to  the  Iq  values  of  Moth-2,  Apple  1-8  and  Smoky-15,  which  were  classified  as 
light  dust  events  (see  Figure  4.11).  We  also  note  that  Fits  1 and  2,  Figure 
4.11  nearly  coincide  over  the  total  range  of  Iq  values  and  this  is  the  same 
situation  for  Fits  1 and  2,  Figure  4.10. 

(U)  The  next  curve,  Figure  4.14,  is  a fit  to  the  data  obtained  on 

Grable-10,  Moth-2  and  Bee-6  (Fit  2).  Two  of  the  Iq  data  points  do  not  signi- 
ficantly alter  the  curve,  however,  as  for  example.  Fit  1 for  Moth-2  and  Bee-6, 
Figure  4.15.  The  curves  of  Fit  2 (Figure  4.14)  and  Fit  1 (Figure  4.15)  are 
about  the  same.  Since  these  two  points  did  not  necessarily  influence  the  fit 
we  decided  to  fit  the  data  points  by  eye.  The  reason  for  this  is  that  we 
assumed  that  Bee-6  could  be  considered  to  produce  heavy  dust-laden  blast 
waves.  The  desert  surface  on  Bee-6  was  a fine,  soft  sandy  soil  similar  to 
Apple  11-13.  However,  the  four  data  points  on  Bee-6  do  not  necessarily 
validate  this  assumption. 
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(U)  The  next  set  of  cur</e* , Figures  4.16,  4.17  and  4.18  deal  with  near- 

ideal blast  waves.  The  average  HOB  for  these  curves  range  from  127  metres  up 
to  240  metres.  The  HOB  * 127  metres  Is  represented  by  Yuma-4  (HOB  = 109 
metres),  and  Wasp  Prime-9  (HOB  ■ 145  metres)  shown  in  Figure  4.16,  Fit  1.  We 
do  not  feel  that  this  curve,  Fit  1,  is  representative  for  the  Iq  values  at 
this  average  HOB  for  several  reasons.  First,  we  mentioned  previously  that  we 
questioned  the  Iq  values  on  Yuma-4;  second,  the  scatter  of  the  data  is  large 
between  those  two  events;  and  third,  the  large  difference  In  HOB  (36  metres) 
between  these  two  events  could  be  a factor  in  the  noted  scatter.  We  feel, 
therefore,  that  in  spite  of  the  lack  of  data,  the  fit  by  eye  to  the  Wasp 
Prime-9  event  is  more  reasonable.  The  remaining  two  curves,  the  combined 
events,  Encore-9  and  Wasp-1,  Figure  4.17  and  Encore-9,  Figure  4.18,  do  not 
show  any  substantial  difference.  The  maximum  difference  between  these  curves 
is  less  than  9 percent  and  that  falls  close  to  the  scatter  for  each  curve 
(see  Table  4.2).  Because  of  the  uncertainty  of  the  data  on  Wasp-1  (vehicle 
orientations  and  small  displacements  which  poorly  determined  Iqavg  values)  we 
consider  the  Encore-9  curve  only  to  be  more  reliable. 

(U)  To  recap,  we  place  high  confidence  on  the  Iq  data  for  the  surface 

bursts  (ideal/near-ideal)  and  for  the  HOB  range  from  about  40  to  70  metres 
(light  dust).  There  is  less  confidence  in  the  height  of  burst  at  HOB  from 
109  metres  out  to  240  metres.  But  this  is  only  because  of  the  small  amount 
of  data  available  at  these  burst  heights.  Nevertheless,  the  Iq  values  seem 
reasonable,  except  for  Yuma-4.  The  heavy  dust  events,  represented  by  Bee-6, 
Met-12  and  Apple  11-13,  indicate  Iq  values  to  be  2 to  3 times  the  values  for 
the  light  dust  events  (compare  Figure  4.5,  Fit  1,  Heavy  Dust  with  Figures 
4.2  and  4.7).  The  surface  conditions  for  the  heavy  dust  events  producing 
these  blast  environments  were  unique,  particularly  the  Frenchman's  Flat  area. 
Only  under  these  conditions  would  a similar  blast  environment  be  produced  and 
these  surface  conditions  would  be  rare  to  find  in  many  other  parts  of  the 
world.  We  tend  to  believe  the  validity  of  the  measurements  (there  are  some 
exceptions  of  course)  but  the  Iq  which  combine  measurements  and  Iqavg  from 
vehicle  displacements  should,  in  general,  be  more  reliable  which  is  what  we 
wanted  to  show. 

4.2  (U)  DYNAMIC  PRESSURE  IMPULSE,  Iq,  HOB  CHARTS 

(U)  We  presented  various  scaled  Iq  versus  scaled  range  curves  for  three 
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types  of  blast  environments;  ideal/near-ideal,  light  dust,  and  heavy  dust 
(see  Section  4.1).  The  procedures  for  classifying  each  event  or  a group  of 
events  into  one  of  the  above  categories  were  outlined.  Use  was  made  of  these 
curves  for  each  type  of  blast  category  to  prepare  Tables  which  list  the  iq 
values  as  a function  of  range  for  different  burst  heights.  The  scaled  Iq 
values  in  these  Tables  extend  from  0.4  kPa-sec  out  to  10  kPa-sec.  The  scaled 
range  for  each  Iq  value  was  selected  by  eye  from  the  fitted  curves.  Tables 
4.7,  4.8,  and  4.9  give  this  information  for  ideal/near-ideal,  light  dust  and 
heavy  dust  blast  environments,  respectively.  The  Iq  values  given  are  for  a 
height  of  0.91  metres  (3  feet)  above  the  ground  surface. 

(U)  The  iq  versus  range  data  for  the  Ideal/near-ideal  blast  waves  are 

given  for  five  burst  heights  (Table  4.7).  Two  of  these  burst  heights  include 
the  Iq  values  from  Yuma-4  (HOB  = 127  metres)  and  Wasp-1  (HOB  * 224  metres) 
which  yielded  questionable  data  (see  Section  4.1).  However,  a comparison  of 
the  ranges  for  several  iq  values  for  HOB  * 127  metres  and  145  metres  shows 
the  percentage  difference  Is  less  than  9 percent.  The  exceptions  are  for  iq 
of  5 kPa-sec  and  4 kPa-sec.  Here,  the  differences  are  about  21  percent  for 
iq  = 5 kPa-sec  and  14  percent  for  iq  * 4 kPa-sec. 

(U)  The  next  set  of  iq  versus  range  data  is  for  the  light  dust  events 

(Table  4.8).  The  averaged  HOB  in  this  set  of  data  extends  from  41.5  metres 
out  to  66.9  metres.  If  we  neglect  the  extrapolated  values  of  iq  = 10  kPa-sec, 
the  maximum  percentage  difference  in  range  for  given  values  of  Iq  is  about  17 
percent  for  these  burst  heights  and  this  difference  is  due  primarily  to  the 
extrapolated  values  for  the  lower  iq  values.  If  we  further  neglected  all  of 
the  extrapolated  values  then  the  percentage  difference  would  be  about  10 
percent.  This  indicates  that  these  data  could  be  consolidated  into  one  aver- 
aged burst  height  and  least  squares  fitted  without  introducing  large  errors. 
(U)  The  heavy  dust  Iq  versus  range  data  spans  the  averaged  HOB  from 

43.4  out  to  73.2  metres  (Table  4.9).  In  this  set  of  data,  there  are  some 
events  that  we  classified  as  light  dust  which  have  been  combined  with  the 
heavy  dust  events.  In  Section  4.1,  we  indicated  the  difficulties  encountered 
in  attempting  to  distinguish  between  light  dust  and  heavy  dust  events  over 
the  total  Iq  versus  range  curves.  The  events  which  we  feel  clearly  fall  into 
the  heavy  dust  category  based  on  the  measures  used  for  distinguishing  light 
and  heavy  dust,  are  Met-12  and  Apple  11-13.  The  surfaces  for  both  these 
events  were  fine,  loose  sandy  soil  and  the  magnitudes  of  dynamic  pressure 
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impulse  were  2 to  3 times  greater  than  for  the  light  dust  events.  Bee-6  j 

could  fall  Into  this  category  as  well.  Similar  surface  conditions  existed 
but  the  maximum  difference  In  Iq  was  only  about  1-1/2  times  greater  than  for 
light  dust  and  this  was  over  a small  portion  of  the  range  (compare  Bee-6  and 
Moth-2,  Figure  4.15).  The  set  of  Iq  versus  range  data  from  Table  4.9  that 
could  then  be  Interpreted  to  be  representative  of  heavy  dust  events  are  only 
two,  HOB  * 44.6  metres  and  73.2  metres.  Between  these  two  burst  heights  the 
difference  in  range  for  constant  Iq  is  small,  less  than  8 percent.  Even  over 
this  whole  span  of  HOBs,  neglecting  the  extrapolated  values,  the  range  differ- 
ences for  constant  Iq  values  decreases  from  about  29  percent  down  to  10  j 

percent  with  decreasing  Iq. 

(U)  For  constructing  the  HOB  Charts,  the  Iq  versus  range  data  from  the 

three  types  of  blast  conditions  have  been  combined  and  placed  into  two  cate- 
gories — light  dust/near-ideal  and  heavy  dust/near-ideal  blast  waves.  The 
HOB  Charts  are  shown  in  Figures  4.19  and  4.20,  one  for  light  dust/near-ideal  ! 

and  the  other  for  heavy  dust/near-ideal  blast  environments,  respectively. 

The  burst  heights  selected  for  constructing  these  charts  were:  three  for  the  ' 

near-ideal,  0,  145,  and  240  metres;  two  for  the  light  dust,  49.2  and  61.1 
metres  (Fit  1);  and  two  for  the  heavy' dust,  44.6  (Fit  1)  and  73.2  metres  (eye- 
fit).  The  selection  of  only  two  HOBs  for  the  light  dust  events  was  mainly  to  ! 

minimize  clutter  of  data  on  the  charts.  As  we  pointed  out  for  the  light  dust  j 

events  there  were  small  differences  in  Iq  versus  range  curves  over  the  span  j 

of  HOB  from  41.5  to  66.9  metres.  As  mentioned  previously,  only  two  of  the  \ 

Iq  versus  range  curves  could  be  interpreted  to  be  heavy  dust  events.  Thus,  j 

the  five  burst  heights  selected  from  each  combination  were  sufficient  to,  at  ] 

least,  determine  the  trend  of  the  iso-Iq  contours.  The  error  bands  have  not  I 

been  placed  in  the  contours  but  these  can  be  readily  obtained  from  the  statis-  I 

tical  data  derived  by  least  squares  fit  procedures  (see  Table  4.5). 

(U)  Upon  inspection  of  these  Charts  we  note  the  effect  of  HOB  on  range 

for  equivalent  Iq  values.  As  the  burst  height  increases,  the  range  initially 
increases  and  then  decreases,  but  for  some  low  values  of  Iq  the  range 

increases  with  increasing  HOB.  Moreover,  the  iso-Iq  are  not  everywhere  para-  * 

11  el.  We  also  note  that  from  about  HOB  = 40  metres  to  70  metres  the  range  is 
reasonably  constant  for  given  Iq  values  for  either  the  light  dust  or  the 
heavy  dust  events. 
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(U)  Although  the  Iq  values  from  the  heavy  dust  events  (HOB  a 44.6 

metres)  were  2 to  3 times  higher  than  the  light  dust  events  (HOB  a 49.2 
metres)  the  comparison  of  Iq  as  a function  of  range  shows  that  the  range 
difference  Is  about  27  percent  greater  for  the  heavy  dust  events  at  the 
higher  Iq  than  for  the  light  dust  events.  The  range  difference  then 
decreases  to  about  10  percent  at  the  lower  Iq.  As  the  burst  height  Is 
increased  up  to  about  70  metres  the  range  difference  between  these  two  types 
of  blast  environments  becomes  smaller,  down  to  about  6 percent.  These  curves, 
heavy  dust  and  light  dust,  converge  at  the  lower  Iq  values. 

(U)  Here,  we  are  canparing  the  heavy  dust  data  from  three  events  which 

are  unique  and  one  of  these  (Bee-6)  has  only  limited  data  with  data  from  five 
light  dust  events  which  we  feel  are  reasonably  reliable.  Under  these  circum- 
stances and  since  the  range  differences  are  not  that  large,  our  tendency  is  to 
place  more  importance  on  the  HOB  Chart  which  Includes  the  light  dust  events 
than  the  chart  which  includes  the  heavy  dust  events.  Other  observations,  but 
somewhat  qualitative,  which  further  lead  us  toward  this  trend  are  the  surface 
conditions  of  the  light  dust  events,  the  pressure-time  measurements,  and  the 
agreement  between  Iq  derived  from  the  light  dust  events  and  the  asphalt  sur- 
faces of  Bee-o  -<id  Met-12.  The  light  dust  events  were  detonated  over  surfaces 
which  consir  ed  of  rough  gravel  mixed  In  with  sand  and.  In  some  cases,  vegeta- 
tion was  Interspersed  throughout  the  area.  For  all  of  these  events  except 
Smoky-15,  only  BRL  mechanical  self-recording  gages  were  used  along  the  blast 
line.  Cable  ditches  or  any  major  obstruction  were  not  required  to  install 
the  gages.  The  surface  conditions  were  more  or  less  undisturbed.  Additionally, 
a cursory  examination  of  the  pressure-time  measurements  made  on  the  light  dust 
events  shows  less  high  frequency  oscillations  than  for  the  heavy  dust  events. 
This  implies  that  the  influence  of  dust  on  these  measurements  (light  dust) 
was  minimal.  We,  therefore,  believe  that  the  main  influences  on  the  derived 
Iq  for  the  light  dust  events  was  the  result  of  precursor  action  and  not 
necessarily  the  dust.  This  is  somewhat  borne  out  by  the  Iq  values  obtained 
over  the  asphalt  surfaces  on  Bee-6  and  Met-12.  The  asphalt  Iq  were  combined 
with  Iq  from  the  light  dust  events  and  both  are  represented  reasonably  well 
by  a single  curve  for  that  specific  averaged  burst  height  (see  Figures  4.2 
and  4.13).  No  significant  amount  of  dust  was  present  over  the  asphalt  sur- 
face. For  these  reasons  we  tend  to  consider  the  HOB  Chart  which  includes  the 
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light  dust  to  be  more  representative  of  the  expected  blast  conditions  over 
various  surfaces  where  precursors  (non-ideal)  would  be  formed. 

(U)  Since  the  surface  conditions  (disturbed,  micro-size  particles)  of 

the  heavy  dust  events  such  as  on  Frenchman's  Flat  area  were  unique,  the  blast 
environments  produced  were  also  unique.  These  unique  surface  conditions  gave 
rise  to  heavy  dust  loading  of  the  blast  wave  (precursor)  producing  much  higher 
Iq  than  for  the  light  dust  events.  We,  therefore,  tend  to  place  less  impor- 
tance on  the  heavy  dust  events  and  perhaps  even  to  disregard  them. 
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TABLE  4.1  (U)  Dynamic  Pressure  Impulse  versus  Range, 

Final  Averages  — Ideal/Near-Ideal  Bias*  Waves,  Surface  Bursts 


UNCLASSIFIED 


QPERATION/EVENT 
Castle,  Koon-3 

Castle,  Nectar-6 


Redwing, 

Lacrosse-1 

Redwing,  Zunl-3 


ACTUAL 

SCALED 

RANGE 

RANGE 

Metres 

Metres 

1981 

372.2 

2499 

469.6 

3322 

624.2 

4206 

790.3 

2316 

193.6 

2408 

201.3 

2591 

216.6 

2743 

229.3 

2920 

244.1 

3990 

333.6 

762 

223.3 

844 

247.3 

1021 

299.2 

1189 

348.4 

1334 

390.9 

2135 

140.3 

2530 

166.2 

3170 

208.3 

3566 

234.3 

4206 

276.3 

5029 

330.4 

!qo  Iqavg 

kPa-sec  kPa-sec 

2.83  2.29 

1.26  1.62 

0.647 
0.346 

33.9  22.1 

29.4  21.6 

24.7  21.8 

20.7  19.8 

17.7 

7.73  5.38 

6.18  8.05 

4.78  5.58 

2.89  2.92 

1.82  2.22 

1.34  1.49 

77.4  62.5 

51.4  30.4 

28.3  37.4 

19.8  29.5 

12.3  14.6 

6.71  4.79 


ACTUAL  SCALED 

Iq  Iq 

kPa-sec  kPa-sec 


2.56 

0.494 

1.43 

0.278 

0.647 

0.125 

0.346 

0.067 

28.0 

2.41 

25.4 

2.18 

23.0 

1.96 

20.3 

1.74 

17.7 

• 1.52 

6.56 

0.563 

7.12 

2.14 

5.18 

1.56 

2.90 

0.873 

2.02 

0.608 

1.42 

0.427 

70.0 

4.72 

40.9 

2.76 

32.8 

2.21 

24.6 

1.66 

13.4 

0.903 

5.75 

0.389 
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TABLE  4.2  (U)  Dynamic  Pressure  Impulse  versus  Range, 
Final  Averages  — I deal /Near- I deal  Blast  Waves/Air  Bursts 


CONFIDENTIAL  - FORMERLY  RESTRICTED  DATA 


ACTUAL 

SCALED 

OPERATION/EVENT 

RANGE 

RANGE 

Metres 

Meters 

Upshot- Knothole, 

756 

245.2 

Encore-9 

872 

282.9 

933 

302.7 

*198 

388.6 

1329 

431.1 

1704 

552.8 

1996 

647.5 

Teapot,  Wasp-1 

354 

317.3 

393 

352.2 

434.5 

389.4 

* 

460.8 

413.0 

607.3 

544.3 

Teapot,  Wasp 

237 

152.3 

Prime-9 

311 

199.8 

465 

298.8 

618 

397.1 

Redwing,  Yuma-4* 

76 

132.8 

107 

187.0 

122 

213.3 

Iqo 

Iqavg 

ACIUAL 

Iq 

SCALED 

Tq 

kPa-ser. 

kPa-sec 

kPa-sec 

kPa-sec 

3.52 

3.17 

3.38 

1.24 

3.00 

2.84 

2.92 

1.07 

2.75 

3.01 

2.88 

1.06 

1.82 

2.45 

2.13 

0.781 

1.51 

1.64 

1.58 

0.579 

0.830 

0.870 

0.850 

0.312 

0.680 

0.759 

0.719 

0.264 

0.917 

0.979 

0.948 

0.978 

0,779 

..1.12 

0.951 

Q.981 

0.643 

0.80 

0.722 

0.745 

0.567 

1.05 

0.808 

0.834 

0.315 

0.688 

0.502 

0.518 

8.63 

9.24 

8.94 

6.83 

6.75 

8.38 

7.57 

5.78 

1.82 

-- 

1.82 

1.39 

0.871 

-- 

0.871 

0.665 

2.38 

3.23 

2.81 

5.07 

1.30 

1.85 

1.56 

2.82 

1.03 

1.85 

1.44 

2.60 

(U)  *Vield  from  the  above  information  makes  this  event  CONFIDENTIAL  - 
FORMERLY  RESTRICTED  DATA 
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TABLE  4.3  (U)  Dynamic  Pressure  Impulse  versus  Range, 
Final  Averages  — Non-Ideal  Blast  Waves/Air  Bursts 

UNCLASSIFIED 


ACTUAL 

RANGE 

SCALED 

RANGE 

■<10 

^avg 

ACIUAL 

iq 

SCALED 

Iq 

OPERATION/EVENT 

Metres 

Metres 

kPa-sec 

kPa-sec 

kPa-sec 

kPa-sec 

Upshot-Knothole, 

736 

287.7 

5.13 

5.66 

5.39 

2.37 

Grable-10 

844 

329.9 

2.69 

2.53 

2.61 

1.11 

(Desert  Rock) 

1335 

433 

— 

0.686 

0.686 

0.302 

Teapot,  Moth-2 

287 

204.1 

6.75 

6.75 

5.39 

366 

260.3 

3.35 

3.86 

3.61 

2.88 

411 

292.3 

2.42 

2.05 

2.24 

1.79 

457 

325.0 

1.88 

1.47 

1.68 

1.34 

518 

368.4 

1.35 

— 

1.35 

1.08 

Teapot,  Turk-4 

716 

193.0 

20.7 

20.7 

6.51 

914 

246.4 

11.8 

10.5 

11.2 

3.52 

1030 

277.7 

7.69 

9.81 

8.75 

2.75 

1128 

304.1 

5.66 

7.77 

6.72 

2.11 

1280 

345.1 

2.32 

— 

2.32 

0.730 

Teapot,  Bee-6 

549 

263.7 

8.10 

9.36 

8.73 

4.79 

(Desert) 

610 

293.0 

4.17 

4.44 

4.31 

2.37 

701 

336.7 

2.32 

2.59 

2.46 

1.35 

777 

373.2 

1.9! 

2.66 

2.28 

. 1.25 

Teapot,  Apple  1-8 

527 

205.6 

17.3 

. - — 

17.3 

7.95 

622 

242.6 

9.10 

-- 

9.10 

4.18 

902 

351.9 

1.95 

2.62 

2.29 

1.05 

991 

386.6 

1.20 

1.77 

1.48 

0.68 

1128 

440 

0.38 

0.76 

0.57 

0.262 

Teapot,  Met-12 

610 

208.9 

43.0 

41.0 

42.0 

16.4 

(Desert) 

686 

235.0 

37.6 

41.1 

39.3 

15.3 

762 

261.0 

16.8 

15.7 

16.3 

6.36 

838 

287.0 

13.8 

15.2 

14.5 

5.66 

914 

313.0 

3.51 

3.57 

3.54 

1.38 

Teapot,  Apple  11-13 

518 

160.3 

41.9 

_ _ 

41.9 

15.4 

625 

193.4 

34.0 

— 

34.0 

12.5 

808 

250.0 

19.3 

-- 

19.3 

7.08 

914 

282.8 

10.3 

7.81 

9.06 

3.32 

1006 

311.3 

5.54 

3.66 

4.60 

1.19 

1128 

349.0 

3.98 

3.23 

3.61 

1.32 

1219 

377.2 

2.55 

2.81 

2.68 

0.983 

Plumbbob,  Smoky- 15 

457 

122.6 

138 

138 

44.7 

726 

194.8 

28.5 

36.8 

32.6 

10.6 

789 

211.7 

23.7 

— 

23.7 

7.68 

841 

225.6 

16.0 

14.1 

15.1 

4.90 

897 

240.7 

12.5 

7.70 

10.1 

3.27 

1038 

278.5 

6.90 

7.22 

7.06 

2.29 

1181 

316.9 

3.75 

3.59 

3.67 

1.19 
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TABLE  4.4  (U)  Dynamic  Pressure  Impulse  versus  Range, 

Final  Averages,  Asphalt  Surface  — Non-Ideal  Blast  Waves/Air  Burst 


UNCLASSIFIED 


ACTUAL 

SCALED 

RANGE 

RANGE 

OPERATION/EVENT 

Metres 

Metres 

Teapot,  Bee-6 

549 

263.7 

610 

293.0 

701 

336.7 

777 

373.2 

Teapot,  Met- 12 

610 

208.9 

686 

235.0 

762 

261.0 

838 

287.0 

914 

313.0 

^0 

Iqavg 

ACIUAL 

Iq 

SCALED 

kPa-sec 

kPa-sec 

kPa-sec 

kPa-sec 

6.59 

6.97 

6.78 

3.72 

4.16 

2.51 

3.38 

1.85 

1.60 

2.10 

1.35 

1.02 

1.02 

1.24 

1.13 

0.62 

19.0 

18.9  . 

19.0 

7.42 

13.7 

14.8 

14.2 

5.54 

8.40 

8.34 

8.37 

3.27 

5.58 

4.38 

4.98 

1.94 

1.72 

1.37 

1.54 

0.60 
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TABLE  4.5  (U)  LEAST  SQUARES  FIT  RESULTS 
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TABLE  4.6  (U)  Dynamic  Pressure  Impulse  versus  Range 
for  Met-12  and  Apple  11-13  — Measured  Values 

UNCLASSIFIED 


ACTUAL 

Iqo 

SCALED 

RANGE 

RANGE 

Metres 

kPa-sec 

MET-12 

Metres 

610 

48.0 

208.9 

686 

37.2 

235.0 

762 

17.2 

261.0 

838 

11.5 

287.0 

914 

3.4 

313.0 

*%> 

kPa-sec 


18.7 

14.5 

6.71 

4.49 

1.33 


APPLE  11-13 


518 

41.9 

625 

34.0 

808 

19.3 

914 

12.1 

1006 

7.4 

1128 

4.8 

1219 

2.48 

160.3 

15.4 

193.4 

12.5 

250.0 

7.08 

282.8 

4.44 

311.3 

2.71 

349.0 

1.76 

377.2 

0.910 
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TABLE  4.7  (U)  Helght-Of-Burst  For  Dynamic  Pressure  Impulse  Versus 
Ground  Range  At  0.91  Metres  Above  Surface 

Ideal /Near-Ideal 
UNCLASSIFIED 


Iq 

(kPa-sec) 

lOudIMI 

240 

10 

87* 

-- 

— 

-- 

— 

5 

128* 

158 

191 

— 

-- 

4 

143 

184 

210 

— 

— 

3 

174 

218 

237 

— 

— 

2 

212 

260 

277 

— 

1.5 

240 

300 

307 

— 

— 

1.2 

261 

325 

330 

258 

262 

1.0 

281 

346 

350 

317 

301 

0.8 

306 

372 

372 

382 

353 

0.6 

341 

402 

402* 

458 

422 

0.4 

395 

445* 

445* 

553 

522 

* Extrapolated 
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. TABLE  4.9  (U)  Height-Of-Burst  For  Dynamic  Pressure  Impulse 

Versus  Ground  Range  At  0.91  Metres  'hove  Surface 

Heavy  Dust/ Non -I deal 
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FIGURE  4.3  (U)  Dynamic  Pressure  Impulse  vs.  Range  - Average  HOB  = 

43.4  Metres,  Scaled  to  1 KT 
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FIGURE  4.4  (U)  Dynamic  Pressure  Impulse  vs.  Range  - Average  HOB  = 

44.1  Metres,  Scaled  to  1 KT 
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FIGURE  4.7  (U)  Dynamic  Pressure  Impulse  vs.  Range  - Average  HOB  = 46 

Metres,  Scaled  to  1 KT 
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FIGURE  4.9  (U)  Dynamic  Pressure  Impulse  vs.  Range  - Average  HOB  = 52.6 
Metres,  Scaled  to  1 KT 
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FIGURE  4.10  (U)  Dynamic  Pressure  Impulse  vr.  Range  - Average  HOB 
Metres,  Scaled  to  1 KT 
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FIGURE  4.11  (U)  Dynamic  Pressure  Impulse  vs.  Range  - Average  HOB  * 61.1 

Metres,  Scaled  to  1 KT 
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FIGURE  4.19  (U)  Dynamic  Pressure  Impulse  vs.  Ground  Range  at  0 91  Metres  Above 
Surface  — Light  Dust/liear-Ideal  Surface,  Scaled  to  1 KT 


FIGURF  4.20  (U)  Dynamic  Pressure  Impulse  vs.  Ground  Range  at  0.91  Metres  Above 
Surface  — Heavy  Dust/Near-Ideal  Surface,  Scaled  co  1 KT 
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SECTION  5 

(This  Section  is  Unclassified) 

SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

5.1  SUMMARY  AND  CONCLUSIONS 

5.1.1  Dynamic  Pressure  Impulse  versus  Displacement 
Displacement  versus  Dynamic  Pressure  Impulse 

Two  types  of  functions  have  been  used  for  least  squares  fit  to 
dynamic  pressure  impulse  versus  displacement  and  vice  versa.  The  power  func- 
tions were  used  to  separately  fit  the  data  for  the  WWII  and  the  M38A1  1/4  ton 
trucks.  After  further  analysis  in  which  we  pointed  out  that  these  two  vehicles 
could  be  treated  as  one  type  of  truck,  we  combined  the  data  for  the  1/4  ton 
trucks  as  well  as  for  the  2-1/2  ton  trucks  (M35  and  M135)  and  used  quadratic 
functions  for  least  squares  fits.  Also,  the  Iq  values  for  the  non-ideal 
events  were  iterated  and  used  in  the  latter  least  squares  fits.  The  resultant 
equations  for  Iq  versus  displacement  are: 

Combined  Vehicles:  M38A1  and  WWII  - 

SO  in  Iq  = 0.133  + 0.376  £n  D + 0.0566  (In  D)2 

FO  In  Iq  = 0.732  + 0.465  In  D + 0.0220  (In  D)2 

Combined  Vehicles:  REO  M35  and  GMC  Ml 35  - 

SO  In  Iq  * 0.429  + 0.333  In  D + 0.0275  (in  D)2 

and  for  displacement  versus  Iq  the  equations  are: 

Combined  Vehicles:  M38A1  and  WWII  - 

SO  In  D » -0.223  + 2.222  in  Iq  - 0.214  (in  Iq)2 

FO  in  D = -1.362  + 2.081  in  Iq  - 0.0998  (in  Iq)2 

Combined  Vehicles:  REO  M35  and  GMC  Ml 35  - 
SO  in  D = -1.365  + 3.339  in  Iq  - 0.4259  (in  Iq)2 

We  find  good  correspondence  between  dynamic  pressure  impulse  and  dis- 
placement, in  spite  of  the  large  scatter  of  data.  A number  of  uncontrolled 
variables  present  in  the  experiments,  which  we  described,  contributed  to  the 
scatter.  These  variables  are  statistically  averaged  out  in  the  use  of  least 
squares  fit.  One,  therefore,  can  calculate  with  reasonable  confidence  the  Iq 
from  known  vehicle  displacements  and  calculate  displacements  (such  as  onset 
of  motion)  for  given  values  of  dynamic  pressure  impulse. 
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5.1.2  Displacement  versus  Scaled  Ground  Range 

For  different  combinations  of  events  we  fitted  the  displacement 
versus  scaled  ground  range  using  various  values  of  n (the  exponent  of  W used 
in  scaling  ground  range).  The  best  value  of  n is  somewhat  dependent  on  the 
burst  height  and  surface  conditions.  The  values  of  n for  best  fit  of  data 
were  from  n = 0.4  to  n ■ 0.43.  Our  reason  for  these  fits  was  to  determine 
the  feasibility  of  scaling  the  ground  range  for  displacement  since  we  correla- 
ted vehicle  damage  with  displacement  (Reference  17).  Thus,  for  known  displace- 
ments as  a function  of  scaled  ground  range  damage  can  be  then  predicted  for 
different  yields.  The  utility  of  this  scheme  needs  further  evaluation. 

5.1.3  Dynamic  Pressure  Impulse  (Unsealed  and  Scaled)  versus  Scaled 

Ground  Range 

We  again  fitted  the  unsealed  dynamic  pressure  impulse  versus  scaled 
ground  range  using  different  values  of  n.  The  values  of  n which  best  fitted 
the  data  were  n = 0.42  or  n = 0.43.  These  were  then  compared  to  the  scaled 
dynamic  pressure  impulse  versus  scaled  ground  range.  We  concluded  from  an 
overall  viewpoint  there  is  little  difference  in  reliability  between  the  two 
sets  of  curve  fits  to  the  data.  TJie  scaled  or  unsealed  impulse  versus  scaled 
ground  range  curve  fits  may  be  used  with  equal  confidence. 

5.1.4  Surface  Effects  on  Dynamic  Pressure  Impulse,  Displacement  and 

Damage 

Comparing  the  data  from  water,  asphalt,  and  Desert  Rock  for  the 
same  events,  (Bee-6  and  Met-12)  we  conclude  the  following: 

1)  The  Iq  for  asphalt  are  lower  than  those  for  water  and  both  are 
lower  than  for  the  desert  line,  (Met-12).  On  Bee-6,  again,  the  Iq  for  asphalt 
are  lower  than  for  the  desert  line. 

2)  Displacements  of  vehicles  on  the  water  line  are  similar  to 
the  desert  line,  (Met-12)  while  the  displacements  on  asphalt  are  smaller  than 
on  desert  surface.  The  vehicle  displacements  on  the  Desert  Rock  are  about  the 
same  as  on  the  asphalt  surface  (Bee-6  and  Met-12). 

3)  Comparison  of  displacement  versus  dynamic  pressure  impulse 
for  desert  and  non-desert  cases  shows  that  there  is  no  substantial  difference 
among  desert,  water  and  asphalt  surfaces. 

4)  Damage  on  the  water  surface  was  somewhat  lower  than  on  the 
desert  surface  while  the  damage  on  the  asphalt  surface  on  the  average 

was  somewhat  higher  than  on  the  desert  surface.  The  damage  on  Desert  Rock 
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for  a given  displacement  was  higher  than  that  on  the  desert  line.  Thus, 
for  a given  Iq  value,  displacement  is  at  least  approximately  independent  of 
the  surface  on  which  the  vehicle  is  placed  but  for  a given  displacement  damage 
is  greater  on  the  harder  surfaces. 

5.1.5  Scaled  Dynamic  Pressure  versus  Scaled  Ground  Range 

We  used  the  once-iterated  Iq  values  (ideal/near-ideal,  blast 
waves)  and  twice-iterated  Iq  values  (non-ideal  blast  waves)  for  constructing 
scaled  Iq  versus  range  for  several  combinations  of  yield  and  burst  heights. 

For  most  of  the  events  we  were  able  to  distinguish  the  ideal/near-ideal  blast 
environments  from  non-ideal  blast  environments.  The  non-ideal  blast  environ- 
ments were  further  divided  into  light  dust  and  heavy  dust  events. 

We  have  high  confidence  in  the  scaled  Iq  versus  range  data  for  the 
surface  bursts  (ideal/near-ideal ) and  for  the  HOB  span  from  about  40  to  70 
metres  (light  dust  events).  There  is  less  confidence  on  the  Iq  data  above  70 
metres  but  this  is  only  because  of  the  small  amount  of  data  available  above 
this  burst  height. 

We  concluded  that  the  surface  conditions  for  the  heavy  dust  events 
were  unique  and  similar  surface  conditions  would  be  rare  in  other  parts  of  the 
world..  We  also  concluded  that  dynamic  pressure  measurements  were  reasonable. 
More  importantly,  combining  the  Iq  results  from  gages  and  those  calculated 
from  vehicle  displacement  results  in  greater  reliability  and  accuracy  in  the 
Iq  values.  We  also  conclude  that  the  Iq  data  obtained  on  the  asphalt  surface 
(Bee-6  and  Met-12)  fit  reasonably  well  with  data  from  the  light  dust  events 
having  similar  burst  heights. 

5.1.6  Dynamic  Pressure  Impulse,  Iq,  HOB  Charts 

Two  HOB  charts  have  been  constructed,  one  for  light  dust/near- 
ideal  and  one  for  heavy  dust/near-ideal.  The  light  dust  and  heavy  dust  portions 
of  the  charts  apply  to  HOB  from  40  to  70  metres  and  the  near-ideal  apply  to 
surface  bursts  and  to  HOB  above  70  metres. 

A comparison  of  heavy  dust  and  light  dust  events  shows  that  the 

Iq  values  for  heavy  dust  are  from  2 to  3 times  the  Iq  values  for  light  dust 

events.  In  terms  of  range,  however,  for  equal  Iq  values  the  heavy  dust  is 
greater  by  about  27  percent  at  the  higher  Iq  than  for  light  dust  and  then 
they  both  converge  to  about  the  same  range  for  lower  Iq  values. 

In  the  HOB  region  from  40  to  70  metres  the  range  for  iso-Iq 

seem  to  be  about  the  same  for  either  the  light  dust  or  heavy  dust  events. 
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Because  of  unique  surface  conditions  for  the  heavy  dust  events 
the  blast  environments  produced  were  unique.  Our  tendency  is  to  place  less 
importance  on  the  heavy  dust  events  and  perhaps  even  to  disregard  them.  We 
consider  the  HOB  charts  which  include  the  light  dust  events  to  be  more  repre- 
sentative of  the  expected  blast  conditions  over  various  surfaces  where  pre- 
cursors (non-ideal)  would  be  formed.  We  feel  that  the  resultant  Iq  values  for 
these  events  (light  dust)  were  caused  principally  by  the  precursor  action 
(based  on  comparison  of  asphalt  and  light  dust  Iq  values).  The  conclusion 
is  that  the  light  dust/near-ideal  HOB  chart  (nomenclature  here  may  now  be  a 
misnomer)  is  the  preferred  chart  for  application  to  other  targets. 

5.2  RECOMMENDATIONS 

1)  A procedure  similar  to  that  presented  should  be  used  for  other  drag- 
type  targets  to  extend  the  region  for  Iq  to  higher  values  than  10  kPa-sec 
and  for  Iq  to  lower  values  than  0.4  kPa-sec.  To  obtain  the  higher  values, 
heavier  and  harder  drag-type  targets  such  as  tanks  & howitzers  can  be  used. 

For  the  lower  values,  lighter  and  softer  targets  such  as  antennas  & coirmuni ca- 
tion equipment  can  be  used.  Additionally,  different  drag-type  targets  were 
exposed  to  nuclear  blast  at  burst  heights  other  than  those  used  here. 

The  evaluation  of  this  information  would  allow  greater  confidence  in  the 
construction  of  the  iso-Iq  contours. 

2)  There  is  a need  to  develop'  predictions  of  dynamic  pressure  wave- 
forms for  the  non-ideal  blast  environment.  Current  waveform  predictions  are 
based  on  the  data  obtained  on  heavy  dust  events  which  we  consider  to  be 
unique.  Application  of  waveform  predictions  is  needed  not  only  for  deter- 
mining loading  and  response  of  drag-type  targets  but  diffraction  type  targets 
as  well. 

3)  After  a thorcgh  evaluation  of  data  for  increasing  the  confidence  in 
constructing  HOB  charts  for  Iq  there  is  a need  to  construct  composite  HOB 
charts  representing  overpressure,  dynamic  pressure,  static  pressure  impulse 
and  dynamic  pressure  impulse.  Various  types  of  targets  are  sensitive  to 
different  blast  mechanisms  causing  failure.  For  targeting  and  vulnerability 
evaluations  an  HOB  chart  of  this  nature  would  show  the  advantages  or  disad- 
vantages against  an  array  of  targets,  particularly  when  constructed  by  colla- 
teral damage  (non-damage  or  limited  damage  to  friendly  targets). 
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APPENDIX  A 

(This  Section  is  Unclassified) 

I DEAL/NEAR- I DEAL  DYNAMIC  PRESSURE  IMPULSE 

The  procedure  selected  to  provide  us  with  initial  dynamic  pressure  impulse 
values  for  ideal/near-ideal  blast  environments  to  be  correlated  with  vehicle 
displacements  was  based  on  theoretical  calculations  and  analytical  techniques 
along  with  the  use  of  measured  blast  parameters  (Reference  1).  From  the  theo- 
retical calculations  and  the  analytical  techniques  we  obtained  the  ratio  of 
static  pressure  impulse  to  dynamic  pressure  impulse,  Ip/Iq,  as  a function  of 

overpressure,  P...  The  theoretical  calculations  consisted  of  the  HULL  code 

1 /3 

computations  for  Dice  Throw  ard  for  1 KT  nuclear  detonation  at  HOB  = 60  W ' . 
The  analytical  technique  used  was  that  of  Ethridge  (Reference  4)  for  a 500 
ton  HE  hemispherical  charge.  These  ratios  of  Ip/Iq  as  a function  of  over- 
pressure are  shown  in  Figure  A.l. 

We  pointed  out  in  Reference  1 that  the  best  measurements  of  blast  para- 
meters which  were  reasonably  reliable  were  overpressure  and  static  pressure 
impulse.  Thus,  for  each  ideal/near-ideal  blast  event,  we  obtained  the  average 
values  of  the  measured  overpressure  and  static  pressure  impulse  at  each  vehi- 
cle location.  These  values  were  then  used  in  conjunction  with  Figure  A.l  to 
obtain  the  dynamic  pressure  impulse. 

The  curve  designated  as  1 KT  60  W^3  (Figure  A.l)  was  used  to  derive  the 
ratios  of  Ip/Iq  as  a function  of  Ps  for  Upshot/Knothole,  Encore-9  and  Redwing, 
Yuma-4.  The  scaled  height  of  burst  for  these  two  events  were:  Encore-9, 

240  W 1//3  metres  and  Yuma-4,  109  W^3  metres.  The  curve,  NHE,  was  used  for 
Castle,  Koon-3  and  Nectar-6;  Redwing,  Lacrosse-1  and  Zuni-3;  and  for  the 
Canadian  20  Ton  and  100  Ton  events.  Both  of  the  Canadian  events  were  hemi- 
spherical and  burst  at  ground  surface.  The  nuclear  events  were  surface  bursts 
and  it  is  to  be  noted  that  the  NHE  curve  is  a good  approximation  of  a 1 KT 
nuclear  weapon  burst  on  the  surface.  And,  of  course,  for  Dice  Throw  the  HULL 
calculations  were  used. 

The  values  of  dynamic  pressure  impulse  derived  by  this  procedure  were 
checked  using  the  Rankine-Hugoniot  and  empirical  relations.  The  percentage 
difference  between  these  two  procedures  in  the  majority  of  cases  was  about  10 
percent,  although  some  were  as  large  as  25  percent. 

The  initial  estimates  of  dynamic  pressure  impulse  for  the  nine  events 
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are  shown  in  Figures  A. 2 through  A. 10  (solid  lines).  We  also  plotted  the 
first  iteration  of  Iq  on  these  Figures  (see  Section  4.1).  We  note  again  that 
IqaVg  which  was  combined  with  Iq^  to  obtain  Iq,  was  derived  from  vehicle  dis- 
placements using  the  quadratic  fits  (Section  3. 1.2. 4).  Overall  the  agreement 
between  the  initial  estimates  and  Iq  are  good  with  some  exceptions.  One  of 
the  exceptions  is  Yuma-4  which  we  discussed  previously  (Section  4.1).  Other 
exceptions  are  1 or  2 points  such  as  on  Nectar-6  and  Zuni-3  where  the  Iq  is 
+18  to  22  percent  from  the  initial  estimates. 
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APPENDIX  B 

(This  Section  is  Unclassified) 

DETERMINATION  OF  DYNAMIC  PRESSURE  IMPULSE  FOR  NON-IDEAL  SHOTS 

B.1  DATA  REDUCTION  PROCESS  # 

The  shots  Involved  In  this  section  were  all  part  of  Operation  Teapot. 

We  used  two  basic  procedures  and  some  variants  of  these  to  determine  dynamic 
pressure  Impulse  from  the  measured  pressure-time  histories.  The  measurements 
consisted  of  total  head  overpressure  and  static  overpressure  versus  time,  the 
latter  at  two  locations:  one  at  ground  level  and  one  at  the  same  height, 
usually  0.91  metres  (3  feet),  as  used  to  measure  total  head  pressure.  In 
addition,  previously  conducted  reduction  of  this  data  (using  only  the  two 
gages  at  the  same  height,  not  the  ground  level  measurement)  provided  dynamic 
pressure  versus  time  waveforms  (not  corrected  to  account  for  the  perturbation  to 
the  pressure  undergoing  measurement  by  the  gage  used  to  make  the  measurement). 

In  the  following  we  describe  the  correction  process  for  the  dynamic  pressure 
provided  and  an  independent  data  reduction  method  we  used  to  calculate  dynamic 
pressure  impulse  from  the  original  gage  data. 

We  use  the.  foil  owing  nomenclature  (Reference  3): 

Ps  - absolute  free  stream  static  pressure 
Pj  - absolute  free  stream  total  pressure 
Pp  - absolute  total  head  Pitot  pressure 
P0  - ambient  pressure 

APS  = ps  “ po  'free  stream  static  overpressure 

APy  * Pj  - P0  — free  stream  total  overpressure 

APp  = Pp  - Po  — total  head  overpressure 

q - dynamic  pressure 

Q-j  - dynamic  pressure  impulse;  subscript  i denotes  ith  method  of 
obtaining  result  (4  values  of  i used) 

Y - ratio  of  specific  heat  at  constant  pressure  to  specific  heat 
at  constant  volume  for  air 

M - local  free  stream  Mach  No.  of  flow  behind  blast  front 
At^  - ith  time  interval  (in  a pressure-time  waveform) 

G(M)  - Mach  correction  factor 

la(Wy  - average  value  of  G(M)  over  pressure-time  history 
* There  was  one  exception  to  this:  Plumbbob,  Smoky-15. 
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Wf  - weighting  factor  for  1th  time  interval 


To  determine  dynamic  pressure  Impulse  when  q(t)  Is  provided,  without 
Mach  correction,  we  must  calculate  the  Mach  correction  and  integrate  the 
corrected  dynamic  pressure  over  time.  This  involves  some  of  the  operations 
used  in  calculating  the  dynamic  pressure  impulse  from  the  raw  data  (pressure 
waveforms),  so  we  start  with  the  latter. 

It  can  be  shown  (Reference  3)  that 


(1) 


and  that  for  y * 1.4  and  M < 2.24  this  can  be  replaced  by 


PT(t)  - Ps(t) 

(2) 

q(t)  = 

1 J 

GTMT 

where  G(M)  = 1 + 1 M2  + ^ M4  + 

1 M6  + 

1600  M * ■ * 

(3) 

_X_ 

Also  PT/P$  = (1  + M2)7"' 

(4) 

For  M < 1 , Pp  = PT  . 

(5) 

However,  for  M > 1 , Pp  f Pj  so  that  to  determine  PT  (from  the  gage  result 

P„)  we  must  first  find  the  Mach  No.,  M. 

P 

To  do  this,  following  Reference  3,  we  have 


Thus,  for  M 1,  we  can  determine  M from  Equation  (4)  since  Pj  and  P$  are 
known  (as  functions  of  time).  Then  G(M)  is  determined  from  Equation  (3)  and 
q(t)  from  Equation  (2).  Breaking  up  the  waveforms  into  time  intervals  (the 
same  set  of  intervals  for  Pj  and  both  Ps  gages)  and  using  this  method  then 
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yields  a set  of  values  of  q (tj).  Then  Q * /q  dt  s Eq(t^)At1 


(7) 


For  M > 1,  we  can  determine  M from  Equation  (6),  then  from  Equat’  (4) 
since  Pp  and  P^  are  known  (as  functions  of  time).  The  same  procedure  as  above 
then  determines  Q.  (y  Is  known  as  a function  of  the  static  overpressure,  see 
Reference  11.) 

In  practice,  the  overpressure- time  rather  than  the  absolute  pressure- 
time histories  are  the  data  provided.  Therefore,  we  wished  to  transform  the 
above  relations  to  a form  suitable  for  using  the  measured  pressures  directly. 
To  do  this  we  rewrite  Equations  (4)  and  (6): 


PT 

pj  - (M,  y),  M < 1 


(4a) 


- f2  (M,  y),  M > 1 


and  note  that  Equation  (2)  can  be  written: 


q(t)  = 


APT(t)  - AP$(t) 


~gTmT 


Equation  (4a)  may  now  be  written: 


(6a) 


(2a) 


4pj  + P0 

' > 


(M,  y) 


or,  APt  * f1  (M,  y)  APs  + P0  (M,  y)  - 1],  M < 1 (8) 

and  similarly,  APp  = f2  (M,y  ) AP$  + P0  [f2  (M,  y)  - 1],  M > 1 . 

Therefore,  for  fixed  values  of  M,  y and  P0  the  desired  relations.  Equations 

(3) ,  are  straight  lines.  We  plotted  these  relations  using  APp  as  ordinate, 
AP<j  as  abscissa,  calculating  f^  and  f2  from  the  right  hand  sides  of  Equations 

(4)  and  (6)  for  assumed  values  of  M and  y.  We  made  separate  linear  plots  for 
each  decade  of  static  overpressure,  APS>  This  is  convenient  since  y can  be 
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obtained  as  a function  of  static  overpressure  (Reference  11).  For  1 <_  Ps  < 

100  psl  we  took  y * 1.4.  For  100  < £ 1000  psl  we  took  y 3 1.35. 

It  turns  out  that  there  were  few  Instances  of  concern  to  us  for  which 

Ps  > 100  psl.  Further,  It  turns  out,  as  will  become  apparent  when  we  con- 
sider Inaccuracies  Involved  In  the  data  reduction  process,  that  changes  in 
both  y and  PQ  induce  changes  In  the  calculated  results  which  are  negligibly 
small  compared  with  unavoidable  inaccuracies  In  the  data  reduction  process. 

We  assumed  P0  * 14  psl.  Also  for  M 3 1.3,  the  difference  between  Pp  and  PT 

Is  completely  negligible  compared  with  Inaccuracies  involved  in  data  reduction. 

We  plotted  Equations  (8)  In  the  manner  indicated  for  values  of  M between 
0.2  and  2.0  in  increments  of  0.2.  However,  in  reducing  the  Operation  Teapot 
data  we  only  rarely  encountered  a value  of  M > 1.2,  so  that  the  Mach  correc- 
tions owing  to  perturbation  of  the  flow  by  the  measuring  instrument  are,  in 
general,  quite  modest*  See  Equation  (3). 

Prior  to  initiating  the  data  reduction  process  we  carefully  examined 
the  pressure-time  waveforms  (Reference  12)  in  order  to  devise  a systematic 
method  with  as  good  accuracy  as  can  be  achieved.  We  noted  several  sources 
of  Inaccuracy:  the  waveforms  are  small , .especially  on. the  pressure  axis, 

typically  0.5"  peak  value;  the  areas  under  the  pressure-time  curves  are 
small,  usually  less  than  1 square  inch,  this  being  especially  true  for  the 
(already  reduced,  but  not  Mach  corrected)  dynamic  pressure-time  curves  where 
the  areas  under  the  curves  are  typically  in  the  range  of  0.2  to  0.4  square 
inches  (with  a few  smaller  and  a few  larger  than  this);  in  many  cases  the 
waveforms  are  highly  irregular,  some  of  the  irregularities  probably  being 
real,  i.e.,  physically  present  at  the  gage,  and  some  probably  being  instru- 
mental; there  are  some  problems  which  appeared  to  be  due  to  baseline  drift 
during  recording  and  some  which  may  well  have  been  due  to  gage  clogging  (with 
dust). 

The  previous  data  reduction  process,  which  provided  the  available  dynamic 
pressure-time  histories,  used  a point-by-point  method;  i.e.,  values  of  total 
and  static  pressures  were  taken  at  identical  (or  almost  identical)  times  from 
the  raw  data  and  the  dynamic  pressure  calculated  at  each  such  time  (as  the 
difference  between  the  total  and  static  overpressures).  Sometimes  the  dynamic 

* In  reducing  the  data  for  the  Plumbbob,  Smoky-15  event,  Equations  (4)  and 
(6)  were  used.  Plots  of  these  Equations  are  given  in  Figures  B.l  and  B.2. 
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pressure,  as  obtained  In  this  manner,  was  negative  — owing  presumably  to 
the  difficulties  enumerated  in  the  previous  paragraph. 

Any  attempt  to  use  the  pressure-time  waveforms  In  the  form  in  which  they 
are  now  available  to  obtain  a dynamic  pressure-time  curve  via  a point-by- 
point  data  reduction  method  seemed  to  us  to  be  doomed  to  failure.  We  con- 
cluded that  better  accurac^  (though  not  nearly  as  good  as  could  be  desired) 
could  be  achieved  by  use  of  a planlmeter  to  measure  areas  under  the  curves 
In  selected  time  Intervals. 

We  used  the  sets  of  parametric  curves  (straight  lines)  corresponding  to 
Equations  (8)  in  conjunction  with  visual  inspection  of  each  set  of  waveforms 
(3  waveforms:  Pp  and  two  Ps,  only  one  of  the  latter  having  been  used  In  the 

previous  data  reduction  process)  to  determine  approximate  average  values  of 
Mach  Ho.  over  trial  Intervals.  Where  average  values  of  M were  relatively 
large  or  rapidly  varying  the  time  intervals  were  taken  fairly  small.  We 
also  looked  visually  at  the  waveforms  In  doing  this  in  order  to  properly 
average  out  effects  of  irregularities  in  planimeterlng  the  areas  under  the 
curve  for  each  time  interval.  When  no  (or  only  small)  irregularities  were 
present  we  took  the  relative  contribution  of  each  time  interval  tothe  total 
area  under  the  pressure-time  curve  into  account  in  choosing  the  intervals. 
Finally,  since  the  waveforms  had  to  be  considered  in  sets,  if  one  waveform 
had  much  more  irregularity  than  others  In  the  set  of  three,  it  largely 
determined  the  choice  of  time  intervals,  the  same  set  of  time  intervals 
being  used,  however,  for  all  three  waveforms  in  each  set.  The  number  of 
time  intervals  used  for  each  set  varied  from  1 to  6 in  our  data  reduction 
procedure  for  Operation  Teapot.  We  believe  that  the  choices  of  intervals 
in  each  case  are  about  the  best  we  can  make:  using  a greater  number  of 

intervals  decreases  the  areas  associated  with  them  and  results  In  greater 
inaccuracy  in  measuring  the  area  with  a planimeter.  (In  general,  we 
planimetered  each  area  twice,  and  in  the  case  of  the  smaller  areas  three 
times  — without  resetting  the  planimeter;  this  improves  the  overall  accur- 
acy somewhat.  Also,  in  many  cases,  including  all  those  in  which  results 
could  be  questioned  for  any  of  several  reasons,  areas  in  question  were 
replanimetered. ) 

B.2  PROCESSING  OF  REDUCED  DATA 

Having  obtained  the  areas  under  the  various  pressure-time  curves  for 
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each  selected  Interval  and  for  the  total  curve  (value  obtained  separately  and 
checked  against  the  sum  over  the  Intervals)  we  proceeded  as  follows  to  obtain 
(4)  values  of  dynamic  pressure  Impulse.  First  we  calculated  average  values 
of  APp  and  APj  for  each  gage  and  for  each  time  Interval  simply  by  dividing  the 
planlmeter  value  APpAt.|,  for  example,  by  Atj,  the  selected  time  interval.  We 
then  read  the  Mach  No.,  M^,  for  the  time  Interval  from  our  set  of  parametric 
straight  line  plots  of  APp  versus  AP^  corresponding  to  Equations  (8).  Inter- 
polation between  the  lines  of  constant  M,  spaced  at  Intervals  of  0.2,  caused 
negligible  Inaccuracy  here  compared  with  other  Inaccuracies  In  the  data  reduc- 
tion process  (alluded  to  previously  and  to  be  discussed  later).  Then  Equation 
(3)  — or  a plot  thereof  — determines  G(M^),  the  Mach  correction  factor  for 
the  1th  time  Interval.  Considering  now  the  APp  and  one  of  the  gages,  ue 
did  this  for  each  value  of  1 and  then  formed  (APp  - APj)^  At^  for  each  i and 
also  E (APp-APj)^  Atj.  Then  we  associated  a weighting  factor  with  each  time 
interval : 

(AP„  - APS)1  It, 

“i  EliPp.iPj),  At, 

♦ 

This  factor  is  the  ratio  of  the  dynamic  pressure  impulse  associated  with  the 
1th  time  Interval  to  the  entire  dynamic  pressure  Impulse  at  the  gage  location. 
(When  the  ground  static  overpressure  gage  Is  used  for  APS,  the  two  gages  are 
a few  feet  apart;  however,  the  ground  gage  does  not  perturb  the  flow  field.) 

Note  that  EW.  = 1 . (10) 

In  using  (9)  we  have  assumed 

APp  - APt,-  M?1.3  (11) 

since,  as  we  have  previously  noted,  the  difference  between  APp  and  APy  is 
negligible  compared  with  unavoidable  inaccuracies  in  data  reduction.  In  fact 
in  the  Operation  Teapot  data  reduction,  use  of  (11)  Is  justified  in  all  cases. 
(Remember  also  that  by  use  of  time  intervals  rather  than  a point-by-point  data 
reduction  process  we  smooth  out  irregularities  in  Mj  — real  or  instrumental 
and  possibly  made  worse  by  data  reduction  procedure  Inaccuracy  — and  substi- 
tute an  average  value  of  for  the  i^  time  interval;  in  a point-by-point 
method  some  higher,  and  lower,  values  of  M would  be  found  in  each  of  our  time 
intervals  in  which  there  are  significant  waveform  irregularities.) 
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With  the  aid  of  Equation  (9)  we  are  now  in  position  to  determine  an 
average  value  of  the  Mach  correction  factor  over  the  entire  pressure- time 
history.  Denoting  this  by  G(M),  we  have 

ETfi)  - EW1  G(M^)  . (12) 

For  reasons  that  will  be  clear  when  we  discuss  accuracy  of  the  method  In 
more  detail  this  Is  the  only  use  we  made  of  the  planlmetered  areas  of  the 
selected  Intervals  for  each  set  of  waveforms;  in  the  remainder  of  the  data 
reduction  process  we  used  only  total  area  under  the  pressure-time  curves. 

Note  that  G(M)  as  given  by  Equation  (12)  Is  greater  than  the  value  which 
would  be  obtained  If  only  one  "Interval"  — the  entire  waveform  — were  used. 
That  Is,  if  we  were  to  base  the  average  Mach  correction  factor  on  the  areas 
under  the  entire  Pp  and  waveforms,  we  would  get  too  low  a value. 

This  Is  because  the  highest  values  of  M and  G(M)  occur  In  Intervals  of 
time  which  have  the  greatest  differences  between  APp  and  APS,  i.e.,  which 
make  the  greatest  contribution  to  dynamic  pressure  Impulse.  Since  the  Mach 
los.  involved  are  fairly  low  for  Operation  Teapot  events,  values  of  G(M)  are 
also  fairly  low  so  that  the  difference  between  G(m]  as  given  by  (12)  and  the. 
value  of  G(Ml  which  would  be  obtained  based  on  the  entire  pressure  pulse  — 
one  Interval  — never  exceeds  a few  percent.  We  now  can  use  the  value  of 
G(M)  obtained  from  Equation  (12)  in  two  ways,  separately  for  the  two  static 
overpressure  gage  results,  each  used  in  conjunction  with  the  single  total 
head  overpressure  gage  result,  to  obtain  4 estimates  of  the  value  of  dynamic 
pressure  impulse.  We  denote  these  estimates  by  Q^,  1 * 1 . . . 4. 

Q1  - /(APp  - APSl)  At/GlSf)  (13) 

where  the  Integral  is  the  difference  of  the  planlmetered  values  of  the  areas 
under  the  entire  curves  of  total  head  overpressure  and  static  overpressure 
(APp  * APy  for  Operation  Teapot  events).  G(M)  is  found  from  Equation  (12). 

Q.j  is  the  value  of  dynamic  pressure  impulse  associated  with  the  static  pres- 
sure gage  placed  at  the  same  location  as  the  total  head  overpressure  gage 
[and  G(M] ) is  also  so  associated]. 

Q2  is  given  by  a similar  equation,  the  only  difference  being  that  the 
planlmetered  value  of  AP$  At  and  the  value  of  G(M)  are  associated  with  the 
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ground  level  static  overpressure  gage;  we  use  the  subscript  2 to  connote  this. 
Then 


Q2  - /(APp  - APs2)  At/Gfi£)  . 

04) 

Next 

Q3  - /q(t)  At/G(M7) 

05) 

and 

Q4  > /q(t)  At/G(M^  . 

(16; 

Here  /q(t)At  Is  the  planlmetered  value  of  the  area  under  the  dynamic  pressure- 
time curve  provided  by  the  previous  (polnt-by-polnt)  data  reduction. 

G(m-j ) and  3(M2 ) rarely  differ  by  more  than  a few  percent  In  this  report, 
so  that  Q3  and  Q4  likewise  rarely  differ  by  more  than  (the  same)  few  percent. 
Since  all  4 estimates  are  dependent  on  a single  gage  result,  and  2 of  the  4 
estimates  are  based  on  one  static  overpressure  gage,  2 on  another,  the  4 esti- 
mates are  certain’y  not  Independent:  however,  as  will  soon  be  clear.  Inaccura- 
cies in  obtaining  dynamic  pressure  impulse  stem  not  only  from  the  raw  data  but 
also  from  the  data  reduction  process  itself. 

Thus  q-j  and  q2  are  both  dependent  upon  the  single  total  head  overpressure 
gage  result  but  are  otherwise  independent.  Q3  and  Q4  are  almost  completely 
dependent  upor.  one  another,  differing  only  because  of  differences  between 
G(M] ) and  G(M2),  but  from  a data  reduction  standpoint  are  almost  independent 
of  Q-j  and  Qj-.  Put  in  another  way,  Q-|  and  Q3  are  only  dependent  upon  one 
another  through  G (M-j ) , which  has  the  same  effect  on  both,  and  Q2  and  Q4  are 
only  dependent  upon  one  another  through  G(M2l  in  carrying  out  the  data  reduc- 
tion process. 

In  many  cases,  after  the  estimates  of  dynamic  pressure  impulse  had  been 

obtained,  one  or  more  of  the  four  estimates  were  either  missing  (gage  failure, 

or  missing  AP$  gage  at  the  ground  range  in  question)  or  obviously  bad.  When 

a reading  was  obviously  bad  a reason  for  this  was  almost  always  apparent. 

Most  usual  reasons  were:  / AP.At  and  /AP«-At  both  considerably  larger  than  the 

P ^ 

differences  between  them,  which  greatly  magnifies  errors  in  the  data  reduc- 
tion process;  baseline  drift  or  instrumental  noise  problems  associated  with 
one  or  more  of  the  gages;  a 'ea  under  curve  too  sinall  and/or  irregular  for 
accurate  planimeter  values  to  be  obtained  especially  in  the  case  of  /q(t)At. 

Rather  than  attempt  any,  necessarily  partly  subjective,  weighting  of  the 
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results  obtained  we,  therefore,  simply  averaged  the  acceptable  Q values  (from 
1 to  4 of  them)  In  each  case.  We  also  made  a comparison  of  results  obtained 
with  values  of  dynamic  pressure  Impulse  Inferred  from  displacement  (based  on 
data  from  Ideal  shots.  Reference  1).  While  these  data  have  considerable 
scatter,  values  of  dynamic  pressure  Impulse  Inferred  from  side-on  and  face-on 
dynamic  pressure  Impulse  versus  displacement  curves  are  almost  always  within 
a factor  of  1.5  of  one  another.  We  did  not  use  any  of  these  results  to  adjust 
any  values  of  Q-f  obtained  from  the  pressure  gage  data.  (See,  however,  Section 
2.2.4  for  a discussion  of  displacement-inferred  values  of  dynamic  pressure 
Impulse;  here  we  have  recourse  to  a larger  data  base  than  In  the  previous  work. 
Reference  1 . ) 

Averaging  the  acceptable  Q values  obtained  gives  a somewhat  higher  weight 
to  Q3  and  Q4  than  to  either  of  Q]  and  Q2;  Q3  and  Q4  correspond  to  the  previous- 
ly reduced  data  (point-by-point)  with  our  Mach  correction  applied.  As  a 
comparative  check  of  the  results  we  plotted  Q2  versus  Q],  Q-|  versus  Q3,  and 
Q2  versus  Q4.  Ideally  the  results  should  then  lie  on  45°  lines.  We  found 
that: 

1)  For  q1  versus  Q3  the  points  are  scattered  uniformly  about  a 45°. 

line; 

2)  For  q2  versus  Q-j  and  Q2  versus  Q4  plots  (with  Q2  as  ordinate) 
the  points  are  scattered  about  a slightly  greater  than  45°  line; 

3)  The  scatter  of  the  points  is  about  the  same  on  all  3 plots. 

This  shows  that  the  values  are  of  equal  reliability.  The  fact  that  Q2  is, 
on  the  average,  a little  higher  than  the  other  values  may  be  due  to  a 

physically  real  difference  resulting  from  the  APp  and  AP$  9age  locations,  or 
it  may  be  associated  with  the  fact  that  the  ground  level  aP$  gage  does  not 
perturb  the  flow  whereas  the  flow  is  perturbed  at  the  other  aP$  gage  location. 
Thus,  we  cannot  say  that  Q2  is  less  reliable  or  less  likely  to  be  correct 
than  the  other  values  of  Q^. 

B.3  ERROR  ANALYSIS 

Finally  we  must  discuss  the  question  of  accuracy  in  more  detail.  We 
start  with  the  planimeter  result.  When  we  measure  the  area  under  the  entire 
pressure-time  curve,  a value  expressed  in  square  inches  is  obtained.  This 
must  be  converted  to  pressure  X time  (psi-msec  in  this  case  although  we  con- 
vert the  final  results  to  kPa-sec).  So  in  addition  to  the  planimetering 
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error,  there  1$  an  error  In  the  scale  length  conversions  on  the  raw  data  plots, 
psi  per  inch  on  the  ordinate  and  msec  per  inch  on  the  abscissa.  We  estimate 
typical  measurements  errors  to  be  about  1?  In  the  time  scale,  2%  in  the 
(smaller)  pressure  scale,  and  ~3£  in  planlmetering  the  small  Irregular  areas, 
in  spite  of  going  around  the  areas  more  than  once  which  diminishes  part  of 
the  planlmetering  error;  (It  doesn't  necessarily  diminish  the  part  of  the 
latter  error  associated  with  Inaccurate  following  of  the  contour,  but  It 
diminishes  any  planimeter  resetting  and  vernier  reading  errors).  In  measuring 
areas  under  the  curve  for  Intervals  within  a pressure-time  waveform  there  is 
a further  error  of  placing  the  vertical  lines  which  subdivide  the  waveform 
into  intervals.  That  Is,  there  Is  an  error  In  At^  which  in  turn  causes  error 
in  APpAtj  and  AP^At^ , as  determined  with  the  planimeter;  since  At.,-  for  the 
APs  gage  waveform  will  not  be  precisely  equal  to  the  measured  Atj  for  the  APp 
waveform  another  small  error  Is  Introduced,  perhaps  ~1X.  We  note  also  that 
the  area  within  an  interval  Is  smaller  than  the  total  area  under  the  waveform, 
usually  by  a factor  of  2 or  3 for  the  most  Important  Intervals  within  a wave- 
form, so  that  the  percent  error  in  planimetering  the  interval  areas  is  larger 
than  for  the  entire  area.  (This  is  the  reason  for  not  using  a larger  number  of 
intervals  than  we  have  used.) 

Thus  before  any  calculations  are  undertaken,  several  small  errors  are- 
introduced  (in  addition  to  Instrumental  errors  In  obtaining  the  raw  data 
with  which  we  begin.)  The  root  mean  square  value  of  the  errors  discussed  is 
typically  -4%  for  the  entire  waveform  and  perhaps  about  5%  for  an  important 
interval  within  a waveform  — for  each  of  APpAt  and  APjAt.  The  absolute  error 
bounds,  however,  can  be  anywhere  from  zero  to  6 or  8%,  for  the  waveform  and 
an  interval,  respectively,  depending  on  whether  the  several  separate  errors 
tend  to  cancel  or  reinforce  one  another. 

The  principal  reason  for  the  much  greater  inaccuracy  in  the  final 
results  achieved  is  due  to  the  magnification  of  this  error  in  the  data  reduc- 
tion process.  Dynamic  pressure  impulse  is  dependent  upon  the  difference 
between  APpAt  and  AP^At  (for  an  interval  or  the  whole  waveform).  With  the 
estimate  of  typical  error  inputs  given  in  the  previous  paragraph  we  can  make 
a tabulation  of  expected  and  (approximate)  upper  bound  errors  as  a function 
of  the  ratio  of  APjAt  to  APpAt.  We  obtain  for  an  entire  waveform: 
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APgAt/APpAt 

0.25 

0.50 

0.75 


% ERROR 

RMS  Approximate  Upper  Bound 


5 

9 

20 


10 

18 

40 


For  an  important  Interval  within  a waveform,  the  entries  in  the  last  two 
columns  should  be  increased  :.y  25  to  40%  (of  the  listed  % error).  For  less 
important  Intervals,  l.e. , those  that  contribute  small  amounts  to  the  total 
dynamic  pressure  Impulse,  still  larger  percent  errors  would  be  obtained.  In 
each  case  the  lower  bound  of  the  percent  error  Is  essentially  zero,  owing  to 
possible  cancellation  of  Input  errors.  Random  addition  or  cancellation  of 
errors  was  undoubtedly  Involved  In  the  elimination  of  obviously  unacceptable 
results  in  some  Instances. 

We  have  also  combined  the  above  analysis  with  Inaccuracy  caused  by  our 
use  of  Pp  * Pj  (for  M < 1.3),  fixing  P0  a 14  psi,  and  interpolating  between 
the  parametric  curves  (straight  lines)  given  by  Equations  (8)  to  determine  M. 

The  inaccuracies  involved  here  are  completely  negligible  compared  with  the 
unavoidable  inaccuracies  just  discussed.  Likewise  errors  in  G(M1-)  and  G(M) 
owing  to  the  assumptions  of  convenience  just  mentioned  are  negligible  compared 
with  the  unavoidable  inaccuracies.  . 

The  errors  we  have  discussed  so  far  relate  to  determination  of  Q-\  and 
Q2*  The  errors  we  introduce  in  determining  Q3  and  Q4  stem  from  two  sources: 

(1)  planimetering  of  the  dynamic  pressure-time  waveforms  (already  provided 
by  a previous,  point-by-point  data  reduction  process.  Reference  (2);  (2)  deter- 
mination of  G(M). 

The  planimetering  error  involved  here  is  larger  than  the  error  in  plani- 
metering APg  and  APp  since  the  waveforms  provided  are  smaller.  However,  there 
is  no  magnification  of  this  error  — this  magnification  has  already  taken 
place  in  forming  APp  - AP$  In  the  previous  data  reduction  from  which  we  have 
q(t)  waveforms.  We  cannot  assess  the  detailed  sources  of  error  involved  in 
obtaining  q(t).  We  roughly  estimate,  however,  that  errors  in  Qj  and  Q4  range 
from  7 to  10%  when  the  area  under  q(t)  curve  is  >0.5  square  inch  to  perhaps 
30%  when  the  area  under  the  curve  is  -0.1  square  inch.  From  our  plots,  pre- 
viously mentioned,  of  Q2  versus  Q-j , Qj  versus  Q4,  and  Qi  versus  Q3,  which 
show  about  equal  point  scatter  in  all  cases,  we  conclude  that  the  inaccuracy 
in  our  method  based  on  planimetering  areas  and  the  inaccuracy  in  the  point- 
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by-point  method  are  roughly  equal. 

The  fact  that  the  errors  are  unavoidably  as  large  as  they  are,  and  are 
still  larger  for  individual  Intervals  within  a waveform.  Is  the  reason  we 
used  the  areas  under  the  total  waveforms,  Mach  corrected  with  G(M),  [with 
G(M)  based  on  Interval  planlmeterlng  as  previously  described]  to  determine 
Q-j  and  Q2  — rather  than  using  the  planlmetered  areas  APpAt-j  and  AP$Atj  along 
with  G(M^)  to  obtain  the  contribution  of  each  interval  to  Qj  and  Q2,  then 
summing  over  the  Intervals  to  obtain  Q-j  and  Q2.  To  show  that  the  method  we 
used  was  the  more  accurate  of  the  two  possible  methods  we  simply  explain 
why  It  Is  that  we  can  determine  G(M)  accurately  from  planimeterlng  values 
within  Intervals,  even  though  the  latter  values  are  themselves  rather  Inaccu- 
rate. The  reasons  for  this  are  as  follows: 

1)  Intervals  with  sufficient  planlmetered  area  to  have  high  weight- 
ing factors,  Wj,  associated  with  them  can  be  planlmetered  much  more  accurately 
than  small  Intervals; 

2)  G(M)  is  only  "large"  for  high  values  of  M and  these  are  always 
associated  with  relatively  large  planimetered  values; 

3)  the  APpAtj  and  APsAt-j  values- will  not  be  nearly  equal  for  any 
interval  with  a large  value  of  W-j  [except  where  M is  small  for  the  entire  wave 
form,  in  which  case  G(M)  Is  not  much  greater  than  unity,  i.e.,  correction  for 
Mach  effect  is  small]; 

4)  our  weighting  factor  is  normalized  so  that  EW-j  * 1 so  that  if 
one  or  more  values  of  W-j  are  too  high,  others  are  necessarily  too  low  thus 
providing  some  cancellation  of  errors;  also  in  the  important  parts  of  a wave- 
form G(M)  varies  less  than  over  the  total  waveform  and  an  error  in  one  mea- 
sured interval  is  most  likely  to  be  compensated  in  an  adjacent  interval,  thus 
augmenting  the  error  cancellation  effect  in  determining  G(M) . 

Had  we  decided  to  use  the  Individual  interval  planlmetered  results,  thus 
determining  q(t^)Atj  for  each  interval,  some  of  the  aforementioned  cancella- 
tion would  also  have  taken  place.  However,  since  use  of  the  individual 
interval  planimetered  areas  provides  an  accurate  determination  of  gTm),  the 
most  accurate  we  can  achieve,  as  previously  explained  — the  method  we  have 
used  yields  as  accurate  a result  as  we  can  obtain  without  need  for  further 
use  of  the  planimetered  results  for  the  individual  intervals. 
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APPENDIX  C 

(This  Section  Is  Unclassified) 

THE  LEAST  SQUARES  PROCEDURES 

C.l  INTRODUCTION 

All  of  the  least  squares  fits  In  this  report  were  obtained  with  the 
TI-59  calculator  using  library  programs  supplied  by  Texas  Instruments. 

All  of  the  fits  were  either  of  the  form 

y 3 mx  + b (1) 

or  y - aQ  + a]x1  + a2x2  (2) 

where  y was  the  natural  logarithm  of  the  desired  quantity;  x was  either  the 
Independent  (controlled)  variable  or  its  logarithm  and  the  least  squares 
variables  Xp  x2  Involved  only  the  special  case  x-j  = x and  x2  3 x . 

The  least  squares  equations  are  obtained  by  minimizing  the  sum  of  the 
squares  of  the  deviations  of  the  data  from  the  curve.  If  the  experimental 
points  have  a variable  scatter  for  a given  small  range  of  x,  weights  should 
be  applied  in  the  procedures  so  that  the  least  squares  equations  correspond- 
ing to  Equation  (1),  for  example,' are  obtained  by  setting  the  derivatives  of 

S 3 EW-jCyj  - (mxi  + b)]2  (3) 

with  respect  to  m and  b equal  to  zero,  thus  obtaining  two  linear  algebraic 
equations  which  can  be  solved  for  m and  b.  Here  S Is  the  weighted  sum  of  the 
squares  of  the  residuals  about  the  fitted  curve  and 

Wi  = a2/^2 

where  a is  a constant  (to  be  determined  from  the  deviations  of  the  data  from 

o 

the  fitted  curve)  and  is  a measure- of  the  expected  deviation  from  the 
true  value  for  an  observation  y^  (of  unit  weight). 

In  the  Texas  Instruments  programs  3 1.  In  individual  cases  to  be 
discussed  we  will  point  out  that  the  programs  used  are  nonetheless  quite 
adequate  for  our  needs.  This  is  mainly  due  to  the  fact  that  in  all  cases, 
we  minimized  the  sum  of  the  squares  of  the  logarithmic  deviations  of  the 
desired  function.  This  is  equivalent  to  minimizing  the  sum  of  the  squares 
of  the  percentage  deviations  from  the  curve  since 

Attn  y)  = f (4) 
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That  is,  for  small  deviations  from  the  curve,  the  fractional  error  is 
equal  to  A (in  y).  For  large  deviations  this  is  not  exact  and  In  our  tabula- 
tion of  results  we  show  both 

Z(£n  y.  - in  y )2  and 
1 1 c 1 yc 

where  yc  ■ the  value  of  y(x^)  corresponding  to  the  curve  fit* 

The  statistical  analog  of  Equation  (4)  (Reference  13)  Is 

Var  (in  y)  * -L  Var  y (Var  * variance)  . (4a) 

y 

Now  In  all  of  the  cases  In  which  we  fitted  the  logarithm  of  the  desired 
function  y,  the  values  of  y covered  a large  range.  If  we  consider,  for 
example,  a displacement  measurement  of  100  metres  to  be  In  error  by  2%,  the 
square  of  its  deviation  from  the  true  value  is  4;  If  we  consider  a displace- 
ment measurement  of  1 metre  to  be  In  error  by  202!,  the  square  of  its  devia- 
tion from  the  true  value  is  0.04.  So  even  though  the  percent  errors  in  our 
data  are  not  necessarily  uniform  (as  best  these  errors  are  known),  assuming 
them  to  be  uniform  is. much  closer  to  reality, than  any  other  assumption  we  can 
make.  From  the  example  just  given  we  readily  see  that  the  coefficients  in  the 
least  squares  fits  would  be  determined  almost  entirely  by  the  data  with  large 
values  of  y if  we  were  to  minimize  absolute  rather  than  percent  deviations 
from  the  fitted  curves. 

C.2  DYNAMIC  PRESSURE  IMPULSE  VERSUS  DISPLACEMENT  (a) 

DISPLACEMENT  VERSUS  DYNAMIC  PRESSURE  IMPULSE  (b) 

In  application  of  our  results  we  need  fits  of  the  data  for  cases  (a) 
and  (b),  i.e.,  with  each  quantity  used  in  the  role  of  independent  (controlled) 
and  of  dependent  variable. 

In  each  subcase  of  this  case  the  data  are  fitted  very  will  by  straight 


lines  on  log-log  plots. 

Thus  for  fitting  on  the  TI-59  calculator  the  fit 

is 

of  the  form 

y = mx  + b 

(1) 

with 

y = in  Iq,  x = in  D 

for  Case 

(a) 

and 

y = in  D,  x 3 in  Iq 

for  Case 

(b) 

* In  most  cases  we  omit 
surrmations  are  taken 

the  subscript  i denoting  the  ith 
over  the  data  points  in  all  cases 

data  point;  the 
so  that  omission 

of 

the  subscript  i will  not  lead  to  any  confusion. 
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and 


b - In  B . 


(5) 


Equation  (1)  can  then  also  be  written  [for  case  (a)],  Iq  - BDm  (la) 
with  a similar  relation  for  Case  (lb).  However,  after  we  have  fitted  case  (a) 
we  can  calculate  (b)  directly  (without  fitting  It)  provided  we  recall  certain 
data  summations  from  the  machine  memory.  For  Equation  (1)  we  obtain  from  the 


TI-59  fit: 


slope  - m * 


£xy  - Exly/N 
Sx2  - (Zx)2/N 


(6) 


v . b . Ey  - ms* 

^Intercept  N 


(7) 


and  the  correlation  coefficient,  r * m (-^jr)  • (8) 

ay 

r2  is  called  the  coefficient  of  determination.  The  various  summations  can  be 
recalled  from  the  machine  memory.  ox2  and  ay2  are  the  variances  of  the  x-array 


and  y-array  data  and  are  given  by  * 

ax2  * 1 E(xi  - x)2  (9) 

a/  = J sty  - y)2  00) 

with  x = ^ Ex  (11 ) 

and  y * ^ Ey  O2) 


N is  the  number  of  data  points. 

Equations  (6)  to  (8)  are  not  symmetrical  in  x and  y and  the  algebraic  in- 
verse of  the  least  squares  fit  Iq  * BDm  (la) 

is  C = (i)^mIn^m  but  this  is  not  a least 

D H 

squares  fit  to  the  data  with  the  role  of  Iq  and  D reversed. 

However,  manipulation  of  Equations  (6)  to  (10)  leads  to  the  equation 


* Bar  over  a quantity  indicates  an  average  value. 
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r 


- J (Ex)2][Ey2 


- ^ Uy)2] 


(8a) 


which  Is  symmetrical  in  x and  y,  i.e. , Is  unchanged  when  x and  y are  Inter- 
changed. 

Thus,  if  we  have  obtained  the  least  squares  fit  Iq  « BDm  we  know  that 
the  correlation  ratio  r has  the  same  value  with  the  roles  of  Iq  and  D 
reversed.  It  can  be  shown  (by  taking  the  origin  at  the  data  centroid),  (Refer- 
ences 13  and  14),  that  the  least  squares  fits 


BD' 


AI, 


m 


(la) 

(lb) 


have  values  of  i and  m which  satisfy  i m * r2  . (13) 

Thus,  having  determined  m (and  r)  In  the  least  squares  fit  to  case  la 
we  obtain  i in  the  least  squares  fit  to  case  (lb)  from  Equation  (13).  The 
constant  A in  Equation  lb  Is  then  determined  by  the  fact  that  both  least 
squares  fits  (la)  and  (lb)  pass  through  the  centroid  of  the  data  distribution. 
(Compare  Reference  15.)  Thus,  from  Equation  (lb)  we  have 


& D * In  A +l  ln  Iq 

‘ (lb) 

and  in  particular 

in  0 ■ In  A +£  In  1^ 

(14) 

where  from  Equations 

(11)  and  (12) 

in  D = Jj-  E(in  D) 
in  Iq  * Jj-  Z(in  Iq) 

(15) 

and  these  values  are  available  from  the  machine  memory  after  running  case 
(la).  Thus  from  Equation  (14)  we  determine  in  A,  hence  A,  and  the  least 
squares  fit  to  case  (lb)  is  fully  determined. 

Having  obtained  the  least  squares  fits,  the  coefficient  of  determination 
provides  a measure  of  the  goodness  of  fit.  By  manipulation  of  Equation  (3), 
with  Wf  =1,  and  Equations  (6)  through  (12)  it  can  be  shown  that 


In  this  report  we  are  concerned  only  with  positive  correlation,  so  that  the 
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nearer  r Is  to  +1,  the  smaller  Is  S,  the  sum  of  the  squares  of  the  deviations 
from  the  fitted  curve,  and  hence  the  better  the  fit  to  the  data.  The  expres- 
sion for  1 - r2  can  also  be  put  In  the  form 


1 


2 . ' »c>2 

r ■ j- 

Hy  - y) 


(17) 


o 

thus  expressing  1-r  as  a ratio  of  the  sum  of  the  squares  of  the  deviations 
from  the  curve  to  the  sum  of  the  squares  of  the  deviations  from  the  mean  (of 
the  data).  The  denominator  may  be  regarded  as  a normalizing  factor.  It  Is 
1i. dependent  of  the  functional  form  used  for  the  fitting  function.  Without 
such  a denominator  1-r  would  tend  to  increase  with  the  number  of  data  points 
— even  If  the  data  were  excellent. 

Now  letting  S^n  j and  S^n  D be  the  sum  of  the  squares  of  the  deviations 
from  the  fitted  curve  when  In  Iq  and  In  D,  in  turn,  play  the  role  of  dependent 
variable,  we  have  from  Equation  (16) 


(1  - 


r2) 


I 


A-  = 


No1 


In  D 

r~ 


In  I, 


No 


In  D 


Here 


JL 

In  Iq 
2 

aln  D 


fcstfn  Iq)2  - 1 (lln  Iq)2] 
l ft(tn  D)2  - J-  ( Un  D)2] 


(18) 


(19) 


in  accord  with  Equations  (9)  through  (12).  Since  we  have  seen  that  r is 
unchanged  when  the  roles  of  D and  Iq  are  reversed.  Equation  (18)  shows  that 
the  sums  of  the  squares  of  the  deviation  from  the  curves  are  not  the  same  for 
the  two  corresponding  least  squares  fits  (to  the  same  data)  but  that 


*ln  D 


2 

°ln  Iq 

~T 

aln  D 


(20) 


In  all  cases  treated  in  this  report 
less  than  S^n 

In  an  application  to  be  made  of  the 
know  the  standard  deviation  of  J ie  value 


we  found 


t to  be  considerably 

xq 


results  of  this  case  we  wish  to 
of  the  slope  of  the  curve  as  well  as 
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the  standard  deviation  of  our  observation.  We  use  the  following  notation 
with  the  subscript  c referring  to  the  curve  fit  In  each  Instance: 


°in  Iq,C 


variance  of  an  observation  of  in  Iq  relative  to  the 
curve;  (note  that  we  have  assumed  the  variance  of 
the  percent  error  to  be  independent  of  the  abscissa 
value) 


» 


variances  of  slopes  of  curve  fits. 


Then  It  can  be  shown  that  (Reference  16) 


0 a . r “ S»  t / N-2* 

I q f C Iq/ 


aln  D,C  Sin  a/"'2 


2 . Sin  d/1'2 

°i  ' „J?  ■ 


No' 


1-r 


2 o' 


in  If 


in  D 
1 in  Ir 


2 , 'in  Iq 


/fl-2 


'm  mT 
Na£n  D 


1-r2  °in  Ir 


Jin  D 


(21) 


(22) 


Using  Equation  (8)  twice,  first  with  D,  Iq  being  Independent  and  dependent 
variable  and  then  with  their  roles  reversed,  we  also  see  that 

°j,  =,  L 
o„  m 
m 


(23) 


In  the  tabular  results  for  the  Displacement  - Dynamic  Pressure  Impulse 
data  curve  fits,  we  list  for  each  subcase  the  quantities  N,  S^n  j^,  o^n  j c, 
S D,  o^n  D c,  i,  m,  a^,  am  In  addition  to  r and  the  fitting  functions. 

To  gain  further  insight  into  the  reliability  of  the  data  fits  we  have 


* Here  there  are  N-2  degrees  of  freedom;  tne  two  degrees  of  freedom  lost 
correspond  to  the  number  of  regression  coefficients  (two). 
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also  tabulated  the  following  additional  quantities  defined  as  follows: 

S’  = Z(^)2 
1 Dc  } 

ln-lar  2 

Sl  = 

lo  Aqc 


(24) 


o,  1/2 

' - |D  A 

aD,C  " 

SIn  1/2 

erf  r = (-^) 

Iq»c  N-2 

and  finally  the  root  mean  square  percent  error  which  is 


(25) 


ED  = (So/N)1/2  x 100% 

Et  = (Si  /N)1/2  x 100% 
*q 


(26) 


The  o'  quantities,  Equation  (25)',  are  standard  deviations  of  the  fractional 
error  (percent  error  apart  from  a factor  of  100).  Thus  when  the  number  of 
data  points  N is  fairly  large  (large  enough  so  that  the  fit  is  fairly  reli- 
able) the  values  of  a^,  aj  are  almost  the  same  as  the  values  of  and  E^, 
respectively,  (aside  from  the  factor  of  100  expressing  the  latter  quantities 
in  percent). 

The  reason  we  have  chosen  to  list  these  additional  quantities,  in  this 
and  in  other  data  fits  to  be  discussed  shortly,  is  that  in  many  cases  the 
data  scatter  is  rather  large  (irrespective  of  the  functional  form  selected 
for  fitting  the  data).  This  means  that  a deviation  A£n  y in  In  y for  a given 
data  point  may  differ  considerably  from  Ay/y.  To  illustrate,  suppose  the 
curve  fit  and  data  point  values  of  a displacement  are  10  metres  and  6 metres. 
The  contribution  of  this  point  to  the  sum  of  the  squares  of  the  deviations 
of  In  D is  then  (In  10-  In  6)2  = 0.261  (irrespective  of  which  value  is  the 
curve  fit  values  and  which  is  the  data  point  value). 

The  contribution  to  Si  is  0 

■ 0.160 
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if  the  curve  fit  value  is  10  and  the  data  point  value  is  6;  if,  however,  the 
curve  fit  value  is  6 and  the  data  point  value  is  10,  the  contribution  to  Sp 

(I^i)2  = 0.444. 

Since  the  least  squares  curve  tends  to  pass  through  the  data  region  with 
a fairly  uniform  distribution  of  points  on  either  side,  on  the  average  the 
curve  fit  will  be  below  the  data  point  about  as  often  as  above  it.  Averaging 
the  above  results  yields  0,160  *-0,4i-  = 0.302  which  is  to  be  compared  with 
the  value  0.261  for  the  logarithmic  deviation.  Had  we  chosen  the  values  10 
and  9.8  instead  of  10  and  6,  the  two  compared  numbers  would  be  virtually 
identical . 

It  is  inherent  in  the  nature  of  the  data  we  are  dealing  with  in  this 

report,  as  is  clear  from  the  discussion  of  the  data,  that  there  is  a fairly 

large  data  scatter  in  many  cases.  Also  in  many  cases  the  number  of  data  points 

is  rather  fewer  than  we  desire.  In  cases  where  the  scatter  of  the  data  is 

small  we  find  that  Sg  is  fairly  close  in  value  to  p and  similarly  is 

fairly  close  to  S£n  j '.  Usually  S£n  D < S^,,  as  in  the  case  illustrated,  and 

S»  • t < S|  ; the  'situation  is  sometimes  reversed,  however,  usually  because 
In  lq  iq 

one  or  two  aata  points  with  large  deviations  from  the  curve  fit  are  below  the 
curve.  Conversely,  for  a case  in  which  one  or  two  data  points  are  far  above 
the  curve  it  can  happen  that  Sy  » y since  there  is  no  limit  to  the  con- 
tribution such  a point  can  make  to  Sy,  y here  indicating  any  quantity  whose 
In  has  been  fitted.  (For  a data  point  below  the  curve  fit  the  maximum  con- 


tribution to  Sy  is 


y-yr 

Hr*)  - ( 


yc  yc 


o-yc  2 
c)  - l 


Finally,  quantities  such  as  Oj  ^ and  Ej  are  directly  associated  with 

the  plotted  data  and  are  easily  visualized.  §ince  they  pertain  to  fractional 

or  percent  errors  they  apply  equally  well  at  all  points  along  the  curve  fit 

(although  the  curve  fit  has  greater  predictive  accuracy  in  the  vicinity  of 

the  data  centroid  than  toward  the  ends  of  the  data  range).  On  the  other 

hand  a^n  j refers  to  the  natural  logarithm  of  the  plotted  data  (ordinate) 

and  the  la?ter  may  be  quite  small  or  large  (or  even  negative),  so  that  while 

a„  T is  a standard  measure  of  the  data  scatter  with  respect  to  the  fitted 
In  Iq  . 

curve,  it  is  not  as  simple  to  interpret  as  are  the  additional  quantities 

which  we  have  tabulated. 
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In  sum,  the  tabulated  measures  of  goodness  of  fit  provide  some  insight 
into  the  assessment  of  the  reliability  and  utility  of  the  data  and  the 
curve  fits  and  help  with  the  interpretation  of  results. 

The  quantities  described  here  are  also  tabulated  for  various  curve  fits 
of  data  other  than  displacement  - dynamic  pressure  impulse  data.  We  discuss 
the  least  squares  treatment  of  the  remaining  data  in  the  following  sections. 

C.3  ADDITIONAL  LEAST  SQUARES  FITS 

Displacement  versus  Scaled  Ground  Range  (a) 

Dynamic  Pressure  Impulse  versus  Scaled  Ground  Range  (b) 

Scaled  Dynamic  Pressure  Impulse  versus  Scaled  Ground  Range  (c) 

In  each  case  the  fitted  curve  is  described  by 


In  y = a0  + a-j  + a2x2  (27) 

where  x = scaled  ground  range.  For  cases  (a)  and  (b)  several  scaling  factors 
were  used  and  an  approximate  optimum  obtained. 

In  running  these  cases  on  the  TI-59  calculator  we  feed  in  the  natural 
logarithm  of  the  y coordinate  of  each  data  point.  Again,  the  machine  gene- 
rates a fit  and  all  of  the  relevant  summations  over  the  data  are  available 
from  the'  machine  memory.  For  this  type  fit  (tri variate) , however,  the 
machine  provides  a quantity  R2  rather  than  r as  a measure  of  the  goodness  of 

fit.  R is  the  multiple  linear  correlation  coefficient  between  y and  the  other 

o o 

least  squares  variables  x-j  = x and  x2  = x ; R^  is  called  the  coefficient  of 
determination  (see  Reference  16)  and  is  given  by 


The  subscripts  1 and  2 refer  to  x-j  and  X2;  r -j  > ry2,  and  r^?  are  given  by 
equations  analagous  to  Equation  (8a). 

r 1 - [xxy  - ^([ix2  - ][sy2  - ®£}\n 


.2...  / „ ,...2x2  o ,..,\2  \1/2 


ry2  - [Zx2y  - Sq^]([zx4  - %^y2  - ^L] 


2x2  v-1/2 


r12  - [XX3  - 5fi](W  - ^][XX4  - 
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(The  superficial  lack  of  symmetry  between  r^  and  results  from  the  fact 
that  the  least  squares  variables  xi  and  as  used  here  are  x and  x2,  respec- 
tively.) The  quantities  r^ , r^*  and  are  the  coefficients  of  correlation 
(also  called  simple  correlation  or  zero-order  coefficients)  between  y and  x-j , 
y and  X2,  and  X]  and  X2»  respectively.  The  important  point  here  is  that 

2 z(y-y  )2 

(1  - R2)  V 07a) 

£(y-yr 

so  that  R2  provides  a measure  of  the  goodness  of  fit  precisely  similar  to  that 
2 2 

provided  by  r in  the  previous  discussion:  R is  a measure  of  the  closeness 

of  fit  of  the  regression  plane  (in  In  y,  xi , X2  space)  to  the  data  points. 

There  is  one  further  change:  replacement  of  N-2  by  N-3  in  Equations  (21) 

and  (25)  so  that 

Gli  Iq,C  = Sln  IqA*3 


Jln  D,  C 


Sln  d/1-3 


(21a) 


and 


°b,C  = 


1/2 


aiQ,c  = (sia/N"3) 


1/2 


(25a) 


We  list  the  same  quantities  in  the  tabular  results  relating  to  the 
fitted  curves  in  this  subsection  as  in  the  previous  subsection. 

Some  additional  fits  were  required  in  the  cases  treated  in  this  subsec- 
tion owing  to  the  form  of  Equation  (27).  This  functional  form  entails  either 
a maximum  or  a minimum  (which  may  fall  within  the  data  range,  near  one  end 
of  the  data  range,  or  at  a point  remote  from  the  range  of  the  data).  In  the 
case  of  Dynamic  Pressure  Impulse  versus  Scaled  Ground  Range  and  Scaled  Dynamic 
Pressure  Impulse  versus  Scaled  Ground  Range,  the  experimental  data  are  such 
that  the  curves  are  concave  upward  and  have  a minimum  usually  approximately 
at  the  large  x (=  scaled  ground  range)  end  of  the  data  range.  In  these 
instances,  in  order  to  have  a fit  of  the  proper  shape  for  large  x,  and  capable 
of  extrapolation  with  no  more  than  the  usual  uncertainty  in  extrapolating  a 
monotonic  well-behaved  function  somewhat  beyond  the  range  validated  b.y  data, 
we  generated  additional  fits  (of  the  same  functional  form  and  with  the  same 
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scaling  In  each  Instance)  over  the  data  contained  in  the  large  x end  of  the 
data  range.  Thus,  in  these  cases  we  have  two  fits  of  the  form  of  Equation 
(27),  one  over  the  entire  data  base  and  one  over  the  portion  of  the  data  base 
relevant  to  large  scaled  range. 

The  need  for  the  additional  fits  stems  from  the  fact  that  there  are  many 
targets  susceptible  to  damage  at  lower  levels  of  dynamic  pressure  impulse  than 
are  the  vehicles  of  interest  in  this  report. 

There  are  possibly  other  methods  of  accomplishing  this  purpose,  includ- 
ing use  of  functions  other  than  provided  by  Equation  (27)  for  fitting  the 
data.  We  note,  however,  that  the  data  itself  is  concave  upward,  so  that  any 
function  conforming  to  this  shape  and  consisting  of  more  than  two  terms  all 
either  increasing  or  decreasing  monotonically  would  seemingly  have  at  least 
one  positive  and  one  negative  term  and  therefore  a possibility  of  a minimum 
somewhere  (though  not  necessarily  near  or  within  the  range  of  the  data).  In 
an  attempt  to  avoid  this  we  used  a one  term  expression  for  the  fit  in  one 
case: 

y = a exp[-b(x-x0)m]  . 

We  chose  a and  x0  so  that  the, curve  was  forced  to  pass  through  a selected 
point:  The  exponent  m allows  some  flexibility  in  adjusting  the  curvature  as 

a function  of  x as  compared  with  a simple  exponential  and  y decreases  mono- 
tonically  as  x increases  for  x > X(j.  (x  = xq  was  chosen  as  the  lowest  value 
of  x in  the  data  range.)  Then  with  a change  of  variable  Y = £n(^-)  and  X = 
x-Xq  the  fitting  equation  takes  the  form  In  Y = In  b + m In  X.  We  fitted 
this  equation  to  the  data  for  one  case.  It  did  not  fit  well  at  all.  (At 
least  part  of  the  reason  for  this,  however,  is  that  the  fitting  equation  is 
ill-conditioned  to  fitting  the  data  without  use  of  a weighting  function  — 
large  percentage  errors  in  y for  small  y have  little  effect  on  the  data;  as 
previously  pointed  out  the  Texas  Instruments  library  programs  for  least  squares 
fitting  do  not  allow  use  of  weighting  functions.) 

Thus,  while  it  is  possible  that  some  other  (probably  not  more  than  two 
term  — however  complicated)  function  or  some  transformation  of  variables  may 
allow  good  fit  to  the  data  without  a minimum  (or  maximum)  in  or  near  the  data 
range  we  believe  it  simpler  and  quite  justifiable  to  use  Equation  (27)  for  all 
of  the  fits  in  this  subsection  in  the  manner  we  have  described.  Equation  (27) 
is  a simple  and  convenient  function  to  use. 
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APPENDIX  D 

(This  Appendix  Is  Unclassified) 

LIST  OF  ABBREVIATIONS,  ACRONYMS,  AND/OR  SYMBOLS 


Ip  - static  pressure  impulse 
Iq  - dynamic  pressure  impulse 

Iqn  - dynamic  pressure  impulse:  initial  estimates  for  ideal/near-ideal 
u events;  first  iteration  results  for  non-ideal  events 

Iq  - final  estimates  of  dynamic  pressure  impulse:  iterated  results  for 

ideal/near-ideal  events;  second  iteration  results  for  non-ideal 
events 

Ps  - peak  overpressure 

D - displacement 

da  - damage 

W,  Y - weapon  yield 

R - ground  range 

X - scaled  ground  range 

HOB  - height  of  burst 

P0  - ambient  pressure 

T0  - ambient  temperature 

sd>  sd»  si»  sp>  st  " scaling  factors 

n - exponent  in  scaling  law 

ao>  a-|,  a2,  b,  l,  m,  xo.  A,  B - constants  in  least  squares  fits 

x,  y - independent  and  dependent  variaoles  in  least  squares  fits 

f»  f]«  f2»  G - functional  forms 

a - standard  deviation 

<j2  - variance 

N - number  of  data  points 

r - correlation  coefficient 

R - multiple  •!  inear  correlation  coefficient 

r^,  R2  - coefficient  of  determination 

S - sum  of  squares  of  deviations  of  data  points  from  fitted  curve 

S'  - sum  of  squares  of  fractional  (relative)  deviations  of  data 
points  from  fitted  curve 


/ S \ 

VN  - number  of  regression  coefficients' 


V2 


= standard  deviation 
of  fractional  error 


E - 


'N  ’ 


x 100  = root  mean  square  percent  error 
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M - mass  of  vehicle 

I - moment  of  Inertia  of  vehicle  (about  designated  axis) 

A - projected  area  of  vehicle 
y - coefficient  of  friction 
V,  Vq  - velocity,  initial  velocity 
F - force 

to  - angular  velocity 

a - angular  acceleration 

T - torque  (vehicle  overturning) 

x - distance  travelled  by  vehicle  (In  equations  of  motion) 
F(V),  g(V)  - functions  of  velocity  (in  equations  of  motion) 
C - constant  in  curve  fits 
K,  K-| , «2  - constants  in  equations  of  motion 


Ps  - absolute  free  stream  static  pressure 
Pj  - absolute  free  stream  total  pressure 
Pp  - absolute  total  head  Pitot  pressure 
PQ  - ambient  pressure 

APs  = ” po  “ free  stream  static  overpressure 

APj  = pT  - p0  — free  stream  total  overpressure 

APp  = Pp  ■ po  — total  head  overpressure 

q - dynamic  pressure 

Qi  - dynamic  pressure  impulse;  subscript  i denotes  ith  method  of 
obtaining  result  (4  values  of  i used) 

y - ratio  of  specific  heat  at  constant  pressure  to  specific  heat  at 
constant  volume  for  air 

M - local  free  stream  Mach  No.  of  flow  behind  blast  front 
Ati  - 1th  time  interval  (in  a pressure-time  waveform) 

G(M)  - Mach  correction  factor 

G(M)  - average  value  of  G(M)  overpressure-time  history 
W-j  - weighting  factor  for  ith  time  interval 
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13  April  1998 


MEMORANDUM  TO  DEFENSE  TECHNICAL  INFORMATION  CENTER 
ATTN:  OCQ/MR  BILL  BUSH 


The  following  reports  have  been  reviewed  by  the  Defense 
Special  Weapons  Agency  Security  Office: 
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DNA-5056F,  AD-C021924,^DTL-80, 0808 
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DNA-582 6F-SUP,  AD-C041417,rf)TL-871167  CfAQ 


The  Security  Office  has  declassified  all  of  the  listed 
reports.  Further,  distribution  statement  "A"  applies  to  all 
the  reports. 


ARDITH  JARRE TT 
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